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FOREWORD 

This report describes part of a comprehensive and continuing' program of re- 
search in multispectral remote sensing of environment from aircraft and satellites. 
The research is being carried out for the NASA Manned Spacecraft Center, Houston, 
Texas, by Willow _ Run Laboratories, a unit of The University of 'Michigan’s Institute 
of Science and Technology. The baste objective of this program is to develop re- 
mote sensing as a practical tool for obtaining extensive information quickly and 
economically.' 

Recently, many new applications for remote sensing have been developed. 

These include agricultural census -taking, detecting of diseased plants, urban land 
studies, mcasurment of water depth, studies of air and water pollution, and general 
assessment of land-use patterns. However, extension of the application of recogni- 
tion techniques to areas other than tho.se used to program the recognition processor 
are-necessary- to make ope rat ional survey systems a reality. The work reported 
herein was directed toward the solution of the problem of extending recognition 
techniques. Four aspects of the work are discussed: (1) simulation of the system- 
atic variations in multispectral data; (2) investigation of techniques for overcoming 
systematic variations during recognition processing of multispectral data; (3) com- 
parison of two different type of likelihood functions; and (4) data collection itself. 

The research covered in this report was performed under NASA Contract NAS 
9-9784, Task C2. The program was directed by R. R. Legaulfc, Associate Director 
of Willow Run Laboratories, and J. D. Erickson, Principal Investigator and Head 
of the Multispectral Analysis Section of the Laboratories. The Willow Run Lab- 
oratories’ number for this report is 31650-74 -T. 
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ABSTRACT 

This report addresses the problem of extending the applicability of recognition- 
processing procedures for multispectral scanner data from areas and conditions used 
for programming the recognition computers to other data from different areas viewed 
under different. measurement conditions. The reflective -spectral region (~0.3 to 
3.0 ^m) is considered. A potential application of such techniques is in conducting 
area surveys. 

Work in three general areas is reported. The first is an examination, from a 
physical basis, of the nature of sources of systematic variation in multispectral 
scanner radiation signals. Calculations with an improved radiative transfer model 
are used to examine the dependence of radiation quantities on a variety of parameters 
and conditions of measurement. The second is an investigation of various techniques 
for overcoming systematic variations in scanner data-. One technique, preprocessing 
with a generalized transformation, is'developed and applied successfully to a set of 
data collected under hazy conditions to remove variations associated with scan 
angle and to extend signatures from one altitude to another. The final study evaluates 
the use of decision rules based upon empirical distributions of scanner signals 
rather than upon the usually assumed multivariate normal (Gaussian) signal distribu- 
tions; It is concluded that the normal assumption is justified and that the time 
required by the increased complexity involved in implementing an empirical rule 
could be better utilized in other ways. 
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INFORMATION EXTRACTION TECHNIQUES FOR MULTISPECTRAL 

SCANNER DATA 

SUMMARY 

~ This report describes part of a comprehensive and continuing program to investigate re- 
mote sensing of the environment from aircraft and satellites. The overall objective of this 
multidisciplinary program at the Willow Run Laboratories of The University of Michigan is to 
develop Earth resources information systems (which use remote sensing) as practical tools to 
provide the planner and decision maker with extensive information quickly andeconomically. 
The major problem discussed in this report is the extension of the application of recognition 
techniques to areas other than that have been used to program the recognition processor. This 
problem must be solved to make operational area- survey systems a reality. Work on four as- 
pects of the problem is reported; the majority of our effort was directed toward the first two. 

Simulation of the systematic variations in multispectral scanner data was carried out by 
means of our substantially improved radiative transfer model. These variations result from: 
■“(1)' changes -in -measurement parameters such as scan angle; (2) the;presence of the atmosphere 
and changes in its properties; and (3) changes in sun position associated with time of day, time 
of year, latitude, and longitude. The model has been developed to the point that close agree- 
ment is achieved between calculated and measured values of sky radiance under light haze con- 
ditions (the only reliable data available that we can use to check the models). We believe that 
we can model the atmospheric effects on radiance 'measurements under cloud- free or uniform 
overcast conditions, and we have developed a reliable basic model upon which we can add re- 
finements (e.g., the effects of scattered clouds). In order to model multispectral scanner data 
more exactly, bidirectional reflectance functions for surface materials should be included in- 
stead of the diffuse conditions that have been assumed to date. 

Techniques were also investigated for overcoming systematic variations during recogni- 
tion processing of multispectral scanner data. A generalized version of earlier preprocessing 
methods was developed and successfully tested on a set of data collected under hazy conditions. 
The procedure was used to eliminate the measurable variation of the data as a function of scan 
angle. Although variations undoubtedly exist after correction, the intrafield and interfield vari- 
ations masked the residual scan angle variations in the data set used for the testing. Bidirec- 
tional reflectance properties of ground covers did not limit the correction to any measurable 
extent because the correction procedure has the inherent ability to remove an average bidirec- 
tional reflectance function. This same correction procedure was used to extend signatures 
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between different data sets, i.e., sets collected at two different altitudes. An a terna P - 
cedure for signature extension shown to he less restrictive involves changing the ^l.hood 
function parameters rather than (or in addition to) preprocessing the data. 

we compared the use of two different typos of likelihood functions in a maximum iikehhood 
decision rule. Decisions were made on data subsets from individual fields With he use o - 
etston rules that assumed muitivariate normal signal distributions and empirical htstogr . 
shmal distributions, respectively. Although the normal likelihood (unction does exactly 
describe data collected by the multispectral scanner, its use does not seriously deg.a e jecogm 
' tion performance for data sets analyzed. Our tests show that removal of systematic effec . 
will result in greater improvements in recognition performance than those that can be o atne 

by use of empirical likelihood functions. « , ^ 

Data were collected for future analysts, development, and testing techniques for extending 
recognition procedures. The flight lines were chosen so that problems associated with arge 

area surveys can be explored. 
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INTRODUCTION AND STATEMENT OF PROBLEM 

Remote sensing with multichannel optical/mechanical scanners is a potential method for 
obtaining area-survey information more rapidly, more accurately, more economically, and/or 
more completely than is possible with methods currently in use. The desired information can 
be derived from recognition maps (or their equivalents) produced by electronic computers de- 
signed or programmed to perform recognition functions on data collected by multispectral scan- 
ners. 

. The essence of airborne and spaeeborne remote sensing is discrimination, on the basis of 
radiation received by the sensor, between various surface materials of interest and their back- 
grounds. Because the radiation emanating from different materials has distinct spectral char- 
acteristics, multispectral remote sensing techniques can make use of observations in more than 
one wavelength interval to discriminate between classes of materials. 

Discrimination over large areas would be much simpler if the spectral radiance received 
from a given class remained constant, but it does not. Throughout the day, the sun changes 
position; the properties of the atmospheric path change with viewing angle and position along the 
flight line; and clouds come and go. Furthermore, the spectral reflectance properties of the 
surface mater mis within any class are subject to variation for physical and/or biological 
reasons. 

The major problem discussed in this report is extension of the application of recognition- 
processing procedures for data from areas and conditions used for programming recognition 
computers (i.e., to set the decision- rule parameters) to other data from different areas viewed 
under different measurement conditions. In other-words, the problem is to develop recognition 
techniques that continue to obtain satisfactory recognition results when conditions change, with- 
out complete reprograming of the recognition processor and without the costly ground observa- 
tions that such reprograming requires. The term signature is frequently used to describe the 
training-set statistics used for recognition, so the procedure of extending recognition-pro- 
cessing techniques is sometimes called signature extension. 

The paragraphs that follow describe, generally, the problem, our approach, and the areas 
in which we have concentrated our efforts. Sections 1.1 and 1.2 are more detailed. The first 
presents and elaborates on our perspective or point of view for recognition studies. The sec- 
ond summarizes the work accomplished during the reporting period. 

We must consider the nature of the variations present in scanner signals in order to de- 
sign methods to overcome their effects. In this report, we restrict ourselves to short- wave- 
length signals in the spectral region of approximately 0.3 to 3.0 jim. We find variations asso- 
ciated with the sun position, the atmosphere, the materials being observed, and the scanner. 
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A multivariate probability density function describes the distribution of incoming signals 
from each class of material, and its shape is broadened by the presence of sources of both 
random and systematic variations. Statistical decision rules are usually employed in recogni- 
tion processing to account for the variations found in the multispectral scanner data; such rules 
are most appropriate for random variations. These decision rules use approximations of the 
signal- probability density function (i.e "-likelihood functions). Therefore, the likelihood func- 
tion is a logical choice to serve* as a focal point for studies of multispectral recognition. 

The analysis and implementation of multispectral recognition techniques aie i elatively 
easy for samples of populations having multivariate normal (Gaussian) density functions. 

Optimal statistical decision rules then can be derived, based on maximum likelihood or related 
principles. For recognition processing of multispectral scanner data, we usually assume that 
the signals have multivariate normal distributions, even though our studies have shown that this 
is not always the case. One possible reason for the non-normality of the data is the presence of 
systematic variations, such as those associated with scan angle (i.e., associated with the length 
of the atmospheric path and the viewing geometry). Experience has shown that, when the scan 
angle is ignored, the probability density functions of different material classes often overlap, 
and poor recognition performance is achieved. 

One way to improve recognition, performance if 1 suc ft .cases _ is^to d efine sey e r al. subc lag gfe g r 
to represent different scan angle intervals of each class of material. This approach can lead 
to many training sets and increased processing requirements. An alternative approach, which 
we have been pursuing, is to preprocess or transform the data to remove the systematic varia- 
tions before applying the recognition rules. In general, both multiplicative and additive adjust- 
ments must be made on the data as continuous functions of scan angle. 

In addition to the scan angle or cross- track dimension of image data, we should also con- 
sider the along-track dimension. Slow changes may occur in the properties of the atmosphere 
along an extended flight path. Even if the atmosphere were to remain constant, changes in the 
position of the sun would result in changes in the received signal characteristics. The approaches 
used to make scan angle corrections in the across-track dimension can be adapted to carry out 
continuous along-track corrections. In addition, one might use feedback from the recognition 
results to modify the preprocessor characteristics to extend the initial training- set data in time 
and space away from the areas used for training. 

Continuous flight, data might not always be available for computation of along-track correc- 
tions. In that case, one must resort to a discrete or discontinuous correction, such as in a flight - 
to-f light or day-to-day processing situation. The generalized preprocessing transformation 
described here provides a mechanism for making corrections to link two data sets. The trans- 
formation can use information about the similarity of material samples observed in both data 

sets or can make use of auxiliary information, or both. 
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A preprocessing operation transforms all signals in exactly the same way, regardless of 
their material class. Therefore, preprocessing is an efficient operation for computation. If 
all materials do not have exactly the same dependence on the parameter controlling the pre- 
processing (e.g., the scan angle), then the correction function represents an average correc- 
tion which minimizes a selected error criterion associated with the sample population used to 
generate the correction function. 

Preprocessing might not be the most effective way to account for the variation associated 
‘wittTa g" iverF parameter', if distinctive information -for various materials is associated with that 
parameter. Alternatively, then, one could adjust the statistical decision- rule parameters in- 
dependently for each material, as a function of the parameter. 

For processing based on the multivariate normal assumption, the mean vectors and co- 
variance matrices of the materials are the decision- rule parameters of interest. The mean 
vector is the more important of the two, and its adjustment has received most of the attention 
to date. In the most general case, adjustments of the covariance matrices should also be con- 
sidered. The relative magnitudes of system noise and radiation signal variance will affect 
whether or not we adjust covariances and the best method to use for the adjustment. 

Insight into the physical processes that result in the observed radiation signals and the sig- 
nals 1 dependence on various factors can be gained through the use of simulation models. The basic 
- model required is one which describes radiative transfer processes within and through Earth’s 
atmosphere. Other models of the sensors, the data processors, and the target materials 
complete the overall system model. The simulation model is particularly useful in explaining 
and characterizing systematic variations in scanner data. 

1.1. PERSPECTIVE - 

In classification or recognition processing of multispectral scanner data, a decision is made 
regarding which of several possible classes of ground surface materials was observed in each 
spatial resolution element. When material class, m, is observed, the primary sensor (the 
multispectral scanner) produces a signal vector, x. Other signal vectors (a, b, and c) from 
auxiliary sources may be introduced to augment the primary signals. The decision process is 
best described with likelihood functions; therefore, likelihood functions are used as the focal 
point for our discussion of all steps in the recognition process. 

At this point, we shall consider only the primary scanner signal vector, x. The likelihood 
function, L(xlm), is mathematically identical to the probability density function, p(xlm),of the 
n-channel signal vector, x = x^, . . . , x^, conditional on material class, m; i.e., L(x(m ) = p(x i m). 
Conventionally, when the probability density function is considered, m is regarded as fixed and 
x is regarded as variable. When the likelihood function is considered, x is regarded as fixed and 
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m is regarded as variable. While acknowledging this distinction, we use the symbol p to pre- 
serve a single notation and use the terms likelihood and conditional density interchangeably. 

1.1.1. SIGNALS AND MEASUREMENT PARAMETERS; MEASUREMENT CONDITIONS 
AND AUXILIARY INFORMATION 

A remote sensing instrument measures the radiation field that reaches it from the scene 
being viewed. The intensity of the measured radiation depends on several measurement pa- 
rameters (wavelength; spatial dimensions, polarization; and time), over which the operator has 
control, and on the prevailing measurement conditions (surface weather; atmosphere constituent 
distribution; sun position; and cloud distribution), over which he has little or ntr control. 

Of necessity, the measuring instrument samples the continuous domain of the measurement 
parameters. For example, every sensor has a specific wavelength passband(s), a restricted 
field of view, a direction of view, a definite polarization characteristic, and a definite interval 
of time over which a measurement is made. One general problem in remote sensing is choice 
of the proper sampling scheme (i.e., designing the sensor and the data collection procedures) 
to allow the desired discrimination between signals from different scene materials to be carried 
out most effectively. 

In multispectral remote sensing, both with cameras and optical/mechanical scanners, 
sampling in the spectral parameter domain has been implemented in discrete wavelength pass- 
bands, i.e., spectral channels. Ideally, these passbands should be chosen to maximize the 
differences between the signals received from various materials on the ground. 

In the spatial domain, the scanner has an optical aperture function and electrical filters 
which limit the spatial frequencies represented by the signals. Furthermore, each scanner 
scene is imaged from a fixed line in space, so that the observation geometry varies continuously 
as a function of scan angle. The direction of scan relative to the sun's azimuth position and to 
the sensor altitude also is an important spatial measurement parameter. 

The polarization parameter usually has been ignored in remote sensing applications, but 
might be used more in the future to minimize the effects of specular reflections or to take ad- 
vantage of certain physical effects. 

Time is a measurement parameter insofar as the starting time, and thus, the initial mea- 
surement conditions, of a collection mission can be selected. However, once a mission has be- 
gun, changes in measurement conditions linked to time are uncontrollable even though they 
might be predictable (e.g., the position of the sun). 

Two conditions of measurement which effect the amount, spectrum, and spatial distribution 
of the radiation that reaches the ground surface are the distribution of atmospheric constituents 
(e.g., haze particles) and the distribution of clouds. The reflecting properties of surface ma> 
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terials depend in part on the surface weather conditions (e.g., temperature, humidity, wind speed 
and direction, and recent precipitation history). 

If information about the conditions of measurement is to be used in recognition processing, 
it must be received through auxiliary information channels. This auxiliary information might 
come from special processing of signals received from the prime sensor, from an auxiliary 
sensor aboard the sensor platform, or from ground-based instrumentation and observations. 

We have identified the likelihood function as the key to an understanding of recognition pro- 
cessing. Likelihood functions may be made simple or complex, according to the number and 
combination of channels considered. These channels may include both the primary signal chan- 
nels (the multispeetral scanner outputs) and the auxiliary signal channels, or their subsets. 

When no auxiliary, information is being used, p(xlm) is the likelihood function. In most rec- 
ognition-processing operations, training sets of data are used to obtain empirical estimates of 
the likelihood function for each class of material. If the conditions which prevailed during the 
collection period were relatively constant, the estimated density functions will be concentrated 
in relatively small regions in signal space (i.e., x space), and decision rules based on them can 
be successfully applied to data observed under similar conditions. The density function esti- 
mated under these relatively constant conditions is not actually p(x|m), the density function over 
alLmeasurement conditions. Rather, it is the density of x, given both m and the specific con- 
ditions of measurement, p(x|m plus the specific measurement conditions). The overall density 
function would be more spread out in signal space and, if applied in recognition processing, 
would result in increased numbers of errors. 

When auxiliary information is used, the definition of the likelihood function is expanded to 
include the additional signal channels. For example, the likelihood function might become 

p(x, a, b, c [ m) 

where x = the signal vector from the multispeetral scanner, a function of wavelength, space, po- 
larization, and time 

a = a signal vector from auxiliary radiation sensors carried on the aircraft (e.g., sun 
sensor) 

b = a vector to describe the measurement conditions, i.e., information from sources other 
than those on the aircraft (e.g., weather observations and ground-based measurements) 
and/or from estimates based on airborne measurements 
c = a vector of measurements (or calculated estimates) of the solar geometry and the 
viewing geometry 
m = the class of material present 
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The expanded signal space now has dimensions equal to the sum of the dimensions of the 
vectors x, a, b, and c. When the data from training sets are taken over the expected variety 
of measurement conditions, the projection or marginal distribution of the cluster of observa- 
tions from any particular class of material onto the x- signal space will be spread out, as illus- 
trated by the x-axis projection of Fig. 1. In the expanded signal space, the cluster of observa- 
tion occupies a much smaller fraction of the total volume, thereby improving chances of dis- 
crimination. For example, joint use of the auxiliary input, a, permits one to narrow the region 
in x space, where the projected observations cluster for any given measurement condition, ag. 

The choice of auxiliary inputs for use in recognition processing will depend upon how easily 
the inputs can be measured and upon how well they represent those measurement conditions 
which produce the greatest variability in the primary scanner signals. The chance of making 
good choices for auxiliary inputs will be improved by development of a better understanding of 
the dependencies of both the multispectral scanner signals and the possible auxiliary signals on 
the variation of the measurement conditions. 

In summary, the problem is to extend the application of recognition- processing rules over 
a variety of measurement parameters and conditions. The empirical use of auxiliary informa- 
tion for this purpose consists of adding information related to the conditions of measurement to 
the signal space and estimating the parameters of the likelihood function. Estimations are based 
on training data drawn from a variety of measurement conditions. Relationships developed 
through theoretical modeling and analyses can be used to augment, structure, or replace the 
empirical procedures. 

1.1.2, DECISION RULES 

Without discussing the details of decision rules, it is sufficient to say that optimum rules 
are rules that depend on the likelihood function of the original or a transformed version of the 
signals. We can implement a practical recognition processor by: (1) choosing a convenient 
mathematical form with which to approximate the true likelihood function for each material 
class; (2) using a training set of the signal data, and corresponding data on measurement con- 
ditions, to estimate the parameters of this mathematical form fob each material class; (3) cal- 
culating the likelihood of each input sample for each material class, using the approximate 
form with its estimated parameters; and (4) deciding to which class the point should be assigned 
based on a criterion function of the likelihoods. A practical processor must accept and process 
signals with statistics described by the true likelihood function, even though the processor classi 
ties signals based on some approximation to the true likelihood function. 

The multivariate normal density function is commonly chosen to approximate the likelihood 
function. The parameters of this form are the mean vector and covariance matrix, and its im- 
plementation in recognition processing requires only the calculation of a quadratic form. In 
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spite of these advantages, the multivariate normal approximation is merely an assumption, the 
validity of which is subject to test. The average probability of mis classification is a convenient 
criterion for deciding whether to use the multivariate normal approximation, some higher- order 
approximation, or, possibly, some simpler form. 

Systematic trends in multispectral data result in large variances and covariances for vari- 
ous classes of materials. In addition, they can cause the multivariate density functions to de- 

— pari-from a normal form. We can remove the effects of systematic trends in the data by 
changing the parameters of the likelihood function according to changes in auxiliary information 

.signals and/or measurement parameters. Alternatively, the average effect of the systematic 
trends on the means of signals can be removed by preprocessing the data, again as a function 
of the auxiliary signals and/or measurement parameters of the prime signal channels. 

1.1.3. PREPROCESSING TRANSFORMATIONS 

Preprocessing has several advantages, even though it is not the most accurate method of 
accounting for systematic variations in signals. First, preprocessing is efficient because only 
one set of corrections is applied to the data. Second, preprocessing separates .the. correction 
process from the decision process. Finally, preprocessing is a proven technique;. its useful- 
ness has been demonstrated for several applications in feasibility studies carried out at The 

- University of Michigan. _ 

There are several reasons for implementing preprocessing transformations. These include: 
(1) decreasing the time of digital recognition computations and/or reducing the complexity of 
analog computing hardware by reducing the number of information channels used in decision cal- 
culations; (2) removing systematic variations from the data and thereby improving recognition; 
and (3) producing variables that are more directly interpretable in terms of the physical pro- 
cesses being observed. In pattern- recognition literature, feature extraction is frequently 
used to describe some of these functions. 

1.1.4. SYSTEM MODELING 

Although the procedures described above might be approached from a purely empirical 
basis, we believe that a balanced program of empirical and theoretical studies provides the 
most reasonable approach to defining future remote sensing systems because the two types 
complement each other and permit more efficient and effective use of developmental resources. 
Such a program would use a model for simulating the entire system from the generation of 
radiation signals through recognition processing. The accuracy required of the model can vary, 
depending on the type of system being studied and other factors. 

An overall system model serves several functions. It provides a mechanism for gaining 
a greater understanding of the physical nature of the problem being studied. For example, it 
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permits one to perform studies to determine the sensitivity of received signals or recognition 
performance to changes in the conditions or parameters of measurement. This type of informa- 
tion can be useful in the design of sensors and processors, and in the evaluation of operational 
system concepts. Finally, the system model can be used for testing and evaluating empirical 
techniques by simulating scanner signals and suggesting modifications of these techniques. 

1.2. SUMMARY OF WORK-P ERF 0RM ET); D U RING THE YEAR 

During the past year of the contract period, progress has been made toward understanding 
and solving the problems of improving multispectral recognition performance and extending 
the applicability of recognition procedures beyond areas used to program (i.e., to set the de- 
cision rule parameters of) the recognition computers. 

Work in four areas is reported herein: (1) simulation of the physical processes that pro- 
duce the radiance signals sensed by a multispectral scanner and examination of the systematic 
variations in these signals; (2) investigation of techniques for overcoming the effects of these 
systematic variations during recognition processing operations, including a generalized de- 
velopment and application of a generalized preprocessing transformation; (3) completion of an 
ihvesfigat'ioh"(begun last year) of the suitability of the multivariate normal likelihood function 
for recognition decision rules; and (4) collection of multispectral data for future analysis and 
testing of processing techniques in an operational area-lsurvey context. " 

1.2.1. SIMULATION AND EXAMINATION OF SYSTEMATIC VARIATIONS IN DATA 

During the past year, we have improved and used a radiative transfer model which char- 
acterizes the spatial, spectral, and temporal distribution of electromagnetic radiation on 
Earth's atmosphere. (The initial development of the model is described in Ref. [1] .) Since we 
are primarily concerned with radiation in the visible and near infrared portion of the electro- 
magnetic spectrum, we have neglected absorption. For tins model, the atmosphere is assumed 
to be plane -parallel, homogeneous, and bounded by a uniform, perfectly diffuse surface. Simi- 
lar models of radiative transfer through planetary atmospheres have been devised by other 
authors but their models usually suffer from one or more of the following restrictions: (1) an 
atmospheric state that is highly idealized and not representative of realistic conditions' (2) a 
radiation field that is determined only at the extremities of an atmosphere; and (3) surface re- 
flectance conditions that are not realistic. 

Our model has removed the first restriction by taking into consideration quantitatively the 
haze content of the atmosphere. We can calculate the spectral transmittance, spectral irra- 
diance, and spectral radiance in any hazy or clear atmosphere for which a relationship between 
optical depth and altitude for each wavelength can be specified. We have eliminated the second 
restriction by finding a simple solution of the radiative transfer equation for radiance within 
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the atmosphere. Thus, our model can be used to calculate irradiance and radiance at any alti- 
tude in the atmosphere. 

In calculations of the radiometric quantities, irradiance, sky radiance, and path radiance, 
one usually assumes a surface albedo independent of wavelength. In our model we consider a 
spectrally dependent surface albedo. Although at the present time, we do not include bidirec- 
tional reflectance in the model, such an effect will be included later. 

One of our major innovations this year has been the inclusion of time as an independent 
parameter. We can now simulate actual flight conditions, taking into consideration the change 
in the solar zenith angle and the position of the sensor over the surface. By analyzing com- 
puter-generated graphs of radiance and irradiance, we can optimize flight conditions for the 
least amount of atmospheric variability, thereby reducing the burden on the processing system. 
Varieties of graphs have been produced and are presented in Section 2. 

The graphs were selected to demonstrate the systematic variations present in scanner data 
and to illustrate the dependence of radiation quantities, such as transmittance, irradiance, total 
radiance, path radiance, and sky radiance, on a variety of parameters and conditions of mea- 
surement. The ultimate goal is to help develop processing techniques that are insensitive to 
systematic variations. The measurement parameters considered include wavelength, scan angle, 
azimuth angle relative to the sun, altitude, time of day, time of year, latitude, and longitude. .. „ 
Measurement conditions were varied primarily by changing the haze content and distribution in 
the atmosphere; various standardized atmospheres, labeled according to their horizontal visual 
range at ground level, were used for most of the calculations. Time was also varied to simulate 
the changes in date that result when data are collected over a period of time, as in an area- survey 
operation. 

A question which has not been answered satisfactorily in remote sensing problems is the 
degree to which one surface element affects a neighboring element. In a haze, there is con- 
siderable scattering of radiation, and analysis of experimental data does seem to show an in- 
fluence of one element on another. A partial solution to this rather difficult mathematical prob- 
lem has been given in Section 2.3.3, in which single scattering in a haze is considered. An 
empirical method of data correction is described in Section 3,1.4. 

The verification of the radiative transfer model depends, of course, on how well it agrees 
with experimental data. Well defined experiments with simple, geometrical and physical con- 
ditions, for which this model is applicable, are not always available, and therefore only partial 
verification of the model is possible. Nevertheless, as the computer- generated plots in this 
report show, the agreement between the model and exact calculations and experimental sky- 
radiance data under clear conditions is very good. Any deviation between model calculations 
and experimental data can be attributed to: (1) approximations used in the development of the 
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model: (2) inability to specify accurately the parameters used in the model; or (3) idealized 
surface conditions. 

A simple radiative transfer model which includes multiple scattering in a hazy atmosphere 
has been developed and tested with a limited amount of actual experimental data. Furthermore, 
the model has been used to simulate flight conditions and to study preprocessing transforma- 
tions. For the most comprehensive model, one should include discrete clouds, absorption, and 
a non- Lambertian surface,- - - - - — 


1.2.2. INVESTIGATION OF TECHNIQUES FOR OVERCOMING SYSTEMATIC VARIATIONS 

Simulation of the type described in the preceding section and examination of real scanner 
data have both shown that there are systematic variations in data that degrade recognition per- 
formance. We have studied methods for overcoming these variations in recognition processing. 

Major emphasis during the year was on generalization of some of the preprocessing trans- 
formations that had been developed and tested earlier. A generalized transformation, called 
the U-V transformation, was developed: it applies both additive and multiplicative corrections 
to scanner data. Section 3.1 describes the transformation and Section 4 presents results of its 
application to real and simulated scanner data." The transformation was applied to a data set 
collected under very hazy conditions, first to correct for scan angle effects and then to extend 
'^decision- rule parameters for data’coUected from a 1000-ft altitude to data collected from 500’0 : “--- 
ft. The scan angle correction is continuous, whereas the altitude correction is discrete or dis- 
continuous. Both can be classed as one -dimensional corrections. Two-dimensional corrections 
are discussed In Section 3.1.3. 

Another type of preprocessing, discussed in Section 3.1.4, is a method for removing the 
effects on the radiance signals of differences between the reflectances of the observed surface 
elements and those adjoining it through coupling by the atmosphere. The method can be re- 
lated to the modeling effort mentioned in the previous section. 

In Section 3.2, we consider, as an alternative to preprocessing, the adjustment of decision- 
rule parameters rather than the transformation of input data. -The alternative approach has 
several potential advantages which are discussed. 

The final presentation in Section 3 is of supporting analysis of interdependencies of radia- 
tion quantities related to sensor outputs. In other words, Section 3.3 presents results of radia- 
tive transfer- model calculations in ways that make clear interdependenceis such as: (1) total 
downward irradiance at any altitude versus that at the surface; and (2) the path radiance as a 
function of time of day versus the corresponding total received radiance. 
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1.2.3. INVESTIGATION INTO THE SUITABILITY OF THE NORMAL LIKELIHOOD FUNC- 
TION FOR DECISION RULES 

Likelihood functions are used in classification decision (i.e., recognition) processes on 
multispectral scanner data. These functions are usually represented by a multivariate normal 
(Gaussian) density function, the statistical parameters of which are determined for the various 
decision classes from subsets of the data. Last year, tests were made of the normality of sub- 
sets of data corresponding to single fields, and all were found to be non- normal. This year, the 
question of how much might be gained by use of a processing rule that does not assume multi- 
variate normality was investigated. A comparison of two maximum likelihood decision rules 
was made on the basis of paired receiver operating characteristic curves, one curve for a de- 
cision rule based on a multivariate normal density function and the other for an empirical 
histogram density function. An operating characteristic curve is a plot of Type I versus Type II 
errors (probability of miss versus false alarm). We conclude that, for the cases considered, 
the improvement in performance resulting from using the histogram likelihood function is not 
sufficient to warrant the added complexity' in processing and the development work necessary 
to implement it. 

1.2.4. DATA- COLLECTION MISSIONS PERFORMED 

One primary and two secondary data- collection missions were flown during the summer of 
' 1 9 T 1 to 'ofefain niulfi spe ctr al 's banner data for ‘future use in the testing and continued development 
of processing techniques to overcome the systematic trends introduced into scanner data by the 
atmosphere, sun position, and scan geometry during area survey operations. The test site is 
located in Ingham County, Michigan. 
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2 

EXAMINATION OF SYSTEMATIC VARIATIONS 
IN MULTISPE CTRAL SCANNER RADIATION SIGNALS BY SIMULATION 

As seen on a global scale, a number of factors can affect the radiation from a given field 
- poiftt-that-is incident upon an airborne- or spacecraft remote sensing device,. Thesejactors are 
as follows for reflected radiation (~0. 3-3.0 jzm): 

(1) Astronomical Factors 

(a) The solar emission spectrum in a given wavelength interval 

(b) Variable sun-earth distance (time of year) 

(c) Variation in solar declination angle (time of year) 

(d) Variation of hour angle of sun (time of day) 

(2) Geographical Factors 

(a) Latitude of field point ” 

(b) Longitude of field point 

(c) Elevation (above sea level) of field point 

(3) Geometrical Factors 

(a) Solar zenith angle 

(b) Solar azimuth angle 

(c) Zenith angle of surface normal 

(d) Azimuth of surface normal relative to solar azimuth - 

(e) Altitude of sensor 

(f) View angles of sensor to surface element 

(4) Physical Factors 

(a) Attenuation and scattering by atmospheric gases 

(1) Temperature versus altitude 

(2) pressure versus altitude 

(3) Relative humidity versus altitude 

(4) Ozone density versus altitude 

(b) Attenuation and scattering by atmospheric aerosol particles 

(1) Aerosol particle density versus altitude 

(2) Aerosol particle size distribution 

(3) Relative humidity versus altitude 
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(c) Cloudiness of sky 

(d) Surface reflectance 

Some of the factors can significantly affect the radiance, whereas other factors have only 
a small influence on the radiance received by a sensor. The overall effect on remote sensing 
is that Earth’s atmosphere can seriously and systematically alter the intrinsic and observed 
radiation of surface objects.- Hence, -we need' to gain an understanding of these phenomena and to 
develop methods to correct the received radiance for the systematic variations caused by the 
atmosphere or account for these variations in recognition processing. A radiative transfer 
model of the atmosphere is an important tool in the process of developing this understanding. 

It can be used to simulate mean radiance signals for multispectral scanners and to study 
parametrically the dependence of radiation quantities on the parameters and conditions of 
measurement. 

2.1. THE RADIATIVE TRANSFER MODEL 

A radiative transfer model has been developed for remote sensing studies. A "detailed 
description of it is presented in a technical report produced last year under this contract [1] 
and in The University of Michigan Engineering Summer Conference Notes [2] . In this section, 
the key aspects of the model are noted. . the , improvements made J.n the-moriel during the past 
year discussed, and a verification of the model is made by comparison with experimental 
measurements. 

2.1.1. ASPECTS OF THE MODEL 

The important aspects of a radiative transfer model, briefly summarized in this section, 
are reviewed in more detail in Appendix I. The reader also may wish to refer to Appendix TV. 

In radiative transfer theory, the most fundamental radiometric quantity is the spectral 
radiance for a particular state of polarization. In a vacuum, spectral radiance is invariant 
with changes in distance. However, when traversing a medium, radiation is attenuated as a 
result of its electromagnetic interactions with the particles that compose the medium. In the 
visible and near-infrared regions of the electromagnetic spectrum, only. (elastic) Rayleigh 
scattering by atoms and molecules and (inelastic) Mie scattering by aerosol particles are 
important as far as remote sensing applications are concerned. 

The scattering by aerosol particles is highly anisotropic and accounts for most of the 
scattering occurring in hazy atmospheres. The calculation of radiative transfer through an 
anisotropic medium is a difficult problem because one must use or approximate the highly 
anisotropic scattering phase function that describes the way that radiation is distributed fol- 
lowing a .single scattering. 
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One quantity of great importance in estimating the scattering or absorbing properties of a 
medium and in interpreting the re.sults of calculations is the optical depth, r, defined as; 

oO 

t = ^K(z)dz 
h 

where h is some base altitude above sea level, and k is the extinction coefficient for the par- 
ticular polydispersion of aerosol particles. 

Another important parameter, visibility or visual range, is frequently used to describe the 
atmospheric state. Following the convention of Kosehmieder [-3] , the visual range, V is defined 
as that horizontal range at which- the visual contrast between a black object and the sky horizon 
reaches 2%.. In our model calculations, V designates particular atmosphere profiles that have 
been identified by Elterman [ 4] . A theoretical, hazeless atmosphere has a visual range of 
336 km. Elterman identifies an atmosphere with V = 23 km as a typical clear atmosphere 
bordering on one with a light haze. The transition between a dense haze and a fog is usually 
taken to_be V = 1.2 km. 

The spectral* radiance, L(r, p, <f>), received by a sensor in an aircraft or spacecraft which 
is viewing Earth’s J surface* is given by the following very simple formula; 

L(r, p, 4>) = L^p, <MT(t, p) + L p (r, p, <t>) ^ 

where r = the optical depth 

p = the cosine of the view angle 
<p = the azimuthal angle 

Lj(p, <P) is the intrinsic spectral radiance at Earth’s surface and is given by; 

2tt i * ' 

Lj(p, </>)=! jp’pTp, 4>, -p f , 0’)L(r o , -P’, 0’)dp’dtf>’ ( 3 ) 

_0 0 

where p*(p, <P, -p’, <t>') is the bidirectional spectral reflectance function of the surface, and 


* All physical quantities described in this section are spectral in nature, but, for simplicity, 
the notation does not show this dependence explicitly. 
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L(r 0 , -p\ 0’) is the total (diffuse plus solar) spectral radiance at the surface. For the special 
case oT a perfectly diffuse (Lambertian) surface, p'({±, 0, -p', 4> r ) is independent of direction, 
and we have 

Lj(4, <t>) =? E_(t 0 ) = Lj (4) 

where p is the hemispherical spectral reflectance or albedo, and E_(tq) is the total downward 
spectral irradiance on the surface. Hence, in this particular case, the intrinsic radiance is 
independent of direction. 

The quantity T(t, p) is called the spectral transmittance and is given by 

T(r, ax) = e’^ T 0' T ^ M * (5) 

L M, 0), the last quantity in Eq.(2), is the spectral path radiance, i.e., that radiance 

r 

which results from the multiple scattering of nontarget radiation into the sensor’s field of view. 

To help understand var-iations in radiance- caused by the atmosphere, we use the radiative 
transfer model to calculate transmittance, irradiance, and total path radiances as functions of ’ 
the parameters on which they- depend?-*, , - . - —■ : ~ ?/• 

2.1.2. MODEL IMPROVEMENTS DURING THE PAST YEAR 

Several significant improvements have been made in the radiative transfer model during 
the past year. The following changes result in increased reliability and accuracy of our model: 

(1) We have included time as an independent parameter in our analysis so that we can 
now calculate spectral irradiance, spectral path radiance, spectral sky radiance, 
and total spectral radiance in terms of local standard time and location on Earth's 
surface. 

(2) Instead of using approximations to the single-scattering phase functions, we now use 
the exact phase functions, as calculated according to Mie scattering theory, and 
realistic particle siae distribution functions, as determined by Deirmendjian [5J. Using 
these general poly disperse phase functions in Eqs.(72) and (73) in Appendix I, we have 
calculated path and sky radiances and compared the results with other computational 
techniques and with experimental data. As we have already seen, the general agree- 
ment is quite good. 

(3) A correction term, which accounts for a more realistic boundary condition for a diffuse 
(Lambertian) surface, was included in the formulas for irradiance and radiance. 
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(4) Some work has been done on a large scale plotting program. One of the advantages of 
this program, aside from its utility in the production of useful plots, is that we can 
now interpolate in altitude, wavelength, and visual range. Hence, we are no longer 
restricted to discrete data points. 

(5) In order to examine the effect which one surface element has on a neighboring element, 
we have had to consider the two-dimensional radiative transfer problem. The complete 

r . - problem of multiple scattering in the two-dimensional situation is quite complicated 

mathematically, but it can be simplified somewhat by consideration of single scattering 
only. In Section 2.3.3, we present the analysis for solving this simplified problem. 

2.1.3. VERIFICATION OF THE MODEL 

The ultimate validity of any model of the physical world lies in its agreement with experi- 
mental data. Unfortunately, there are very few experimental determinations of irradiances and 
path and sky radiances for which a specification of all the boundary conditions is given. How- 
ever, there are some comparisons that we can make. 

" Making use of Mie scattering theory and a knowledge of particle size distributions for 
typical aerosols, Deirmendjian [ 5] has calculated scattering phase functions for several wave- 
lengths and aerosol types. These theoretically exact phase functions were used in the radiative 
transfer equations (Eqs. (60) and (61) of Appendix I) to calculate the path and sky radiances for 
atmospheric conditions under which Ivanov [6] measured sky radiances for six wavelengths 
from 0.40 pm to 0.70 pm at every 15°, in planes coincident with, and at angles of 45° and 90° 
to, the solar plane. The atmosphere was free of clouds and the optical depth of the atmosphere 
also was measured at each wavelength. The surface of Earth was said to be composed of 
herbage; however, our calculations were made for several different values of the surface albedo. 
Figures 2 through 6 illustrate how well our model calculations agree with the experimental 
data. In all cases, the shapes of the theoretical curves agree quite well with those of the 
experimental data. The reflectance (albedo) necessary to cause each pair of curves to match 
at a zenith angle of 0° (arbitrarily selected) was determined. Using these values, we drew a 
derived spectral reflectance curve for this particular surface for comparison with reflectance 
curves for standard surfaces. Figure 7 shows two reflectance curves, one derived in the man- 
ner just described and the other representing the mean of the reflectance curves for green 
vegetation that are present in ERSIS, the Earth Resources Spectral Information System that 
has been developed for NASA under this contract [ 7] . Except for one wavelength, the shapes 
of the curves are similar. The differences in magnitude might be caused by any one of several 
reasons: (1) some atmospheric absorption has taken place and absorption is not included In 
the model; (2) the actual surface is not well described by the mean curve for vegetation; or 
(3) we arbitrarily chose a scan angle of zero degrees for a matching point, and another choice 
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FIGURE 2. DEPENDENCE OF SKY RADIANCE 
ON SCAN ANGLE (IN THE SOLAR PLANE); 
SOLAR ZENITH ANGLE ^65.8°. Visual 
range = 269 km; wavelength = 0.404 /im; surface 
albedo = 45% 


FIGURE 3. DEPENDENCE OF SKY RADIANCE 
ON SCAN ANGLE (PERPENDICULAR TO THE 
.SOLAR PLANE); SOLAR ZENITH ANGLE = 
65.8°. Visual range = 269 km; wavelength = 
0.404 fim; surface albedo = 45% 



FIGURE 4. DEPENDENCE OF SKY RADIANCE 
ON AZIMUTH ANGLE. Visual range = 261 km; 
wavelength - 0.447 /im; surface albedo = 16%; 
solar zenith angle - 65.7°. 
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FIGURE 5. DEPENDENCE OF SKY RADIANCE 
ON SCAN ANGLE (IN THE SOLAR PLANE); 
SOLAR ZENITH ANGLE = 68.6°. Visual 
range = 205 km; wavelength = 0.553 jim; surface 
albedo = 46%. 
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FIGURE 6. DEPENDENCE OF SKY RADIANCE 
ON SCAN ANGLE (PERPENDICULAR TO THE 
SOLAR PLANE); SOLAR ZENITH ANGLE = 
68.6°. Visual range = 205 km; wavelength - 
0.553 jim; surface albedo = 46%. 
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would bring the reflectance curves into closer agreement. On the whole, however, our model 
does show very good agreement with the limited amount of data available, at least as far as 
the angular dependence of sky radiance is concerned. We find good agreement for atmospheres 
with low haze content; additional confirmation is still needed for atmospheres with high haze 
content. Comparisons with multispectral scanner radiance data have been less successful to 
date because the model does not include bidirectional reflectance effects. 

As an additional illustration of the general validity of our atmospheric radiation model, 
we can compare our sky radiances with those computed by Coulson et al. [8] using the results 
of Chandrasekhar's theory [9 j . This comparison is shown in Figs. 8 through 12 (for a Rayleigh 
atmosphere). Figures 8 and 9 illustrate the sky radiance in and perpendicular to the solar plane, 
respectively, for a solar zenith angle of ~37° and Fig. 10 shows the sky radiance for a solar 
zenith angle of 0°. In all cases, the agreement is excellent although there is some deviation 
in the case of an extremely large solar zenith angle of '-84°. 

2.2. RESULTS OF CALCULATIONS 

“in our radiative transfer model, the transmittance, irradiance, and radiance have the follow- 
ing functional dependences: * 

T = T(h, 9, X, V) . ® 

E = E(h, 9 Qt X, V, p) 
l = L(h, e Q , x, v, p, e , 0 ) 

where h is the altitude in km, 6 is the nadir or zenith scan angle, <f> is the azimuthal angle be- 
tween the scan plane and the solar plane, X is the wavelength of the radiation, V is the visual 
range in km, d ^ is the solar zenith angle, and p is the hemispherical reflectance or albedo. 
Implicitly, irradiance and radiance also depend on time of day and year and on location over the 
Earth's surface, i.e., 

(9) 

6 0 = 0 o (t, d, latitude, longitude) 

where t is the local solar time (hours) and d is the day of the year. 

Overcoming the variation of sensor radiation signals with respect to all the parameters 
listed above constitutes the main problem in extending recognition performance. Perhaps the 
greatest variations are the results of changes in scan angle, visual range, and time of day. 

The remainder of this section presents results of calculations made for a variety of values 
for each parameter listed. Additional results are presented in the reprint that is included as 
Appendix IV. For example, the reprint contains both a discussion of contrast transmittance and 
a number of graphs that illustrate its dependence on a variety of parameters. Methods for 
overcoming systematic variations during recognition processing are discussed in Section 3. 
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FIGURE 7. COMPARISON OF REFLECTANCE 
CURVES — FROM MODEL CALCULATIONS 
AND ERSIS DATA 



SCAN ANGLE (°) 


FIGURE 8. DEPENDENCE OF SKY RADIANCE 
ON SCAN ANGLE (IN THE SOLAR PLANE); 
SOLAR ZENITH ANGLE = 36.9°. Visual 
range = 336 km; wavelength ~ 0.546 
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SCAN ANGLE (°) 


FIGURE 9. DEPENDENCE OF SKY RADIANCE 
ON SCAN ANGLE (PERPENDICULAR TO THE 
.SOLAR PLANE ); : SOLAR ZENITH ANGLE = 
36.9°. Visual range = 336 km; wavelength = 
0.546 fim. 



SCAN ANGLE (°) 


FIGURE 10. DEPENDENCE OF SKY RADI- 
ANCE ON SCAN ANGLE (IN THE SOLAR 
PLANE); SOLAR ZENITH ANGLE =0°. Visual 
range - 336 km; wavelength = 0.546 jim. 
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SCAN ANGLE (°) 


FIGURE IX. DEPENDENCE OF SKY RADI- ' 
ANCE ON SCAN ANGLE (IN THE SOLAR 
PLANE); SOLAR ZENITH ANGLE = 84.3°. 
Surface albedo = 2 5%;~ visual range = 336 km; 
wavelength - 0.546 jum. 



SCAN ANGLE (°) 


FIGURE 12. DEPENDENCE OF SKY RADI- 
ANCE ON SCAN ANGLE (PERPENDICULAR 
TO THE SOLAR PLANE); SOLAR ZENITH 
ANGLE = 84.3°. Visual range = 336 km; 
wavelength = 0.546 pim; surface albedo = 25%. 
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2.2.1. VARIATION OF SENSOR SIGNALS WITH SCAN ANGLE 

During data collection, scan angle is the most rapidly changing measurement parameter. 
Substantial variations associated with scan angle are observed frequently in scanner signals. 

The two principal sources of these variations are the atmosphere and the bidirectional reflec- 
tance characteristics of the surface materials. In this section, we present graphs that illustrate 
the extent of scan angle variations associated with the presence of the atmosphere as calculated 
with the radiative transfer model. As seen in Eq. (2), the transmittance and path radiance are 
the two atmospheric parameters of prime interest. 

Figures 13 and 14 illustrate the dependence of transmittance on scan angle and wavelength 
and on scan angle and visual range, respectively. The transmittance decreases as the scan 
angle increases, and the changes are greatest for low visual ranges. 

The variation of path radiance with scan angle and visual range is illustrated in Fig. 15. 
Increasing values for large scan angles and low visual ranges are evident, as is a peak corre- 
sponding to the antisolar angle, which occurs because the scan is in the plane of the sun. This 
asymmetry can complicate the compensation for systematic variations. If the flight path is 
toward or away from the sun, however, the antisolar peak is not observed and the path radiance 
is a symmetric. function of, scan angle. 

Figure 16 shows the variation of total radiance with scan angle and visual range. It is the 
total radiance that is measured by the scanner. Here, the variation for radiance from a diffuse 
surface observed from 1 km is not as great as it was in the case of path radiance, since the 
transmittance changes tend to oppose those in the path radiance. Nevertheless, a distinct anti- 
solar peak still occurs because the scan is in the plane of the sun for this illustration. 

By definition, there can be no dependence of irradiance on scan angle. However, dependence 
of sky radiance on scan angle and visual range is shown in Fig. 17. Note the high peak near the 
solar zenith angle, particularly in the case of very hazy conditions. The sun sensor at the top 
of The University of Michigan's aircraft responds to the weighted integral of sky radiance in 
the hemisphere above the aircraft in addition to the direct solar signal. 

2.2.2. VARIATION OF SENSOR SIGNALS WITH VISUAL RANGE 

Earth's atmosphere has varying degrees of haziness, depending on the density distribution 
of aerosol particles. In our model, the degree of haziness is Characterized by the horizontal 
visual range at sea level. During a long data-collection flight, slow changes can occur in haze 
content and, from day to day, there are usually discrete changes in haze content. The effects 
of such changes must be overcome in processing for large area surveys. 

Figure 18 shows the dependence of spectral transmittance on visual range and wavelength 
for a nadir scan angle of zero degrees. The greatest change occurs when the atmosphere is 
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FIGURE 15. DEPENDENCE OF PATH RADI- 
ANCE ON VISUAL RANGE AND SCAN ANGLE. 
Solar zenith angle = 30°; wavelength = 0.55 4 m; 
altitude = 1 km; azimuthal angle = 0 ° (in the plane 
of the sun); surface albedo for green vegetation. 



FIGURE 16. DEPENDENCE OF TOTAL RADIANCE 
FROM DIFFUSE SURFACE ON VISUAL RANGE AND 
SCAN ANGLE. Wavelength = 0.55 4 m; solar zenith 
angle = 30°; azimuthal angle = 0 ° (in the plane of 
the sun); altitude = 1 km; surface albedo for green 
vegetation. 
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FIGURE 17. DEPENDENCE OF SKY RADIANCE 
ON VISUAL RANGE AND SCAN ANGLE. Solar 
zenith angle = 30°; wavelength = 0.55 (im; altitude = 
0_1 km; azimuthal angle = 0° (in the plane of the sun); 
surface albedo for green vegetation. 
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FIGURE 18. DEPENDENCE OF TRANSMIT- 
TANCE ON VISUAL RANGE AND WAVELENGTH. 
Scan angle - 0°; altitude - 1 km. 
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quite hazy, i.e. for visual ranges from 2 to 8 km. Transmittance does of course change rapidly 
with altitude, especially for very hazy conditions. This effect is illustrated in Fig. 19 for a 
wavelength of 0.55 pm. 

The irradiance also depends on visual range. The dependence of the diffuse upward irradi- 
ance on visual range and altitude is shown in Fig. 20 for a wavelength of 0,55 pm. As in the case 
~ of transmittance, the diffuse upward irradiance varies most rapidly for low altitudes and small 
visual ranges. The change in the diffuse downward irradiance with visual range and altitude 
is even greater, as can be seen in Fig. 21. The spectral shape of the diffuse downward irradi- 
ance is almost independent of visual range (Fig. 22), but the magnitude roughly doubles as the 
visual range is reduced from 23 to 2 km. 

At a wavelength of 0.55 pm, the overwhelming amount of irradiance results from directly 
attenuated solar radiation and, therefore, a plot of total irradiance (direct solar plus diffuse) 
should show little change with altitude and visual range. This is true, as is shown in Fig. 23, 
although one does see that the irradiance at the surface is lower than that above it, and that 
after 3 to 4 km, the irradiance is essentially constant. 

Radiance also changes with visual range. The sky radiance at the Earth's surface can 
change rapidly with visual range and wavelength, as indicated in Fig. 24. The wavelength 
dependence is modified by the spectral albedo of the surface which, in this case, is green vege- 
. .s tation. Table 1 shows the average reflectance for green vegetation for several wavelengths. ■ ^ 
Path radiance is affected in approximately the same way by the surface albedo as shown in 
Fig. 25. Its variation with altitude and visual range is illustrated in Fig. 26 in which case 
there is little change with altitude above 4 km. 

The wavelength/visual-range dependence of the total radiance from a vegetation surface is 
illustrated in Fig. 27 for an altitude of 1 km. For a large visual range, the path radiance is 
small and the total radiance is dominated by the intrinsic radiance from the surface as modified 
by atmospheric attenuation. Thus, we expect and see a spectral dependence similar to that for 
green vegetation at a visual range of 23 km. The spectral dependence is quite different for a 
visual range of 2 km, in which case the path radiance is a major contributor to the total radiance 
at the short wavelengths. The negative correlation between path radiance and transmittance 
again can be seen by comparing Figs. 25 and 27; in particular, the 0.9 pm total radiance de- 
creases with visual range, whereas the path radiance increases. 

2.2.3. VARIATION OF SENSOR SIGNALS WITH TIME OF DAY 

Irradiance and radiance depend on the sun's position which, in turn, depends on the date 
and time of observation and the latitude and longitude. The passage of time is important in area 
survey operations. For simulating signals from an airborne scanner system and sun sensor, 
it is also necessary to specify the flight direction, the altitude, and some estimate of the amount 
of haze present. A typical flight configuration is depicted in Fig. 28. 
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VISUAL RANGE (km) 


.. . - FIGURE 20. DEPENDENCE OF DIFFUSE UPWARD 
FIGURE 19. DEPENDENCE OF TRANSMIT- IRRADIANCE ON VISUAL RANGE AND ALTITUDE. 
TANCE ON VISUAL RANGE AND ALTITUDE. Solar zenith angle = 30°; wavelength = 0.55 /um; 
Scan angles 0°; wavelength = 0.55 jum. surface albedo for green vegetation. 


TABLE 1. AVERAGE REFLECTANCE FOR GREEN 
VEGETATION COMPILED FROM REAL DATA 
FOR SEVERAL WAVELENGTHS 
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FIGURE 21. DEPENDENCE OF DIFFUSE DOWN- 
WARD IRRADIANCE ON VISUAL RANGE AND 
ALTITUDE. Solar zenith angle = 30°; wavelength - 
0.55 jam; surface albedo for green vegetation. 



FIGURE 22. DEPENDENCE OF DIFFUSE DOWN- 
WARD IRRADIANCE ON VISUAL RANGE AND 
WAVELENGTH. Solar zenith angle - 30°; altitude = 
1 km; surface albedo for green vegetation. 
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FIGURE 23. DEPENDENCE OF TOTAL DOWN- 
WARD IRRADIANCE ON VISUAL RANGE AND 
ALTITUDE. Solar zenith angle = 30°; wavelength = 
0.55 |im; surface albedo for green vegetation. 



FIGURE 24. DEPENDENCE OF SKY RADIANCE ON 
VISUAL RANGE AND WAVELENGTH. Solar zenith 
angle = 30°; scan angle = 0°; altitude - 0 km; azi- 
muthal angle = 0° (in the plane of the sun); surface 
albedo for green vegetation. 
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FIGURE 25. DEPENDENCE OF PATH RADIANCE 
ON VISUAL RANGE AND WAVELENGTH. Solar 
zenith angle = 30°; scan angle = 0°; azimuthal 
angle = 0° (in the plane of the sun); altitude = 1 km; 
surface albedo for green vegetation. 



FIGURE 26. DEPENDENCE OF PATH RADI- 
ANCE ON VISUAL RANGE AND ALTITUDE. 
Wavelength = 0.55 i± m; solar zenith angle = 0°; 
scan angle = 0°; azimuthal angle = 0° (in the 
plane of the sun); surface albedo for green 
vegetation. 
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FIGURE 27. DEPENDENCE OF TOTAL 
RADIANCE ON VISUAL RANGE AND WAVE- 
LENGTH. Solar zenith angle = 30°; azimuthal 
angle = 0° (in the plane of the sun); altitude = 

1 km; scan angle = 0°; surface albedo for 
green vegetation. 
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FIGURE 28. TYPICAL FLIGHT CONFIGURATION FOR AIRBORNE SCANNER SYSTEM 

00 - Solar Zenith Angle 
A^ = Solar Azimuthal Angle 
+ 9 = Right Scan Angle 
- 9 = Left Scan Angle 

0 = Relative Azimuthal Angle Between Scan Plane and Solar Plane 
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We can consider a definite place on Earth’s surface. Southeastern Michigan, on a specific 
date, September 1, 1971 and study the variation of trradiance and radiance as a function of time 
of day. Figure 29 shows the diffuse downward irradiance for the conditions stated above. The 
greatest changes occur early in the morning and late in the afternoon, when the visual range is 
small. Figure 30 illustrates the variation of the total downward irradiance. Here, there is 
very little dependence on visual range because: (1) most of the effect is the result of direct 
solar irradiance, which is still the major contributor of total downward irradiance, even at an 
altitude of 1 km; and (2) there is the negative correlation between path radiance and transmit- 
tance . 

If one scans the local North-South meridian continuously, the sun will cross that meridian 
near local noon. Such an effect is clearly seen in the sky radiance plot (Fig. 31) for a moder- 
ately hazy day. 

Now consider a flight due west on a hazy day, at an altitude of 1 km, over a green field in 
Southeastern Michigan, on September 1, 1971. The path radiance is shown in Fig. 32 in terms 
of scan angle and eastern standard time. Path radiance increases rapidly near the horizon, 
especially when the sun is high in the sky. Also, the antisolar peak is clearly noticeable near 
noon. 

Finally, we can look at the total radiance received by a scanner for the conditions stated 
above. In Fig. 33, we see that this radiance has a less rapid change as a function of scan angle 
than does path radiance, although a distinct antisolar peak is still visible. The slower variation 
in this case is the result of the compensating effect of transmittance, i.e. transmittance causes 
the directly attenuated radiance from the surface to decrease with increasing scan angle, where- 
as path radiance increases with scan angle. Significantly, in Fig. 33, the dominant changes of' 
signal are associated with time and the position of the sun. 

2.2.4. VARIATION OF SENSOR SIGNALS WITH SEASON 

The amount of solar radiation reaching a point on the Earth's surface at any time depends 
on the time of year as well as on the time of day. Figure 34 illustrates the variation of the 
diffuse upward irradiance with local solar time in Southeastern- Michigan for the winter solstice 
(December 21), the vernal equinox (March 21), and the summer solstice (June 21). 

Atmospheric scattering decreases the intensity of the direct solar radiation, forming 
diffuse downward radiation as a result. This diffuse downward irradiance is shown in Fig. 35 
as a function of the local solar time for the four seasons. 

The total (solar plus diffuse) irradiance as a function of local solar time is illustrated in 
Fig. 36 for the same dates as in the previous cases. 
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FIGURE 29. DEPENDENCE OF DIFFUSE DOWN- 
WARD IRRADIANCE ON TIME AND VISUAL 
RANGE. Southeastern Michigan, 1 September 1971. 
Wavelength = 0.55 pm; altitude = 1 km; surface 
albedo for green vegetation. 



FIGURE 30. DEPENDENCE OF TOTAL DOWN- 
WARD IRRADIANCE ON TIME AND VISUAL 
RANGE. Southeastern Michigan, 1 September 1971. 
Wavelength -0.55 pm; altitude = 1 km; surface 
albedo for green vegetation. 
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FIGURE 31. DEPENDENCE OF SKY RADIANCE ON 
TIME AND SCAN ANGLE. Southeastern Michigan, 

1 September 1971. Altitude ~ : 0 km; wavelength = 0.55 
4 m; visual range = 8 km; surface albedo for green 
vegetation. 



FIGURE 32. DEPENDENCE OF PATH RADIANCE ON 
TIME AND SCAN ANGLE. Southeastern Michigan, 1 
September 1971. Wavelength = 0.55 4 m; altitude = 1 
km; visual range = 8 km; surface albedo for green 
vegetation. 
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FIGURE 33. DEPENDENCE OF TOTAL RADIANCE 
ON TIME AND SCAN ANGLE. Southeastern Michigan, 
1 September 1971. Wavelength = 0.55 jim; altitude = 1 
km; visual range = 8 km; surface albedo for green 
vegetation. 
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FIGURE 34. VARIATION OF DIFFUSE UPWARD IRRADIANCE WITH TIME AT AN 
ALTITUDE OF 0 km OVER A SURFACE OF GREEN VEGETATION. Southeastern 
Michigan; visual range = 8 km; wavelength =0.55 jam. 



EASTERN STANDARD TIME 


FIGURE 35. VARIATION OF DIFFUSE DOWNWARD IRRADIANCE WITH TIME AT 
AN ALTITUDE OF 0 km OVER A SURFACE OF GREEN VEGETATION. Southeastern 
Michigan; visual range = 8 km; wavelength = 0.55 jim. 
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We can now consider the relative magnitude of direct solar irradiance, downward diffuse 
irradiance, and total downward irradiance as a function of time of day and time of year. These 
effects are illustrated in Figs. 37 and 38. In June, the sun is higher in the sky and all irradi- 
ances will be greater than those in December. On June 21, from sunrise to 8:30 a.m. and from 
5:00 p.m. to sunset, most of the solar energy goes into scattered radiation (diffuse component), 
whereas on December 21 most of the irradiance is diffuse from sunrise to sunset because the 

sun.is always low in the sky. For longer wavelengths or greater visual ranges, the direct 

solar irradiance is larger. 

2.3. OTIIER SOURCES OF VARIATION 

The preceding section presents calculations that demonstrate the systematic variations 
introduced in multispectral scanner signals by changes in the sun’s position, in the viewing 
geometry, and in the condition of the atmosphere. In the real world, one would also expect both 
random variations in the characteristics of the atmosphere and the presence of clouds to intro- 
duce less predictable changes in scanner signals. 

Three other sources of variation have not been discussed yet. These are: (1) the sensor 
itself; (2) the reflectance properties of the" surface being viewed; and (3) the reflectance 
properties of the surface elements adjacent to the element foeins viewed. 

2.3.1. SENSOR CHARACTERISTICS 

No measurement device is noise-free, and the multispectral scanner can introduce random 
and, perhaps, systematic variations into the signals. It is impossible to distinguish between 
random sensor noise and random variations in the radiation signals. However, potentially, 
instability and drift in electrical circuits can produce serious problems for area- survey opera- 
tions. 

2.3.2. REFLECTANCE CHARACTERISTICS OF ELEMENT BEING VIEWED 

The radiative transfer model discussed above incorporates a simple model for the reflect- 
ance properties of the Earth’s surface; it assumes a diffuse (Lambertian) reflector. There is 
evidence in reflectance measurement and modeling efforts that many materials have bidirec- 
tional reflectance characteristics which must be taken into account [ 10] . In other words, the 
amount of radiation reflected depends on both the direction of the illumination and the direction 
of observation. 

In order to see how a bidirectional reflectance property might come about for an agri- 
cultural field, for example, one can argue heuristically that a sensor looking straight down 
sees a maximum amount of soil and a minimum amount of vegetation. At some angle off the 
nadir, there might be no soil visible and hence a different reflectance spectrum. Therefore, 
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EASTERN STANDARD TIME 


FIGURE 36. VARIATION OF TOTAL DOWNWARD IRRADIANCE WITH TIME AT 
AN ALTITUDE OF 0 km OVER A SURFACE OF GREEN VEGETATION. South- 
eastern Michigan; visual range = 8 km; wavelength = 0.55 jim. 



EASTERN STANDARD TIME 


FIGURE 37. VARIATION OF DIRECT SOLAR, DIFFUSE DOWNWARD, AND TOTAL 
IRRADIANCES WITH TIME AT AN ALTITUDE OF 0 km OVER A SURFACE OF GREEN 
VEGETATION. Southeastern Michigan; June 21; visual range = 8 km; wavelength = 

0,55 fim. 

43 




WILLOW RUN LABORATORIES 



EASTERN STANDARD TIME 


FIGURE 38. VARIATION OF DIRECT SOLAR, .DIFFUSE DOWNWARD, AND TOTAL 
IRRADIANCES WITH TIME AT AN ALTITUDE OF 0 kin OVER A SURFACE OF 
GREEN VEGETATION (DECEMBER 21). Southeastern Michigan; visual range = 8 km; 

wavelength - 0.55 jim. 
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the action of scanning a uniform surface, even in the absence of an atmosphere, can produce 
systematic signal variations if the reflector is nondiffuse. 

A method for incorporating a nondiffuse surface reflectance characteristic into the radiative 
transfer model has been identified and should be implemented. While not exact, the approxima- 
tion should allow more realistic simulations of scanner signals. 

Random variations of two types are introduced by the surface materials being viewed. For 
agricultural problems, we can label these within-field variations and within-class, or field-to- 
field, variations. Frequently, these two types of random variations are confused with each other 
and with systematic variations. 

Within a single field, there are random variations in the reflectance characteristics of a 
crop in addition to systematic trends associated with measurement parameters, such as scan 
angle. The statistical decision rules are designed to accommodate this type of random variation. 
Within-class variations, on the other hand, are not as well understood and are more difficult to 
predict and account for in recognition processing. 

2.3,3. REFLECTANCE CHARACTERISTICS OF ELEMENTS ADJOINING THE ELEMENT 
BEING VIEWED 

In last year's report [ 1], it was shown that both the irradiance on an element and the path 
radiance received when the element is observed depend on the reflectances of the elements that 
surround that element. Consequently, the received signal depends in part on the surrounding 
elements. This year, a formulation has been made of a method for analytical evaluation of the 
magnitude of this effect. 

To understand how the radiation field (path radiance) varies in a medium above a surface 
on which the surface albedo changes spatially, it is necessary to consider the mathematically 
complex problem of radiative transfer in a two- or three-dimensional model. The complete 
problem with multiple scattering and an arbitrary single-scattering phase function has not been 
solved, but one can solve the problem if only one scattering is assumed. We can make the 
following simplifying assumptions for a two-dimensional radiative transfer model: 

a. a perfectly diffuse (Lambertian) surface 

b. a monotonically decreasing phase function 

c. homogeneous atmosphere 

d. only single scattering 

e. surface composed of two half spaces, each having a different albedo 

The basic geometry for radiation being emitted from the surface is illustrated in Fig. 39. 

We shall assume the simple case in which the observer is at point P, receiving radiation which 
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originates at Q and is scattered once at P\ Also, assume that the irradiance on the surface 
is known and is independent of a spatial change in albedo and that the radiance from the atmo- 
spheric scattering can easily be calculated. Hence, the once- scattered spectral radiance from 
the right half space is given by: 


T 7t/2 1 


L(r, (X, = | | Jl(t 0 ) e ^ r °" r ^ P(M, <p, <£ r )e 


t A>r^~(r' ~r)/ v d/i i d ^’dT’ 


( 10 ) 


T o° 0 


where, for our particular geometry p = 1, i.e., 6 =0° and L(tq) is the intrinsic radiance from 
the surface. 


We will consider the atmosphere to be composed of a light haze which scatters the radiation 
in a highly anisotropic manner. A mathematically convenient phase function approximating the 
actual phase function is given by; 


|Ae~ a/ ^ ; 0 £i/ si 

p(y) = A (1!) 

0 ; -1 < v £ 0 


where v is the cosine of the scattering angle. We can determine the constants A and a by fitting 
this function to a realistic phase function at v = 1 and by normalization, i.e . : 


p(l) = Ae a = P 0 


and 




l/2jp(i/)d^ = 1 = ~ 

-1 


e~* /u dv=^\ e — dx=^E 2 (a) 


( 12 ) 


(13) 


where E^fa) is an exponential Integral. Thus, from Eqs. (12) and (13), and the recursion formula 
for exponential integrals, we have: 

P 0 f 1 - aE,(a)I = 2 (14) 

can be solved for the constant a. 
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Equation (10) now becomes: 

7 77/2 1 


if f 


L(t, 1 , <f>) = \ j i *->v ' q 
7 0 ° 0 

Hence, the inner integrals (I) have the form: 
tx/2 1 
1 = 1 


L(T„)e' (r 0 “ 7 ’ ,/ ' J, Ae" a ^ l 'c ' <7 ' ~ t) dfi ' d<#>' dr' 




0 0 


where 


O-Tq-T 


0 - a 


From Fig. 39, we see that p’ = cos X cos and 

dp'd ^ 1 = Jdxd^ 

The Jacobian of this coordinate transformation is -cos X- Thus, 


n/2 tt/2 


-■/ /• 


e -k sec X sec ^ cos xdxdt// 


0 o 

where k = ot + {3. Making the transformation x = sec x, transforms the integral into: 


nr /2 


-f 1 
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but, S(kx) = le ^ Sec ^ is Sievert’s integral 
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Tables exist of S(kx) versus kx and therefore the complete integral I can be integrated 
numerically. Thus, 


CO 



Carrying out this integration will allow one to see how rapidly the radiance decreases with 
distance from the edge of the surface. 
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3 

TECHNIQUES FOR OVERCOMING SYSTEMATIC VARIATIONS IN SCANNER DATA 

Several kinds of systematic variations in multispectral scanner data have been discussed, 
and the nature and extent of these variations have been illustrated. There are also random vari- 
ations in the data which have led to the use .of statistical decision rules in recognition processing. 
In this section, we examine possible methods for overcoming systematic variations before or 
during recognition processing and consider the consequences of these methods, on one's ability 
to handle the random variations. The succeeding sections present results of studies carried out 
to test some of these ideas. 

The decision rules that we consider in this section are those which assume that the signals 
have multivariate normal distributions and use the likelihood function in the decision process. 
This is, by far, the most common approach taken toward the recognition processing of multispec- 
tral scanner signals. All one needs for evaluating the likelihood functions are the mean vector 
and covariance matrix for each recognition class considered. 

We have already seen that the received radiance signal consists of two parts, a desired 
part that originates at the ground surface, and an extraneous part that is a result of the atmo- 
sphere's presence between the surface and the sensor. There is a third part that must be con- 
sidered for any practical use of the detected radiation signal. It is the additive noise -equivalent 
radiance that represents noise added by the detector and subsequent electronics and storage 
media. That is: 


L<9, i, j) = p'.. E.< 0)^(0) 4 L pj (0) 4 L nj (21) 

where: L = the radiance signal, 

i = the class of ground cover being scanned 
j = the spectral channel 

*6 ~ a vector that describes the parameters and conditions of the measurement 
p!. = the bidirectional reflectance of the surface (for a diffuse surface, pi. = p../? i, 
where p.. is the diffuse or directional reflectance) 

Ej = the ir radiance in channel j 
Tj = the corresponding atmospheric transmittance 
L . = the path radiance in channel j 

pj 

L^. = the noise- equivalent radiance in channel j 


♦Note that the definition of Q has been generalized in this section to include all parameters 
and conditions of measurement. 
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Equation (21) has exact physical meaning only when the reflectance of the surface is diffuse 
and, consequently, has no angular dependence. In Section 3.1, we discuss a generalization in 
which angular dependence is allowed. 

We shall assume for simplicity that any systematic variations in the data depend on 0 , and 
that there are only two independent sources of randomness in the data, p!j and L^.. The mean 
vector and the elements of the covariance matrix needed for the likelihood calculations are 
respectively: 

<L(«, i, j)> = <p! j >E j (l»)T.(9) + I Jp .(9) + <L nj > (22) 

and 


<[L(0, i, j) - <L (0, i, ])>} [L(fl, i, k) - <L(0, i, k)>] > = 

E j {9ftf0)J!l k {Bn k {OWp' i . - <p! j >] [p^ - <P^>]> + 


<(L . - <L .> )(L . 
nj nj ' nk 


< L 


nk' 


)> 


(23) 


where <> denotes an ensemble average. If we assume that the mean system noise and there- 
fore also the mean noise- equivalent radiance is zero, we have for the mean: 


<L(0, i, j)> - <p^>E^9)Tji$) + 1^(0) (24) 

and for the spectral covariance matrix element: 

CoVjjj ue, i) = x (e)<[p' iy - <p: j >][p^ - <P' a >]> ♦ <L nj L„ k > (25) 

In the next two subsections, we first explore the implications of preprocessing the data to 
remove systematic effects from each channel and then consider the advantages of alternatively 
changing the decision rule parameters, i.e., the mean and, if necessary, the covariance matrix 
Both empirical methods and a method based on the radiative transfer model are considered. 

The final subsection examines the interdependencies of the various radiation quantities to gain 
insight for the other analyses. 

In Section 4, we describe the use of the U(0, j) and V(S, j) method when applied to scan angle 
correction and extension of recognition from one altitude to another. We are faced with the 
question of how to evaluate the performance improvement. Since the desire is to improve rec~ 
ognition performance, the probability of misclassification seems obvious. In fact, we show that, 
for the data set that was tested, an improvement in performance was easily noticeable. There 
are, however, other waysto visualize and evaluate the improvement. For the scan angle effect, 
the use of an incorrect mean vector in recognition processing can easily be seen to cause de- 
graded performance. Thus, a measurement of the dependence of the data mean on scan angle 
both before and after preprocessing, is an easily measureable criterion for data corresponding 
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to one type of ground cover . While any one field might exhibit some residual variations with scan 
angle, the same variations should not occur in the majority of fields. 


3.1. PREPROCESSING 

Equations 24 and 25 show that the mean vector and covariance matrix of signals received 
from a given ground cover depend on the vector, 0, which describes the measurement parameters 
and measurement conditions. However, if the material reflectance functions are independent of 
9 (e.g., are diffuse or Lambertian reflectors), appropriate preprocessing of the data will remove 
the 9 effects from the signal means. When this is done, however, the covariance matrix of the 
transformed data in general will not be independent of a. 


As discussed in Appendix V, appropriate preprocessing includes both an additive and a 
multiplicative correction function. Our generalized preprocessing transformation, called the 
U-V transformation, meets this requirement. 

The operating principle of the U-V transformation is one of adjusting a signal, L(9, i, j ), 
(measured under any arbitrary set of measurement parameters and conditions that are repre- 
sented by $) to _L( i, j), the value it would have had, had it been measured under a reference 
set of conditions and measurement parameters represented by e Mathematically: 

M^O’ ^ j)L(#j i, j) + V(0, j) ^26) 

where U(0, j) and V(0, j) are functions that can be determined empirically from a data set or cal- 
culated with a radiative transfer model and auxiliary measurements and information. Once 
transformed, we process the data, using a fixed, maximum likelihood decision rule. 


The U-V transformation can be applied to more than just diffuse reflectors. For example, 

if not independent of 9 for all materials, but has a systematic variation that is well approxi- 
mated by 

p ij(®> = p ij K lj te) + K 2j( 9 > . (27) 

where p?\ depends only on the material, andK^.(0) and 1^(0) depend only on 9 , the U-V trans- 
formation will remove the 9 dependence from the signal means. Equation (21) now takes the form: 

L(*, i, i)=Ptf 3l (e)+K 4 .{9) (28 ) 

where K gj (0) and K 4j (0) are functions that can be thought of as being apparent, rather than mea- 
surable, (irradiance x transmittance) and path radiance, respectively. This form also handles a 
simple dependence of the actual path radiance on the surface reflectance. To the extent that the 
actual surface reflectance properties depart from the above approximations, the transformation 
applies an average correction factor which minimizes the errors remaining after correction. 
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We now consider the nature of the TJ and V correction functions for the simple diffuse model 
of the radiance signal (Eq. 21). The functions take the forms: 

E 4 (*n)T A8 n ) 


W> j) = 


~EpTfJW 


W, J) = L pj (fl 0 )- U (9, j)L pj <0) 


(29) 

(30) 


When we combine Equations (21), (26),Hf29), and (30), we have 

, M . 


E(9)T(9) n j 


(31) 


We now need to consider the likelihood function for the set of transformed signals L(0^, i, j). 
For the mean, we have 


<U0 Q , i, j)> = <p; j (0)>E.(^ o )T j (^ o ) + Lp.( 9q) 


(32) 


when we make the usual assumption that <L^> = 0. For the elements of the covariance matrix 
- <[L(0 q , i, j> i, j)> ] [L(i9q, i;-k)- -'<L( 0 o Vi, ; 'k)>] > 

= - <p ;j >)<p fe - <•>& > 


< L .L , > (33) 

nj nk v y 


'WV 

E k ( VW 

■ Ej9)Tj8) 

J J 

E k (^T k (S) 


Examination of Eqs. (32) and (33) reveals that, for a Lambertian reflector with p| (0) - 
< p' l j(Q) > j the mean vector is independent of 9, but the covariance matrix is 9- independent as 
noted earlier. 


3.1.1. POTENTIAL USES AND LIMITATIONS OF PREPROCESSING TRANSFORMATION 
Thus far we have described 9 as a measurement parameter or condition of measurement. 
Let us now examine 9 in more detail, by considering specific identities for 9 and corresponding 
potential uses of preprocessing transformations. We can consider two classes of uses for pre- 
processing, depending on whether we consider 9 to be continuous or discontinuous. An example 
of a continuous function is when 9 is the scan angle. The received radiance is scan angle depen- 
dent because the amount of atmosphere between the ground cover and the multispectral scanner 
depends upon the scan angle. Another example is when $ represents the along- track direction 
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(and time). The received radiance can change gradually because of gradual atmospheric changes 
or sun-position movement. Preprocessing can be used to adjust for such changes. 

If we consider 0 to be discontinuous, we again have many potential applications for prepro- 
cessing. We might let L(3 q) represent signals from sunlit portions of the ground and L(8) repre- 
sent cloud-shadowed portions. We would then have a simple processing (or preprocessing) 
scheme to extend signatures obtained from sunlit areas to the shadowed areas. We could let 
L(5 q) and" L(0) represent data taken at two different times at either the same or different alti- 
tudes. Again we would be extending recognition capability without using separate training areas. 
With this particular extension, we would be capable of making area surveys in which one set of 
training data, in principle, would be used for the entire survey. Another possibility is to let 9 
represent an idealized measurement condition, such as a laboratory condition for reflectance 
measurements; then a transformation could be made to convert the measured radiance values 
into reflectance values. 

The potential uses of tills L {9) and L(fl Q ) preprocessing concept are numerous. But, before 
we become too enthusiastic, let us examine some of the limitations. An obvious limitation is 
the assumption that the ground acts as a Lambertian reflector or that all materials have the same 
bidirectional characteristic. Although this limitation is frequently quoted, there are few re- 
liable measurements with which to evaluate the bidirectional reflectance of the various ground 
covers. Another limitation is that the reflectance of a given type of ground cover is not always 
invariant. For example, the reflectance of a corn field can vary from field to field and even with- 
in any one field, and these variations can affect recognition processing results. An effect of 
field-to-field variations is that the training area may not adequately represent the identification 
class. An effect of witlun-field variations becomes apparent when we look at the derivation of 
the functions U(£, j) and V{$, j). These functions must be derived either directly from the data 
or, when the radiative transfer model is used, from measurements of the condition of the atmo- 
sphere and the geometry of the data collection. When the functions are derived from the data, 
intraiield variations can affect the accuracies of measurement and correction. The actual method 
for deriving U(0, j) and V(9, j) from the data will be discussed presently. 

Another limitation is the dependence of radiance upon the albedo (average diffuse reflectance) 
of the ground at and around the point being measured. Measurements made near boundaries of 
fields indicate that the reflectance of one field can influence the measurements of the adjacent 
fields through atmospheric effects (see Section 2.3.3). This is a measurable effect near the 
boundaries, but it becomes negligible at some distance away from the boundary. In Section 3.1.4, 
and Appendix II, we present a model and a suggested algorithm for removal of this effect. 

The noise in the data provides another limitation. There is always Inaccuracy present in 
measurements of the mean vector and covar lance matrix for each training set. If we are 
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developing U{0, } ) and V( ‘J, j ) from the data, measurement errors will introduce additional in- 
accuracy. In addition, the noise statistics may vary under different conditions of 9. 


3.1.2. PREPROCESSING TECHNIQUES FOR ONE- DIMENSIONAL CORRECTION 
In order to use the preprocessing algorithm in which U (9, j) and V{9, j) are used, it is nec- 
essary to find these functions. One method is to consider Eq. (26) to be a function of U(0, j). If 
we have measurements of two materials under conditions 9 and 6^, we will have two simulta- 
neous linear equations in two unknowns. If we label the two materials i = 1, 2, we can find the 
data: 


L(9 1, j) - L(0 2, j) 


(34) 


Similarly, for the function V(0, j) 


L(0, 1, j)L«J 2, j) - L(fl, 2, j)L<9 1, j) 

Vie, j) . TO 7 l,j)- Uf,2. D ~ 


(35) 


Because U(<9, j) and V(0, j) are independent of the material i, these correction functions can be 
used for all data. Obviously, the reflectances of materials 1 and 2 must be different if one is to 
obtain good estimates of U and V. 

A redundant method of calculating V(9 X j) and V(0, j) is better suited to our purpose, ..since, 
our measurements are noisy and the materials chosen to be i = 1, 2 may actually have variations 
for different values of 9. By averaging, we obtain a more representative value for each of the 
correction functions. If we adapt a linear regression method of averaging, we have: 


V(9, j) = 


n 

1 n 

r- 

1 11 

L 

L( V’ L( V k - j) 

ue,i, j) -{Tub, k, j) 

i=l 

k=l 

L k=l J 


n f 1 n 

T ue, t, j) - Tl( 9, k, j) 
i=l L k=l 


(36) 


v ^ " u ( d > i, j) (37) 

i=l i=l 

For twc fields (n = 2), Eqs. (36) and (37) reduce to Eqs. (34) and (35), as expected. Therefore, 
we will only be considering Eqs. (36) and (37). In the above derivation, there was no assumption 
that the two correction functions were continuous functions of 9. Thus, they apply even when 9 
is discontinuous. In fact, the principle application of these formulas is for discontinuous 0. Even 
when 9 is continuous, the calculated functions can vary considerably for small changes of 9 


55 



WILLOW RUN LABORATORIES 


because of measurement noise. Therefore, it may be necessary to combine Eqs. (36) and (37) 
with a smoothing procedure to obtain usable correction functions. 

One method of deriving smoothed functions, U (0, j) and V(9, j), from the data is to assume 
that they are of known form but depend upon unknown parameters. For example, we can assume 
that the correction functions are polynomial functions of 0, and estimate the coefficients. We 
used this method and second- order polynomials where 

. U(0, j) = 1 + ^(0- 9 Q ) +u 2j ((9- a 0 ) 2 r- - - ( 38 ) 

V(0, j) =v 1 .(0 - e Q ) + v 2 _(e - d Q ) 2 - (3 Q ) 

We found the maximum likelihood estimator for the unknown coefficients, which is equivalent to 
minimizing S(j), where: 


n. , 

s(j) j) " u(0k ’ j)L( V ij j) ' V( V j)]2 (40) 

i=lk=l 

where 9, is the k-th sample point of 0. 

K 

An alternative method of applying parameter estimation is to derive the form of the correc- 
. -tion Junctions. from the radiative transfer, model. This method has not been used yet, but is dis- 
cussed in Section 3.1.4. 


3.1.3. EXTENSION OF TECHNIQUES TO TWO DIMENSIONS 

The uses of the preprocessing method that were mentioned in the last section can be con- 
sidered one- dimensional corrections. We discussed letting Q represent either scan angle or 
along- track measurement conditions. If we were to consider both effects at the same time, 
we could first make a scan angle correction, followed by an along- track correction. Another 
method is to find a pair of correction functions dependent upon the two conditions in such a way 
that it is not possible to correct for first one, and then the other, condition. Such a scheme is 
a generalization of the one- dimensional approach. Intuitively, such a scheme makes sense, since 
one would expect the scan angle correction functions to change as the along-track measurement 
conditions change. 

3.1.4. TECHNIQUES FOR REMOVING SYSTEMATIC VARIATIONS IN SCANNER DATA 
WITH THE USE OF THE RADIATIVE TRANSFER MODEL 

The radiative transfer model, with a diffuse surface, has been used to investigate systema- 
tic variations in scanner data caused by atmosphere. For the calculation of total radiance, how- 
ever, it is necessary to include bidirectional reflectance functions for the surface elements in 
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view. Assuming that these have been included, we can develop the formalism to be used in re- 
moving systematic variations caused by atmospheric effects and bidirectional reflectance. 

Let us consider the total spectral radiance, given by: 

L(U * L 0 U)TU) + L U) (41) 

where L Q U ) is the intrinsic radiance at the surface, TU ) is the transmittance from the surface 
to the scanner, and L^(£) is the path radiance. The independent parameter £ is used to i*epre- 
sent a typical independent variable^ Inverting Eq. (41) allows us to solve for theintrinsic ra- 
diance, i.e., 

L o m = F(£)LU) + GU) (42) 

where one can easily see that the F and G functions are given by: 

FUMTUH' 1 . . . (43) 

G«) = -LpU)/TU) (44) 

which are directly calculable from the model. Thus, by performing the (F, G) transformation 
on real scanner data, L(0» we could use the radiative transfer model to correct for atmospheric 
effects. -According to the model, however, we can also write the intrinsic radiance as 

.. (45) - 

where R(£ ) is a known function which depends on the bidirectional reflectance properties of the 
surface element being viewed, and LJO is the final corrected radiance. One should note that 
L c (£) depends on the magnitude of the surface reflectance but not on its bidirectional reflec- 
tance properties. Hence, one can write the corrected radiance as: 


L C U) = [RUH^LqU) = SUJLqU) (46) 

or, combining Eqs. (46) and (42), we get 

L C (5) = P(5)LU) + Q(«) (47) 

where 

_P(i) = FU)SU) (48) 

QU)=GU)SU) ( 49 ) 


The lack of a reliable model for bidirectional reflectances of natural surfaces has contrib- 
uted to our previous omission of bidirectional surface effects. Recently, however, such a 
model has been developed [11 j and we can now use it to determine the P and Q transformation 
functions. Therefore, by using the atmospheric correction functions (F and G), followed by the bi- 
directional correction function, S, one can remove all of the known systematic variations in real 
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scanner data. One also can use functional forms based on these correction functions in empiri- 
cal techniques like those discussed in Section 3.1.2. 

3.1.5 REMOVING THE EFFECTS OF NEIGHBORING MATERIALS FROM RADIANCE 
SIGNALS NEAR BOUNDARIES 

Earlier sections (e.g., Section 2.3.3) have shown that the reflectances of adjoining surface 
elements can affect the received radiance in two ways: (1) they govern the amount of path ra- 
diance to an extent; and (2) they affect themagnitud e of irradiance on the surface and, there- 
fore, the amount of radiance received directly from the element. Insofar as the reflectances of 
the adjoining elements are different from those of the observed element, we have developed a 
method for removing the effect of the differences in a preprocessing operation. The result of 
applying the preprocessing transformation would be to make each point appear as it would have 
had it been surrounded b 5 ' material of its own reflectance; the transformation would not remove 
the effect of the atmosphere entirely. 

Appendix II presents a derivation of the correction formula. A block diagram that illustrates 
its application is presented here in Fig. 40. The function, K(o>), depends on a function, h{x, y), 
which is measured from the data and describes the spatial spreading effect that the atmospheric 
haze has on the radiation. When it is completed, the modeling effort described in Section 2.3.3 
should provide a functional form for h(x, : r - , ^ .. .. u 

3.2. ADJUSTMENT OF DECISION- RULE PARAMETERS 

The preceding section discusses methods for preprocessing data to remove systematic 
variations prior to recognition processing. There are several reasons why it might be prefer- 
able to change the decision rule parameters instead (i.e., to change the mean vector and, perhaps, 
the covariance matrix of each class of material): 

(1) If the covariance matrix is independent of the controlling measurement parameter, 0, 
changes in the mean vector do not change the variances as they do in preprocessing. 

(2) The correction functions can be different for each material rather than being fitted to 
a few (or averaged over all) materials. 

(3) Corrections can be made at relatively slow rates corresponding to the scan rate of the 
scanner or slower rather than at the video rate (i.e., the rate at which resolution ele- 
ments are scanned). 

In computing the likelihood for an individual data point, the stored mean value for the class 
of interest is subtracted from the input data point and the difference is used in the computation. 
Since the operation of changing the mean value is a linear one, it has no effect on the variance 
of the signal. 
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In a case analogous to preprocessing, the mean value in channel j for a given material class, 
i, can be written in a new form: 

<U0» J)> = A(0, j)<L( 0^, i, j)> + B(0, j) (50) 

where A (0, j) and B(0, j) are correction functions independent of material class, and where 
L(0 q, i, j) represents the radiance under a standard or arbitrarily chosen reference condition 
of measurement, In general, however, the correction functions A and B could depend on the 
material class as well as on the channel. The use of a correction function of the form of Eq.(50) 
would permit a different pair of correction functions for each class of material. 

Data rates and bandwidth requirements are important aspects in the design of processing 
systems. The bandwidth is particularly important for analog processing systems. Slow rates of 
change of decision- rule parameters are ideally suited for implementation in hybrid (digital/ 
analog) processors. 

Just as discussed earlier for the correction functions U and V, the functions A and B might 
be determined empirically from the data, might depend on auxiliary inputs and theoretical rela- 
tionships, or might combine the two methods. We also can introduce a method which can be used 
to improve either of the others; it is the concept, of feedback or adaptive, .processing. With this 
concept, each time a datum-point is classified,, that point is included in the training set for- the- .j 
particular decision outcome and will slightly alter the likelihood function for subsequent deci- 
sions. For example, if we let 9 represent the along-track direction of data collection, the 
mean vector and covariance matrix associated with each decision outcome would gradually 
change as the illumination conditions and atmosphere changed. 

3.3. SUPPORTING ANALYSIS OF INTERDEPENDENCIES OF RADIATION QUANTITIES 

Both the preprocessing techniques (Section 3.1) and the decision- rule- adjustment techniques 
(Section 3.2) depend on the interrelationships between the various physical radiation quantities and 
their relationships with primary and auxiliary sensor signals. The radiative transfer model 
provides a mechanism for simulating the radiation quantities and analyzing these interdepen- 
dencies. 

In Fig. 41, we consider the ratio of the total downward irradiance at an altitude of 5 km to the 

total downward irradiance at the surface versus wavelength for three visual ranges. As expected, 

there is a greater change for the V = 23-km case than for the V = 2-km case because, in the 

-4 

former, we approach pure Rayleigh scattering which has a \ wavelength dependence. In a 
dense haze, the wavelength dependence is much weaker. 

Figure 42 illustrates the variation of the ratio of diffuse downward irradiance to total down- 
ward irradiance with respect to visual range for three altitudes. The haze present at low 
altitudes scatters much of the direct solar radiance and, hence, a large amount of solar 
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FIGURE 40. BLOCK DIAGRAM OF CORRECTION TO REMOVE 
THE EFFECT OF NEIGHBORING MATERIALS FROM RADIANCE 
SIGNALS NEAR BOUNDARIES 



FIGURE 41. RATIO OF TOTAL DOWNWARD IRRADIANCE AT AN ALTITUDE OF 
5 km TO TOTAL DOWNWARD IRRADIANCE AT EARTH'S SURFACE AS A FUNC- 
TION OF WAVELENGTH. Solar zenith angle = 30°; surface albedo = 0.5. 
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energy is redistributed in the form of downward diffuse radiation. For higher altitudes, above 
the majority of the haze, only a relatively small amount of radiation exists in the diffuse state 
and therefore the ratio is independent of visual range. 

We can now simulate the response of a sensor which is measuring the total downward irradi- 
ance at some altitude as a function of the time of day. Such an effect is clearly seen in Fig. 43 
for three visual ranges. For large solar zenith angles there is considerable attenuation in the 
dix'ect solar radiation, and since the attenuation is less at higher altitudes, the ratio of total ra- 
diance at some altitude to its value at the surface is always greater than unity. From Eq. (70) 
one can see that tins ratio is given by: 


E Jr) T 0 " T 
E_(t q )~ 


(51) 


for large solar zenith angles, 0^. 

If a sensor mounted on top of an aircraft were responsive only to diffuse sky radiation, we 
could compare the irradiance at some altitude to its value at the surface. This comparison is 
shown inFig. 44. At the Earth’s surface, the two irradiances are of course equal and we then get 
a 45° line. At the top of the atmosphere, by definition, there is no downward diffuse component. 
For intermediate altitudes, the slope of- the curve simply reflects the attenuated solar spectrum. " 

Now consider the variation of spectral radiance. The equation for total radiance is: 

H = L 0 T + L p = h, + L p (52) 

where Lq is the intrinsic radiance from a surface, T is the transmittance, L is the path ra- 
diance, and = LqT is called the bare radiance. Also, transmittance and path radiance are 
negatively correlated in relation to the scan angle. Figure 45 illustrates the ratio of path ra- 
diance to total radiance, plotted as a function of time of day. As the solar zenith angle increases, 
the intrinsic radiance decreases and, therefore, the ratio of to L t increases. For very hazy 
conditions, the path radiance increases and the transmittance decreases, thereby increasing the 
ratio. The slight maximum near noon is the result of strong backscattering of radiation by the 
haze layer. The important point of Fig. 45 is that the ratio of path radiance to total radiance is 
essentially constant during a major part of the day for the downlooking case illustrated there. 

There exist simplified experimental methods for determining the total transmittance of 
Earth's atmosphere from the ground. Therefore, it is interesting to compare the path radiance 
and transmittance in terms of visual range and wavelength. This dependence is shown in Fig. 46, 
in which we see that a relatively small change in transmittance (~10%) corresponds to a large 
change (^100%) in path radiance, at least for the spectral region 0.50 ^m-0.70^m. Path radi- 
ance is a monotonically decreasing function of transmittance for a constant wavelength. 
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FIGURE 42. RATIO OF DIFFUSE DOWNWARD IR RADIANCE TO THE TOTAL 
DOWNWARD IRRADIANCE AS A FUNCTION OF VISUAL RANGE. Southeastern 
Michigan; solar zenith angle = 30°; wavelength ~ 0.55 jam; surface albedo = 0.5. 



EASTERN STANDARD TIME 


FIGURE 43. RATIO OF TOTAL DOWNWARD IRRADIANCE AT AN ALTITUDE OF 
5 km TO TOTAL DOWNWARD IRRADIANCE AT EARTH'S SURFACE AS A FUNC- 
TION OF TIME. Southeastern Michigan; wavelength = 0.55 fim; surface albedo = 0.5. 
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FIGURE 44. DEPENDENCE OF DIFFUSE DOWNWARD IRRADIANCE 
AT ALTITUDE h ON DIFFUSE DOWNWARD IRRADIANCE AT 
EARTH'S SURFACE. Solar zenith angle = 30°; visual range = 8 km; 
surface albedo =0.5. 



O 
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FIGURE 45. RATIO OF PATH RADIANCE TO TOTAL RADIANCE AT 
AN ALTITUDE OF 5 km AS A FUNCTION OF TIME. Wavelength = 0.55 
Jim; nadir scan angle = 0°; surface albedo = 0.5. 
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Finally, we can compare (Fig. 47) the path radiance to the attenuated intrinsic radiance 
(bare radiaiice) for various wavelengths and visual ranges. From 0.40 |iim to 0.50 /rm, a rel- 
atively large change occurs in bare radiance as a result of the increase of the solar ir radiance 
with increasing wavelength and because transmittance increases with wavelength. Beyond 0.50 
jam, the solar spectrum decreases and transmittance increases, leading to an almost constant 
value for bare radiance, whereas path radiance decreases rapidly with increasing wavelength. 
As in Fig. 46, the path radiance in Fig. 47 decreases monotonically with bare radiance for a 
given wavelength. 
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FIGURE 46. DEPENDENCE OF PATH RADIANCE AT AN ALTITUDE OF 5 km ON 
TOTAL ATMOSPHERIC TRANSMITTANCE. Solar zenith angle = 30°; nadir scan 

angle - 0°; surface albedo =0.5. 



FIGURE 47. DEPENDENCE OF PATH RADIANCE ON BARE RADIANCE AT AN 
ALTITUDE OF 5 km. Solar zenith angle = 30°; nadir scan angle = 0°; surface 

albedo = 0.5. 
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4 

APPLICATION OF PREPROCESSING TECHNIQUES 

A description of preprocessing techniques is presented in Section 3. In this section, we 
apply the U-V transformation to both simulated and real scanner data. 'For purposes of these ^ 
tests” the parameter vector, Q , is used to represent the sensor scan angle in most instances, 
and altitude in one instance. 


4.1. RESULTS OF USING THE U-V TRANSFORMATION ON SIMULATED DATA 

The radiative transfer model was used to generate simulated data to make an initial test on 
the U-V correction procedure. Simulated data can be made noiseless, so only a small amount 
of data need be analyzed in order to obtain a qualitative performance evaluation. 

Data were generated for the following conditions: 


Visual Range (V): 
Surface Reflectance (p): 
Wavelength (X): 

Altitude (h): 

Scan Angle: 

Sun Position (°): 


2, 8, 23 km 
0.1, 0.2, 0.3, 0.4 
0.8, 0.6, 0.4 pm 
500, 1000, 5000 ft- ■ 

-40° to 40°, every 10° 

30° from overhead, in plane of scan 


Figure 48 shows the calculated radiance values for V = 2 km, X = 0.4 pm, altitude = 1000 ft, 
and reflectance values of 0.1, 0.2, 0.3 and 0.4. The solid curves represent the radiance values 
before U-V correction. A considerable amount of scan angle effect can be seen in these curves. 
For example, the radiance received at a scan angle of -30° (the antisolar angle) from a 0.1 
reflector exceeds the radiance from a 0.3 reflector at +30°. 

The dashed curves show the U-V corrected radiance curves as a function of scan angle. 

The U and V correction functions (Eqs. 34 and 35) were derived with the data for p = 0.1 and 
0.4 for each angle and were then applied to all the data. Within the accuracy of the calculations, 
the corrected radiance values are independent of scan angle, as expected, because the present 
radiative transfer model does not include bidirectional reflectance effects. Figures 49 and 50 
show the U and V curves, respectively, for the three different visual ranges. Figure 49 shows 
very little multiplicative correction used on the data, the maximum being less than 1.04. Al- 
most all of the correction comes from the additive correction shown in Fig. 50. In Section 4.2, 
we will find that more multiplicative correction is used when measured, rather than computed, 
data are used. 




66 



WILLOW RUN LABORATORIES 



< "40 -30 -20 -10 0 10 20 30 40 


SCAN ANGLE (°) 

FIGURE 48. CALCULATED 0.4-/xm RADIANCE VALUES 
AS A FUNCTION OF SCAN ANGLE 



-40 -30 -20 -10 0 10 20 30 40 


SCAN ANGLE (°) 

FIGURE 49, U FOR CALCULATED RADIANCE VALUES 


67 



WILLOW RUN LABORATORIES 


We used aquadratic U- V correction (using equations derived from Eq. 40) to reduce the scan an- 
gle effect in Fig. 48. The results are shown in Fig. 51 in which a residual scan angle effect is seen 
The correction was designed to have the corrected radiance, for any reflectivity, be that mea- 
sured at a scan angle of 0°. The deviation of the curves from straight lines is relatively large, 
which is predictable since, in looking at Fig. 48, one would not expect quadratic additive and 
multiplicative corrections to be effective. 

In Fig. 52, we show data corresponding to the parameters used for Fig. 48, except that the 
wavelength is 0.8 pm rather than 0.4 pm. At this wavelength, very little scan angle effect is 
present. 

As one might expect, the scan angle effect became more noticeable with decreasing wave- 
length and visibility and increasing altitude. More explicitly, the following chart lists those 
parameter combinations that showed radiance values at one scan angle and one reflectance to be 
greater than values computed at a different scan angle for a lower value of reflectance: 


X 

Visibility 

Altitude 

(pm) 

(km) 

(ft) 

0.4 

2 

1000 

0.4 

2 

5000 

0.6 

2 

5000 

~o:4' 

‘ "8“ --• 

’’ 5000 

0.4 

23 

5000 


Note that the first parameter combination corresponds to Fig. 48. 

It is premature to use the above list to define conditions under which the scan angle effect 
is severe. For instance, the radiation transfer model does not include bidirectional reflectance, 
but rather assumes Lambertian reflectance. One could expect that atmospheric phenomena alone 
might cause a severe scan angle problem for those parameter combinations in the table and that 
other combinations of parameters might be added when the model is expanded to include bi- 
directional reflectance, 

4.2. RESULTS OF THE USE OF THE U-V TRANSFORMATION ON REAL SCANNER DATA 
Our choice of multispectral scanner data to test the U-V concept was based on several 
considerations. First, we wanted more than one data set gathered over the same ground area 
to enable us to try signature extension from one set of data to another. Thus, we could compute 
the mean and covariance signatures for training areas in one data set, and then find a U-V trans- 
formation which can convert the signatures so they can be used for processing a second data set. 
One of the data sets should be collected from a fairly low altitude (where fields are large in 
angular subtense) in order to use Eqs. (36) and (37). To use these equations for scan angle 
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SCAN ANGLE (°) 

FIGURE 50, V FOR CALCULATED RADIANCE VALUES 



FIGURE 51. CALCULATED RADIANCE VALUES AFTER 
APPLICATION OF QUADRATIC CORRECTION 
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correction, measurements on the same ground covers at both scan angles 0 and 0^ were neces- 
sary. Data taken on a relatively hazy day were desired so that a measurable scan angle effect 
would be present that could be at least partially caused by the atmosphere. We also wanted 
data for which ground truth was available, to test not only the ability to remove the scan angle 
effect but also the change In ensuing recognition performance. There is not a direct relationship 
between the presence of the scan angle effect and the ability to recognize correctly, because 
recognition is limited by other factors, e.g., interfield and intrafield reflectance variations. 

Another consideration in the choice of data was that the data correspond to an agricultural 
area in which the fields were easily delineated and each field was homogeneous (more or less). 
Recent data were desired because continuing improvements in our data collection equipment 
make the more recent data more representative of the present-day state of the art. To be 
confident of the datas' quality we wanted data collected close in time to one of the periodic cali- 
brations. We wanted the day of the flight to be cloudfree, or almost so, since we were not im- 
mediately concerned with the problem of recognizing both shadowed and sunlit ground covers 
in the same data set. Finally, to reduce time and expense, we wanted data already digitized, if 
possible. 

The data set that was chosen seems to have fulfilled our requirements. The data were col- 
lected on 3 September -1969, over a rural area near the Willow Run laboratories. The visibility 
was estimated to have been 4 miles, and the flight log noted that heavy haze was present. A few 
small clouds were present at high altitude. From the photographs, the ground cover we studied 
appears to be free of cloud shadows, although it is possible that shadows might explain some of 
the intrafield variations that were noticed. Two data sets had already been digitized, taken at 
1000 and 5000 ft, respectively. Data from 2000 and 10,000 ft also are available, but not digitized. 

A total of 29 fields were chosen for study. From these fields, 7 were chosen for U-V cal- 
culations and 14 for training for recognition. The 29 fields chosen had recognizable boundaries 
and represented multiple types of ground cover. Each of the 7 fields used for the U-V determin- 
ations encompassed at least half of the total angular coverage. The area overflown was divided 
in approximately the middle of the scan by Willow Road, so three fields on one side of the road 
could be used to find the U-V correction; the other four on the other side could also be used, 
with Eqs. (36) and (37). 

The 29 fields along with the ground truth are listed in Table 2. A map of the fields appears 
in Fig. 12 in [11] . 

The 1000-ft data were processed to compute the U-V functions. Filtering of data reduced 
the number of scan angles and scan lines, each by a factor of four. Figures 53 and 54 show the 
computed U and V functions, respectively, for channel 1 (0.40-0.44 j lm). There is considerable 
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TABLE 2, FIELDS SELECTED FOR TESTING 
U-V TRANSFORMATION 


Field Type 

Number 

Fields Used for 
Determining 
U and V 
Functions 

Fields Used for 
Recognition 
Training 

Corn 

4 

X 

x 


10 




17-18 




40 

X 

X 


52 


X 


66 


- 

Soybean 

67 


X 


3 


X 


15 


X 


20-21 


X 


39 




49 

r 


Grain Stubbie 

55 




9 

X 



23 

X 

X 


25 


X 


42 




53 



Idle 

56 




5 




44 


X 


54 

X 

X 


11 



Bare Soil 

24 

38 

X 

X 


57 



Alfalfa 

16 


X 

Alfalfa Mix 

26 

27 

X 

X 
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variation in the two figures, caused at least partially by the noise in the data and by rapid local 
changes in reflectance of the fields. Some of the changes are caused by the scanning, since the 
instantaneous field of view of the scanner is small enough to produce aliasing at 1000 ft. There- 
fore, we used Eqs. (38) and (39) (the quadratic formulation) to compute the U-V corrections. 

In Fig. 55, we have a plot of the scan angle effect as seen in the means of six soybean 
fields for scanner channel 1 (0.40-0.44 pm). In each field, the along-track averages of the data 
for the 10 channels was computed, and these averages were plotted as a function of scan angle. 

It appears that there is a general trend followed by all of the fields, as well as some field varia- 
tions. In fact, the variation over any one field with scan angle exceeds the field-to-field varia- 
tion at constant scan-angle. 

In Fig. 56, the same fields are shown after the U-V corrections were applied. The isolated 
points are believed to be caused by noise bursts in the data introduced during the A-D process. 
The fields are shown to exhibit much less scan angle effect than before correction. In these 
data, the field-to-field variations are larger than the scan angle effect observed in any one field, 
which is just the opposite of the effect observed in the uncorrected data. Generally, a desirable 
goal for scan angle correction would be to make the scan angle effect somewhat less than the 
field-to-field variations. Usually, refinements in the correction beyond this criterion will not 
be reflected in the recognition results. 

In Figs. 57 and 58, we see for six corn fields the same kind of plots as Figs. 55 and 56. 

The trends observed in the soybean fields are clearly repeated. Our plots show the scan angle 
effect directly, so it is evident that the U-V correction does reduce this effect. A much less 
sensitive method of testing for the removal of the effect is to use recognition results. The intra- 
field and interfield variations tend to obscure the improvement gained by the scan angle cor- 
rection. The U-V correction would not remove all of the limitations to perfect recognition. 

Still, it is desirable to see the effects of the U-V correction on recognition, not as a measure 
of U-V effectiveness, but rather as an indication of how well data can be recognized before and 
after corrections. This knowledge would be useful in deciding whether correction should be 
applied and in predicting the resultant improvement in recognition capability. 

Recognition results were obtained for both the 1000-ft data and the 5000-ft data for the 29 
fields. For the 1000-ft data, an area located near the center of the scan was selected to derive 
the likelihood functions used in the maximum likelihood decision rule, in each of the 14 training 
fields. Recognition processing was performed; no decision was made for any point for which 
the maximum likelihood did not exceed a level calculated to reject one point in 10 (assuming 
normal data). 

This recognition procedure was repeated for the 1000-ft data set in three conditions: 

1) as collected; 2) after filtering; and 3) after both filtering and U-V corrections were applied. 
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SCAN ANGLE (°) 

FIGURE 54. V CALCULATED FROM MEASURED RADIANCE 



FIGURE 55. SCAN ANGLE EFFECT IN SIX SOYBEAN FIELDS 
(UNCORRECTED) 




I 


WILLOW RUN LABORATORIES — | 

i 

S 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 


FRACTION OF SCAN ANGLE COVERAGE ANALYZED 


FIGURE 56. SCAN ANGLE EFFECT IN SIX SOYBEAN 
FIELDS AFTER QUADRATIC U-V CORRECTION 
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FIGURE 571 SCAN ANGLE EFFECT IN SIX CORN FIELDS - 
{UNCORRECTED) 


0.9 


> 


w 

o 


w 

£ 


o 

co 


0.8 


0.7 


Legend 

A 4 


□ 10 
® 37-18 
* 52 
° 66 
- 67 


qPo °oa 






4 /v ^ro „ 

> x *°0 0 0 * x, ‘ 1 

* o 


ao 

G 

A A 


*>A * * A ° 

A AA** a* a ° * A A 

a a a ** 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 

FRACTION OF SCAN ANGLE COVERAGE ANALYZED 


FIGURE 58. SCAN ANGLE EFFECT IN SIX CORN FIELDS 
AFTER QUADRATIC U-V CORRECTION 
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Different decision-rule parameters were used for each of the data conditions, each derived 
from the data set for which the rules would be applied. In addition, the 5000-ft data were pro- 
cessed as collected and after U-V corrections were applied. For the second condition, the 
likelihood decision-rule parameters were obtained from the 1000-ft data, after it had been 
corrected (condition 3). The 7 fields that were used for deriving the scan angle correction 
were used to correct the 1000-ft rules for use on the 5000-ft data. The procedure was to use 
Eqs. (36) and (37) on the mean vectors of the likelihood rules, with and 9 denoting the 5000-ft 
and 1000-ft data, respectively. Tn this manner, we are directly testing the U-V concept for 
signature extension between data sets. The covariance matrices in the likelihood functions were 
used without change. Probably a preferable method would be to change the matrices in accord- 
ance with the gain changes in the scanner amplifiers. Another alternative would be to use 
measured variances in the channels at a constant scan angle. 

A summary of the recognition results is shown in Table 3. The first three columns are for 
the 1000-ft data and the last two for the 5000-ft data. The first three columns correspond to 
data that are: 1) uncorrected; 2) filtered; and 3) filtered and U-V corrected, respectively - 
The final two columns are for 5000-ft data that have been either: 1) scan angle corrected and 
recognized with the use of the 1000-ft likelihood decision functions with corrected mean vectors; 

, T pi;. 2) noncpErec,ted r andTecognized with training sets from 5000-ft data. To obtain the first row, 
we simply counted recognitions in each field. For the second row, a field was assumed to be 
composed of the material with the largest number of recognitions. The third row is the average 
probabilities of no-decisions in the 29 fields, caused by the maximum likelihood function being 
less than the critical value. These numbers reflect the variance changes caused by filtering 
and the variation of the data with scan angle; they also reflect data from material types not 
represented in the training sets. Also, the training areas were all taken near the center of the 
scan. The fifth row is similar to the first, except that only the decisions were used, so that all 
of the no-decision points were ignored. The last row is similar to the first, although only the 
14 training fields were examined. Thus, this row excludes the interfield variations which are 
confusion factors as far as U-V scan angle corrections are concerned. 

As we mentioned earlier, several phenomena in addition to the scan angle effect cause dif- 
ferences in the table. However, it can be seen that the U-V scan angle correction does indeed 
improve recognition capability over the unfiltered condition. We can also see that, to a limited 
extent, signature extension between data sets at two different altitudes has been accomplished. 
This is best seen by comparison of the last two columns of Table 3. 

Tn order to obtain the results indicated, the fields identified as grain stubble, idle, and bare 
soil were grouped together, as were alfalfa and alfalfa mix. This was done because the recogni- 
tion showed that there was a considerable amount of confusion within the two groups, possibly 
caused by the inaccuracy or incompleteness in the ground information. 
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We also used the mean value of the data from one field to represent that field, and performed 
a recognition using the likelihood decision rules developed from the 14 fields. This method in- 
correctly replaces intrafield variations with interfield variations in deriving the decision rule. 
When the tests were run, it was found that a large number of fields were not recognized, in the 

uncorrected data, especially. The experiment was repeated, with a decision made by a criterion 

20 

that would reject five points in 10 if the likelihood functions accurately described the inter- 
field distributions. The results are shown in Table 4. The columns are the same as were used 
in Table 3. The rows represent correct, incorrect, and no decisions, respectively. One might 
expect the second row to be the same as the second row of Table 3 (majority rule), but it is not. 
The difference is caused by the nonlinear scan angle function and by the fact that the recognition 
process is nonlinear. 

4.3. COMPARISON BETWEEN SIMULATED AND REAL DATA 

We can now compare calculations with actual flight data. Model calculations were made for 
a sun position, an atmospheric state, and scanner geometry that matched the flight conditions as 
closely as possible. Figure 59 shows the ca-lihrated spectral radiance as measured by the 
scanner over five soybean fields on 4 September 1969, for an aircraft flying due west at an 
altitude of 1000 ft. The solid line is the result of our calculations for a wavelength of 0.56 pm 
normalized arbitrarily at a nadir scan angle of zero degrees. The fact that the actual data '■ v 
points do not approximate the almost linear condition represented by the model indicates that 
the assumptions used in the model do not completely represent real conditions. One could vary 
the state of the atmosphere, and a variation in the density profile would cause a greater change 
in radiance with scan angle, but it seems more likely that the omission of bidirectional reflec- 
tance in the model accounts for the difference in shape between the calculated and measured 
data. 

Figure 60 seems to indicate a more systematic, uniform variation in measured radiance 
data than that for the data in Fig. 59. This occurs at a shorter wavelength (0.42 pm) where 
more scattering takes place and hence we can expect the atmospheric effects to be more im- 
portant, that is, the scattering tends to mask the intrinsic variations caused by surface reflec- 
tances. 

Nonstandard atmospheric states can be taken into account by changing the vertical density 
profile of the aerosol distribution and various angular functions can be used to simulate bidirec- 
tional reflectance. Therefore, by making minor changes in the mathematical formulation of the 
model, to accommodate bidirectional reflectances, one should be able to remove all the known 
systematic variations in real scanner data. 


78 



—WILLOW RUN LABORATORIE 

TABLE 3. COMPARISON OF RECOGNITION RESULTS 

S — 



IK 

1KF 

DATA 

1KUV 

5K- 5K 

5K- IK 

Number of fields with probability of 
correct recognition less than 0.5 
(out of a total of 29) 

14 

13 

7 

14 

13 

Number of fields with wrong recognition 
using majority rule 
(out of a total of 29) 

4 

4 

1 

7 

5 

Average probability of no decision 

0.12 

0.32 

0.21 

0.01 

0.02 

Average ratio of the probability of an 
incorrect decision to the probability of 
any decision 

0.35 

0.20 

0.18 

0.4 5 

0.50 

Number of fields with more incorrect 
recognitions than correct decisions 
(out of a total of 29) 

12 

4 

4 

15 

13 

Number of training fields with the r . u _ 
probability of a correct recognition 
less than 50% 

(out of a total of 14) 

6 

4 

0 

5 

8 • 


TABLE 4. RESULTS OF THE USE OF FIELD MEAN FOR RECOGNITION WITH 
AN INTRAFIELD LIKELIHOOD DECISION RULE 



IK 

1KF 

DATA 

1KUV 5K-5K 

5K-1K 

Number of correct decisions* 1 

19 

21 

25 

10 

17 

Number of incorrect decisions 

10 

7 

3 

19 

12 

Number of data points not recognized 

0 

1 

1 

0 

0 


* Out of 29 


79 



TOTAL RADIANCE (mW/cm -sr-jim) 


WILLOW RUN LABORATORIES 



South North 

NADIR SCAN ANGLE (o) 


FIGURE 59. COMPARISON BETWEEN CALCULATED AND MEASURED RADI- 
ANCE DATA FOR SOYBEAN FIELDS (WAVELENGTH = 0.56 /um). Altitude = 

1000 ft; visual range = ~6 km. 




FIGURE 60. COMPARISON BETWEEN CALCULATED AND MEASURED RADIANCE DATA 
FOR SOYBEAN FIELDS (WAVELENGTH = 0.42 /urn). Altitude = 1000 ft; visual range = 

~ 6 km. 
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5 

INVESTIGATION lOTO THE SUITABILITY OF THE NORMAL LIKELIHOOD FUNCTION FOR 

RECOGNITION DECISION RULES 

Likelihood functions are used in classification decision (i.e. } recognition) processes on 

multispectral scanner data. These functions are usually represented by a multivariate normal 

(Gaussian) density function, statistical parameters of which are determined for the various 

decision classes from subsets of the data. Last year, tests were made of the normality of 

individual subsets of data corresponding to single fields, and all were found to be non-normal 

2 

at the 1% level of significance when a standard X goodness of fit test was used [ 12] . A review 
of that work, as well as a more detailed discussion which is summarized in the remainder of 
this section, is presented in Appendix III. A comparison of two maximum likelihood decision 
rules was made on the basis of paired receiver operating-characteristic. (ROC) curves, one 
member of each pair representing a multivariate normal decision rule and the other a rule 
based on an empirical density function. For each data set, an alternative hypothesis was 
assumed, and the Type I versus Type II errors (probability of miss versus false alarm) were 
plotted for different decision levels for each decision rule; see Fig. 61 for a typical pair of 
ROC curves. From such curves, a direct comparison can be made between the two likelihood 
function decision rules. 

The choice of an alternative hypothesis is an important consideration. If we consider the 
data to be points in a hyperspace of which each coordinate corresponds to a transformed 
spectral channel, then the question becomes one of where the alternative hypothesis will be 
located in the hyperspace. The location and shape, of the alternative distribution can be expected 
to affect the decision errors. A distribution was chosen that was uniformly located in the hyper- 
space. The use of this distribution corresponds to the use of many different separate distribu- 
tions located uniformly in the hyperspace. Thus, the results correspond to an average of the 
performance that we would obtain using a large number of separate distributions. This choice 
of an alternative distribution has the additional advantage of making it possible to test each data 
set individually. 

Using first one decision rule, then the other, we found the Type I errors by selecting a 
decision level and counting of the percentage of points that were rejected. The Type II errors 
were found by direct calculation. Data points for the alternative were assumed to be located 
uniformly throughout a hyper- rectangular parallelepiped, the dimensions of U’hich were set so 
that 0.9995 of the volume of the Gaussian distributions would be included. 

As a result of our tests, we have decided that the use of the normal likelihood function for 
individual fields is justified for recognition processing of multispectral scanner data. This 
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function is much quicker to generate and use than the histogram function. Also, the improve- 
ment in performance resulting from use of a histogram likelihood function and all of the 
channels is not believed to be significant. More promising approaches are: (1) preprocessing 
of the data to compensate for scan angle and similar systematic effects; and/or (2) compensation 
by changing the decision rule parameters. By using the first approach, we would expect the 
transformed data to have distributions that would be more nearly normal and, at the same time, 
have reduced-variances. . — - 
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6 

DATA -COLLECTION MISSIONS 

One primary and two secondary data-collection missions were flown during the summer of 
1971 to obtain multispeetral scanner data for future use in the testing and continued development 
of processing techniques to overcome the systematic trends introduced into scanner data by the 
atmosphere, sun position, and scan geometry during area survey operations. 

The primary flight was made between two missions flown for the 1971 Corn Blight Watch, 
and the secondary missions were carried out enroute from corn blight missions. Ground instru- 
mentation and reference reflectance panels were deployed during the primary mission. 

The chosen test site is located in Ingham County, Michigan, partly on and partly south of 
the Michigan State University Agi'icultural Farm. These areas were chosen because: (1) in- 
formation on the types of crops present was to be available from other sources; (2) they were 
conveniently located for coverage enroute to or from the Corn Blight flight line in Indiana; and 
(3) space was available for laying out the reflectance panels and setting up ground instrumenta- 
tion. 

The basic cycle of coverage included one pass at 5000 ft altitude over a 1 mile x 8 mile 
area south of Michigan State University (MSU), two adjacent passes over the MSU Agricultural 
Farm and campus, and three short passes,each at 1000 ft over areas of special interest. This 
cycle was flown on the secondary missions on 17 August and 21 September 1971, but was modi- 
fied somewhat for the primary mission. Three additional passes were made on the primary 
mission over the instrumented test site at altitudes of 1000, 2000, and 5000 ft. This expanded 
cycle was repeated three times, roughly at 9 a.m., 10:30 a.m., and 12 p.m. on 6 August 1971. 
Between these groups of passes, six consecutive passes were made over the 1 mile x 8 mile 
test area. These passes will provide data with which to evaluate the time -dependent changes 
occurring in the scanner signals and to campare these changes with model calculations. 

An additional set of passes was made along flight lines that are perpendicular to the 1 mile 
x 8 mile strip and extended over it. Each pass was displaced by one mile from the previous 
one, so that an area coverage was obtained in which each pass contains a portion of the area 
for which ground information is available. 

The ground instrumentation was used to measure spectral irradiance at the surface of the 
ground (0. 4-1.3 pm), spectral radiance from the reflectance panels, sky radiance at three wave- 
lengths, incoming global radiation over the 0.4-to 3 -pm wavelength interval, both outgoing short 
wavelength radiation and net radiation over bare soil and vegetation stands, surface radiometric 
temperatures, and meteorological parameters. 
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7 

CONCLUSIONS AND RECOMMENDATIONS 

7.1. SIMULATION AND MODELING OF SCANNER SIGNALS 

7.1.1. CONCLUSIONS 

The principal conclusions regarding the application ol the radiative transfer model are: 

(1) A model has been developed which accounts for the transfer of radiation through 
Earth's atmosphere and the model has been partially verified by comparison of 
calculations based on the model to experimental sky radiance data. 

(2) The present model is versatile enough to simulate a large number of conditions. Many 
of these conditions, which correspond to realistic flights, have been simulated and are 
discussed in this report. 

(3) The mathematical formalism has also been developed to enable us to correct the actual 
scanner data for atmospheric effects and surface bidirectional reflectance effects. 

(4) A potential advantage of the radiative transfer model as opposed to an empirical tech- 

... unique Ues-iiLthe saying of -process sing time. . .c\o. vtr 

Therefore, we conclude that by using the radiative transfer model one can gain greater insight 
into the geometric and physical aspects of the radiometric quantities involved in the analysis of 
scanner data. 

7.1.2. RECOMMENDATIONS 

The radiative transfer model developed during the past two years can be applied to the 
simulation of sensor response under realistic flight conditions. One can expect some discrep- 
ancies, (as seen in Figs. 59 and 60) to arise as a result of the inherent limitations of the model. 

To remove these limitations and therefore bring the model into closer agreement with actual 
conditions, one should consider several improvements in the extension of our model. 

We have compared radiances calculated from the radiative transfer model for different 
stratospheric aerosol densities; however, the change in radiance, irradiance, and transmittance 
caused by a large variation in the tropospheric vertical density profile has not been explored. 

The particle size distribution and the aerosol density profile should be varied for investigation 
of the corresponding variation in transmittance, irradiance, and radiance, to simulate unusual 
weather conditions. 
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A surface which is perfectly diffuse (Lambertian) in reflectance properties has been as- 
sumed for almost all work done in radiative transfer through 'planetary atmospheres. The 
results of some theoretical and experimental studies on the reflectance properties of natural 
surface conditions indicate, however, that the Lambertian assumption is not always valid. 

Many deviations between mathematical models and actual data could, in fact, be attributed to 
failure to consider bidirectional reflectance in the expanded radiative transfer models. 

The question of how much one surface element affects the radiance from a contiguous sur- 
face element has not been answered satisfactorily, primarily because of the mathematical com- 
plexities in the modeling efforts. Nevertheless, a simplified two-dimensional analysis of this 
problem has been carried out as described in this repo 2 *t. So far we have only considered single 
scattering in the two-dimensional analysis but it might be necessary also to investigate multiple 
scattering. 

The effect of discrete cloud formations has not been studied in a mathematical model of 
radiative transfer through Earth’s atmosphere. Previously, detailed studies of the transmit- 
tance and reflectance properties of clouds have been done at Willow Run Laboratories [ 13] and 
the results of this work should be included in an advanced radiative transfer model. 

For this report, we have investigated the interdependencies of irradiance, transmittance, 
Yadiahce, J atid 'Variafts-'-combinatibhW of the’s&on the parameters, visual range, wavelength,' scan 
angle, time, and surface albedo. We have seen how the various quantities are correlated for a 
particular atmospheric state. By varying the atmospheric condition, we should be able to dis- 
cern the most important features to be used in a comprehensive, unified radiative transfer 
model. 

Finally, for normal or standard hazy atmospheres predominantly composed of water, there 
is little absorption in the visual region of the spectrum. However, in certain industrial areas, 
a considerable amount of absorption by aerosol particles can take place. Therefore, we must 
investigate the overall effects of absorption on the transmittance, irradiance, and radiance in 
both the visible and near infrared spectral regions. 

7.2. TECHNIQUES FOR REMOVING SYSTEMATIC EFFECTS 

7.2.1. CONCLUSIONS 

A straightforward method of scan angle correction and signature extension has been devel- 
oped and tests have been conducted. The correction scheme should practically eliminate most 
systematic data variations caused by the presence of the atmosphere. The limitations seem to 
be noise present in the data and interfield reflectance variations. Theoretically, the 
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bidirectional reflectance patterns of the ground covers should also be a limitation if their angu- 
lar properties strongly depend on material type; in the tests conducted to date, this effect was 
not noticeable. In deriving the correction functions for scan angle correction, we selected 
areas of assumed constant ground cover to derive the correction function. For data set-to-data 
set signature extension, the corrections are derived by location of at least two pairs of areas 
for which one member of each pair is in each data set, the reflectances of each pair can be 
assumed to be identical, but the reflectances of the pairs are quite different. Once the areas 
are located, the scan angle correction and signature extension are straightforward calculations. 

An initial test was made to determine whether we could recognize fields by using the aver- 
age measured radiance for each field. This method shows promise for large area-survey appli- 
cations, although more experience is needed for an evaluation. 

7.2.2. RECOMMENDATIONS 

The studies performed have led to the following recommendations. 

(1) To account for systematic variations in scanner data, the technique of varying the 
parameters of the likelihood functions in the recognition process in correspondence 
with the.data measurement. parameters. and. conditions should be. investigated. Pre- 
processing was used in the reported tests of scan angle correction and it is an efficient 
method when all materials exhibit the same angular effects. This other method pro- 
vides theoretical advantages, when scan angle effects depend on material. It also pro- 
vides greater versatility and, thus, the potential for more sophisticated processing. 

For example, various forms of adaptive processing can be used to aid in recognition 

of intrafield variations, bidirectional reflectance differences between ground cover 
types, and slowly changing interfield variations. Perhaps, combinations of the two 
methods will prove useful. 

(2) The correction methods should be tested on other data, a logical choice being data col- 
lected over a larger area. With this choice, the one-pass multispectral scanner opera- 
tion could be extended to an area survey operation. For a large area, the correction 
also would have to handle gradual changes in the data not connected with scan angle; 
i.e., changes dependent on time and distance. 

(3) Interfield variations and methods for handling them should be studied with a larger 
data set than the one used for this study. After scan angle corrections, the recognition 
capability available at present appears to be limited by such variations. 

(4) The length of time required for recognition processing should be decreased. Using 
presently available digital methods, the time required for the recognition operation 
substantially limits the amount of data that can be processed in technique -development 
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and feasibility studies. Under samp ling, filtering, and channel selection are now com- 
monly used; additional methods are still needed. 

(5) There is a need to study data sets and model calculations specifically to determine 
whether radiance measured from the edges of fields differs sufficiently from that in 
central positions to warrant the immediate development of our proposed correction 
method for minimizing this effect of the atmosphere. 

(6) A model should be developed that would explain the sources of the terms in the co- 
variance matrix in order to increase understanding of the physical basis of the varia- 
tions in the data, as well as to better match the likelihood functions to the data. 

7,3. INVESTIGATION OF DECISION- RULE LIKELIHOOD FUNCTIONS ■ 

The use of the multivariate normal likelihood function for individual fields in recognition 
processing is justified because: (1) it is much quicker to generate and use than are histogram 
likelihood functions; and (2) the improvement that would result from the use of histogram func- 
tions instead is insignificant. More promising approaches are preprocessing of the data to re- 
move systematic effects and/or change of decision rule parameters. 
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Appendix I 

THE RADIATIVE TRANSFER MODEL 

This appendix contains a description of the most important aspects of a radiative transfer 
model for use in correction of remote sensing scanner data for hazy atmospheric conditions. 

1 . 1 . ELECTROMAGNETIC INTERACTIONS 

In radiative transfer theory, the most fundamental radiometric quantity is the spectral ra- 
diance for a particular state of polarization. In a vacuum, spectral radiance is an invariant, 
i.e,, it does not depend on distance. In traversing a medium, however, radiation is attenuated 
as a result of its interaction with the particles composing the medium. Some of the most com- 
mon electromagnetic interactions are the following: 

(1) Elastic scattering 

a. Rayleigh scattering by atoms and molecules 

b. Thomson scattering by free electrons 

c. resonance scattering 

‘ ( 2 ) : ‘ ' *" 

a. ordinary Compton scattering (energy loss) 

b. inverse Compton scattering (energy gain) 

c. Raman scattering by atoms and molecules 

d. Mie scattering by aerosol particles 

e. fluorescence 

(3) Capture 

a. line absorption 

b. photo-ionization 

c. inverse Bremsstrahlung 

d. photodLssociation of molecules 

e. pair production 

f. photonuclear reactions 

In the visible and near infrared regions of the electromagnetic spectrum, only Rayleigh scat- 
tering by atoms and molecules and Mie scattering by aerosol particles are important for remote 
sensing applications. Any aerosol absorption which takes place in this spectral region will vary 
slowly with wavelength as opposed to the rapid change in absorption properties in the case of 
gaseous absorption. 
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I. 2. THE ATMOSPHERE 

Earth's atmosphere consists primarily of gas molecules such as 02? Ng* H^O, ^^2' and ^3 
as well as aerosol particles, i.e., much larger groupings of molecules having a definite shape. 
Typical aerosol particles would be water droplets, dust, smoke particles, or any other semi- 
permanent suspension in the atmosphere. 

For visible and near infrared radiation, the density fluctuations in the atmosphere act as 
scattering centers, and since the wavelength of the radiation'is much greater than the size of a 
scattering center, an oscillating electric dipole radiation is created. This is called Rayleigh 
scattering. The dependence of the intensity of the scattered radiation on wavelength and scat- 
tering angle is as follows: 

Intensity «X 4 (1 4 cos 2 Q) .... (53) 

where X is the wavelength of the radiation and 8 is the angle between the direction of the scat- 
tered and incident radiation. It is this strong dependence on wavelength which led Lord Rayleigh 
to conclude that the scattering of sunlight by scattering centers in the atmosphere accounts for 
the blue color of the sky. Although this type of scattering does account for the general blueness 
of the sky, certain serious discrepancies occur when one attempts to use Eq. (1) to describe the 
sky radiation. in detail. In general, the spectral and ; angular dependence ; of the ; s ky radiation is 
not given by Eq. (1) and the polarization points in the sky do not correspond exactly with those 
as determined from Rayleigh's analysis. Two explanations can usually be given for these dis- 
crepancies: (1) Rayleigh's theory considered only single scattering, but multiple scattering can 
be important, especially in the visible spectral region; and (2) scattering by larger particles 
(aerosols) was neglected. 

During the past 25 years, much work has been done to remove the discrepancy caused by 
consideration of single scattering only. In particular, the efforts of Chandrasekhar [9] and 
Mullikin [ 14] have resulted in the climatic achievement of providing an exact formulation for 
the radiation field in a homogeneous*, plane -parallel atmosphere with multiple scattering 
according to Rayleigh's scattering law. Although the mathematical formalism is rather in- 
volved, Coulson et al. [8] performed a computer analysis of the equations and determined the 
radiation field emerging from the top and bottom of a plane-parallel, homogeneous, Rayleigh- 
type atmosphere. Additional work by Bellman and Kalaba [ 15] and Grant and Hunt [ 16] have 
yeilded results which give the radiation field in an inhomogeneous atmosphere with Rayleigh 
scattering. 

* Homogeneous means that the scattering properties of the atmosphere do not depend on 
altitude. 
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Unfortunately, the atmosphere is never free of particulate matter, and there is evidence 
that the aerosol component is increasing with time [ 17, 18}. Thus, computations based on a 
Rayleigh-type atmosphere cannot accurately represent the real Earth's atmosphere. In fact, 
in the case of hazy atmospheres, most of the scattering is caused by aerosol particles. The 
theoretical analysis of nonRayleigh-type atmospheres has not met with the same degree of 
success as for Rayleigh-type atmospheres, primarily because of the nature of scattering by 
aerosols. Instead of following the angular dependence of dipole scattering as in the Rayleigh 
case, the angular dependence of radiation scattered by particles the sizes of which are com- 
parable to the wavelength of the incident radiation, is a more rapidly changing function of the 
scattering angle. This strong dependence on angle is illustrated in Fig. 62. The high degree 
of anisotropy characteristic of aerosol scattering is very difficult to deal with mathematically 
and certain approximations are needed to solve the problem of radiative transfer in a plane- 
parallel, homogeneous aerosol atmosphere. 

Assuming that aerosol particles are basically spherical in shape, their scattering, absorp- 
tion, and total cross sections and the angular scattering properties (scattering phase function) 
can be calculated according to classical electrodynamics. This is called Mie scattering, the 
mathematical development of which is given by Stratton [ 19] , van de Hulst [20], and Kerker 
[21], If one knows the particle size distribution, then one can easily determine the scattering, 
absorption, and extinction coefficients as well as the complete scattering phase function for the' 
polydispersion of spherical aerosol particles. We shall designate the absorption, scattering, 
and extinction coefficients by a{ z), |3(z), and k (z) respectively, where z is the altitude above sea 
level. 

One quantity of great importance in estimating the scattering or absorbing properties of a 
medium is the optical depth, r. It is defined as 

oo 

r =j*(z)dz (54) 

h 

where h is some base altitude above sea level. The total optical depth of an atmosphere is 
given by 

GO 

T 0 = jjc(z)dz (55) 

0 

Thus, r is a variable optical depth, varying from zero at the top of an atmosphere to its total 
value, Tq, at the bottom of an atmosphere. For Rayleigh atmospheres, lies between 0.1 and 
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1.0 in the visible and near infrared spectral regions, whereas for aerosol atmospheres Tq can 
be considerably greater than 1.0. The physical significance of optical depth is that it is a 
measure of the number of interaction lengths or mean free paths of a photon of a certain energy. 
Therefore, a large optical depth, which occurs at short wavelengths, implies an optically thick 
medium and hence a large number of scatterings. A small optical depth, such as 0.1, which 
usually occurs at long wavelengths, indicates that few scatterings take place. An example of the 
relationship between optical depth and altitude for a standard atmosphere is shown in Fig 63, 
and is shown with the percent deviation in Fig. 64. Further details on the physical of Earth’s 
atmosphere can be found in Bullrich [22], Fleagle and Businger [23], Junge [24], and Robinson 
[25]. 


1.3. VISIBILITY 

To remove the effects of the atmosphere on the radiation received by an airborne detector, 
we need to define the state of the atmosphere. One measure of the atmospheric state is the 
visibility or visual range. 

If a black object is moved away from an observer in a horizontal plane, the contrast be- 
tween that object and the sky horizon background gradually diminishes until a limit is reached 
beyond which the object will no longer be distinguishable. The apparent brightness of the black 
object at a horizontal distance, x, from an observer is 

B Q = B b (l - e' ,fX ) (56) 

where B^ is the brightness of the horizon and k is the extinction coefficient. The contrast (C) 
is defined as 


B, - B_ 

c - 2 . - e ~ KX 


(57) 


The convention followed for the limiting contrast is that given by Koschmieder [4] as 2%, i.e., 
C = 0.02. Thus, the range corresponding to this contrast is 


= V = 


In 50 

K 


or 


(58) 


V - 3912 
K 


(59) 


Equation (59) is then a simple formula relating horizontal visual range (km) to the atmospheric 
extinction coefficient. The clearest day possible would be one in which there were no 
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particulates. At a wavelength of 0.55 pm, the corresponding Rayleigh extinction coefficient is 
- 2-1 

1.162 x 10 km and hence the visual range is 336 km. Actually, however, some haze is al- 
ways present and 23-km visual range is taken by Elterman {5] to represent the line between a 
clear and a hazy atmosphere. Excellent pictures of the effect of haze on seeing ability from an 
aircraft are given by Larmore and Hall [26] for a variety of conditions. The line of demarcation 
between a haze and a fog is usually taken to be 1.2-krn visual range. A more complete treatment 
of visibility is given by Middleton [27], and Johnson [28]. 


1.4. THE RADIATIVE TRANSFER EQUATION 

In this section, we shall discuss briefly the radiative transfer equation which is used to de- 
termine the radiation field within a plane-parallel, homogeneous atmosphere with aerosol scat- 
tering. A more complete discussion and derivation of the equation is given by Turner [ 2] and 
Malila et al. [1]. 

The basic integro-differential equation of radiative transfer for a plane -parallel, homo- 
geneous atmosphere illuminated by solar radiation is given by 


2rrl 


dL 

‘df 


L f f w n 

f “ L ("> \ I P (/*, <t>y p 1 , 4> r )dp’d<£’ - E g (r)p(p, <A, -p QJ <£ Q ) 


0 -1 


- (1 - u)q)B(t) 


(60) 


where 


E g (r) = E 0 e" r//I 0 (61) 

L(t, m, $) is the spectral radiance at optical depth r, zenith angle 9 (of which the cosine is p), 
and azimuthal angle <f>. E s (t) is the solar irradiance at optical depth r and p(p, <f > , p', <£') is the 
single -scattering phase function which describes the function of energy scattered from the di- 
rection p', <f>' into the direction p, B(r) is the Planck radiation function and is called the 
single scattering albedo defined as 


% = K ( 62 ) 

The direct solar radiation enters Earth's atmosphere with a zenith angle the cosine of which is 
Mq and azimuthal angle <j>^. 

For the visible spectral region, B(r) is usually negligible and o.' Q = 1 since there is very 
little absorption by either gases or particulate matter. Nevertheless, Eq. (60) is quite difficult 


95 



WILLOW RUN LABORATORIES 


to solve: exactly and has been done only for isotropic and Rayleigh-type scattering. For realistic 
atmospheres with an aerosol component, the phase function is highly anisotropic and approxima- 
tions must be made to solve the transfer equation. One approximation to simplify Eq. (60) is 
the following: 

p(p, 4>, V', 4>') = F 6(p - p'W* - </>') + B6(/i + + <f> - <*>’) ( 63 ) 

where F is the fraction of energy scattered into the forward direction and B is the fraction 
scattered into the backward direction. Since scattering by aerosols is strongly peaked in the 
forward direction, this approximation seems to be a reasonable one. Use of Eq. (63) in Eq. (60) 
permits a solution in terms of arbitrary constants which are determined from boundary condi- 
tions. The general boundary conditions are, 

(64) 

L(Q. -M. 0) = 0 


277 1 

L(t q , p, <*►) | V-'p'iV, 4>> <#»')[ L ( r 0 > 

0 0 


-p‘, <£’) + L £ (t 0 , -a’, </> f )j<V d<f> f 


(65) 


where Eq. (64) simply states that there is no diffuse radiation entering the top of Earth's at- 
mosphere and Eq. (65) indicates that we must integrate the diffuse and solar radiance with the 
bidirectional reflectance over the hemisphere of incoming radiation. If we deal with a Lam- 
bertian (perfectly diffuse) surface, the boundary conditions reduce to 


E JO) = 0 

W * p[ e - (t o’ + Vo e ‘ To/ "°] 


( 66 ) 

(67) 


where E (0) is the downward diffuse irradiance at the top of the atmosphere, E_(t q ) is the down- 
ward diffuse irradiance at the bottom of the atmosphere, E + (Tq) is the upward diffuse irradiance 
at the bottom of the atmosphere, Eq is the direct solar irradiance at the top of the atmosphere, 
and p is the hemispherical reflectance of the surface. 

Breaking the radiation fields up into two components, an anisotropic one for p = 0, and an 
isotropic one for p * 0, we can solve for the irradiances at any point within the atmosphere, i.e., 


/a M 0 E 0 

*V 7) " p Q + (1 * T])r 0 


. . . 1 + 2(1 - 7])t 

(i - 7]){r 0 - t) + pa 0 


(68) 


^o E o 




(1 - v)(tq “ r) - [ A 0 + (1 - ^Tq]' 


-t/„ 0 + 


1 +2(1 - ij)(l -p)r 


(69) 


0 
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E (t) = 


/i o E o 


+ (1 - vWq 


2M 0 P(1 “ t?)7 


M 0 + (1 - V)(r 0 - r) + ! t 2( i . i| 1(1 - P) r 0 


where V is f/4tt and E (r) is the total downward irradiance. 


(70) 


Having determined the irradianees, we can now find the radiances by using the approxima- 
tion 


E”(r) + E M (r) 


L(r, p, 0) ~ [ E^(t) 6 (p - p. 0 )d(7r + 0 Q - 0) + E^(r)$(p 4 m 0 ) 6 (^ - ^ 0 ) j + 2n~ 


(71) 


where the primed irradianees represent the radiation field with p = 0 and the double primed 
irradianees represent the radiation field with p * 0. The complete spectral path radiance in 
the upward and downward hemispheres are then, respectively 

T , . ^0 

V" ® 47t[p 0 + (! - 7)r 0 ] 

^j(l - T) )Tq{p(p, 0, P 0 , 71 + 0 Q ) + p(p, 0, -P 0 , 0 Q )] + P 0 p(p, 0, -p Q , 0 q ) 


2 PoP 


+ 1 + 2 (I - V )Ti" -~p )Tr 


1 - e 


-(^0 -r)/p 


+ j(l " V)[p{H, <P, P 0 , it + 0 Q ) 


8(1 - ??)pjjp 1 

+ P(P> & -M 0 > ^ 0 )] - r+‘2(l - 7})(1 - p)tq 


-( r o -f)/p 

(t q 4 p)e -(r + p) 


(72) 


L p (r, -)i, t>) = 4,[n 0 + (1 - t))r 0 ] 


(4. - ’))7 0 [p(-M, 0, Pq. »■ + + p(-M, <?>, -M 0 , 0p)] + P 0 P(-M, -P 0 , <?> 0 ) 


(equation continued) 
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+ i+2ft a ~-W^ '} < x - e ' T/M ) - {<> - *»*■<-**. * V * + V 

8(1 - 7J )nh) l /, \ 

+ p(-M, <h -V <%>] - ! + 2(1 - m - pSt’q j <Me +T -<‘>) 

where the single-scattering phase functions are given by 


pO*. 0, h 0 , JT + 0 O ) ^p^/i 0 - Vd " 

- Pq ^cos(0 - 0 O )J 

(74) 

p(p, -V 0 > %) =p[“PP 0 + V( x " ^ 2 )( 

1 - j^)eos(4> - 0 O )J 

(75) 

p(-M, 4>, h 0 , IT + 0 Q ) = p^-pMq - V (l - ^ 2 )( 

1 - ^jcos(0- 0 O )J 

(76) 

p(~P, <$>, -Vq, <Pq) - p[/i/i Q + V(l - M 2 )(l 

- Mq)cOs( 0- 0 o ) ] 

(77) 
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Appendix II 

DERIVATION OF AN EMPIRICAL TECHNIQUE 
FOR CORRECTING HAZE EFFECTS NEAR BOUNDARIES 

Analyses of scanner data and radiative transfer model calculations have shown that the 
material surrounding the material being observed contributes to the received spectral radiance. 

[2]. It is because of the atmosphere that some of the radiation reflected from the surrounding 
materials is scattered directly into the scanner’s field of view and some irradiates the element 
being viewed and is reflected to the scanner. Both components are detected as part of the tar- 
get's radiance and are indistinguishable from it. If a white or grey panel is surrounded by 
green vegetation, the spectral content of the path radiance resembles that of a green material 
and, therefore, the effect has been called green haze. 

This appendix describes a mathematical technique which can be implemented empirically 
to remove the contribution of the surrounding materials to the observed radiance when they 
differ from the material being observed. The first step is to find a boundary between two dif- 
ferent fields. If we choose x, y coordinates, so that the boundary is on the y axis, a plot of the 
measured radiance across the boundary might be that shown in Fig. 65. Here a and b are the 
measured asymptotic values of radiance. We then proceed as follows: 

(1) Form f(x), where: 

f(x) = irb-^[L 0 ( x ’ y) -a] (78) 

(2) Find F(to), the Fourier transform of f(x). 

F(cc) = j*e ^^(xjdx (79) 

(3) Form 

K<«) = G(o>) (80) 

where G(u>) is a low pass filter corresponding to the desired system bandpass. 

(4) Find the inverse transform of K(w). 

CO 

k(x) = ~j ejwXK (cu)da> (81) 


~ oo 
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(5) Transform the data by the equation 

L/(x, y) =J|l(x - u, y - v)k(l/^ +v^)dudv (82) 

The above algorithm has the property that (1) if the radiance can be written as: 

. . (83) 


L(x, y) = ETp(x, y) ± L p + |jETp(x', y T )h(x *x f ,y - y'Jdx’ dy' 


where h(x, y) is function that represents the spatial spreading effect of the atmospheric haze, 
and (2) if the measurements shown in Fig. 63 were of a surface with reflectance 


P 0 (x, y) = p a , x < 0 


= p b , x > 0 


then 


L’(x, y) - [1 + H(0, 0)]ETp(x, y) + L 


(84) 


(85) 


where H(o> , w ) is the Fourier transform of h(x, y). 
- x y 


H(0, 0) = f! h(x, y)dxdy, a constant, 
and h(x, y) is assumed to have radial symmetry, i.e., 


h(x, y) - H(]/j?+y*) (87) 

We can demonstrate the property mentioned above by assuming Eqs. (83) and (84). We 
first write p Q (x, y) in a one-dimensional form: 

p Q (x, y) = (p b - P a )U(x) + p a (88) 

where U(x) is the unit step function. From Eq. (83) we have, for the two steady state values of 
L 0 (x, y): 

a = ETp + L ft + ETp H(0, 0) 
s u s 


( 86 ) 


= ETpJl + K(0, 0)] 4 L p 


(89) 
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b - ETp b [l + H(0, 0)] + L p 

Combining Eqs. (89) and (90) 

b - a =ET(p b -p a )[l +H(0, 0)] 

From Eqs. (83), (88), and (89) we derive: 

L Q (x > y) - a = ETp 0 + L p + ETp Q *h(x, y) - ETpJl + H(0, 0)] - L p 


= ET(p b - P a )[U(x) + U*h x (x)] 
Now from Eqs. (78), (91), and (92): 

6(x) + h^x) 

f < x) = 1 TlTjW" 

Thus, from Eqs. (79) and (80) we find: 


1 +H (0) 

K(w) * l i 


Let us now consider the use of K(oj). From Eq. (82) we see that: 


(91) 


(92) 


(93) 


(94) 


l/ - L l KM < 95 > 

where l} denotes the Fourier transform of L. From Eq. (83): 

l} = ET p X + L 6{ w ) + ETp^! 

P 1 

ETp*(l + H) + L 8 (oj) (96) 

r 

After applying K(u>), we have 

L’ 1 s^ETp^l + H 1 (0)] + L5 (w)]g(w) (97) 

which has the inverse transform of Eq. (85) for G(w) = 1. This completes the development. 

One additional point should be mentioned. The data are unchanged in the middle of large 
fields, since the green haze effect is a local phenomenon. 
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Appendix III 

AN EMPIRICAL INVESTIGATION INTO THE 
SUITABILITY OF THE NORMAL LIKELIHOOD 
FUNCTION FOR PROCESSING MU LTISPECTRAL 
SCANNER DATA 


by 

R. B. Crane 
W. A. Malila 
W. Richardson 

Abstract 

Likelihood functions are usually used in classification decision (Le., recognition) processes 
on multispectral scanner data. These functions can be represented by a multivariate normal 
(Gaussian) density function, the statistical parameters of which are determined for the various 
decision classes from subsets of the data. A comparison was made to determine whether or 
not. improved classification results, could be obtained by use of a different form to represent 
the likelihood functions, namely, an empirical, multivariate probability, density histogram of 
decorrelated variables. First, tests were made of the normality of the individual subsets of 
data, and all were found to be non-normal at the 1 % level of significance when a standard chi- 
square goodness-of-fit test is used. Operating characteristic curves then were generated to 
represent decisions made with each form between each given class and a uniformly distributed 
alternative class; a uniform distribution was chosen because the results are then least depen- 
dent on the choice of the alternative data set. It was found that the probabilities of miselassi- 
fication with the two forms were approximately identical for almost every data set, even though 
a large number of individual data points were classified differently. It was concluded that the 
conventional assumption of a multivariate normal distribution of multispectral scanner signals 
from individual fields is sufficiently accurate to warrant its use in recognition processing 
rather than a more complicated empirical distribution function. 

m.l. INTRODUCTION 

A multispectral scanner is a mapping instrument that collects and stores infrared, visible, 
and ultraviolet images in electrical form on magnetic tape. The scanning itself is similar to 
that employed by an infrared scanner, whereby the circular scan of the collecting optics and the 
forward motion of the scanner provide the two-dimensional scanning pattern. The images col- 
lected by the optics are then separated into wavelength bands by a prism, grating, or 
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interference filters before detection and storage for subsequent recognition processing. Either 
before or after storage, the signals can be digitized by an appropriate sampling procedure, if 
digital processing is to be used. 

The likelihood function plays an important role in the processing of multispectral scanner 
data. It is used in the Bayesian decision rule, as well as more classical rules such as the 
Neyman- Pearson, minimax, and maximum likelihood. The likelihood of an observed set of 
values is the joint density function for continuous data or the joint probability function for dis- 
crete data. In this appendix, we restrict our attention to discrete data. For the data that are 
considered, the likelihood functions are not known, but must be estimated from the data set it- 
self. Subsets of the data are chosen from which various likelihood functions are derived, each 
function corresponding to a separate decision outcome. Such subsets are often known by addi- 
tional information to represent materials to be recognized and are therefore called training 
sets. These tests will not include the possibility of interfield variations which cause differences 
between the statistics of the training sets and those of other fields of the same class. A con- 
venient approximate method of calculating the likelihood function from the data subsets is made 
possible by the assumption that the data are generated by a Gaussian or normal process. The 
normal likelihood function is completely determined by first- and second-order moments of the 
data subsets. Although the normality assumption can sometimes be justified by an application 
^of. the,, central Timit theorem,, -that argument, .cannot be justified here. : - ■ 

The alternative to the normal assumption that we will consider is that of constructing, from 
histograms, an empirical likelihood function for each subset of data. Each data point of a subset 
of data, consisting of the values from the spectral channels for each particular sample, is lin- 
early transformed into a new set of values. The transformation is chosen so that the values in 
any pair of channels are uncorrelated. Histograms are then formed for each of the new chan- 
nels, and these histograms, properly normalized, become the new approximation to the likelihood 
function. 

Practical considerations led to the use of this particular method of deriving a likelihood 
function. The likelihood function should represent the data that are to be analyzed; hence, it is 
derived from the same data. The subset has a finite number of members, as would any subset 
used for multispectral analysis; thus, one would not expect exact agreement between data and 
the underlying likelihood function. Finally, the histogram method can be justified if two assump- 
tions concerning the data are made: 1) the channels are independent when they are uncorrelated; 
and 2) each histogram accurately describes the marginal distribution. Because data from sev- 
eral sources are to be analyzed, the effect of the first assumption, at least for many of the data 
sets, can be expected to be small when compared to the effects of using the Gaussian assump- 
tion. In fact, the first assumption would be true if the data were distributed normally. The 
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second assumption is justified by the limited data analyzed. Since there is going to be an un- 
certainty in the estimation of each marginal distribution, that uncertainty is reflected in the 
number of intervals of the histograms. The histogram describes the distribution function as 
accurately as the measurements permit. 

m.2. DESCRIPTION OF TEST PROCEDURES 

The multispectral scanner data analyzed were obtained from agricultural fields overflown 
in two different years, at different locations. The set of 1966 data was collected over one grow- 
ing season in Indiana, on and near the Purdue University Agricultural farm. The 1969 data set 
was selected from one run over the Imperial Valley in California. Each set of data corresponds 
to one field, with one agricultural crop. Care was taken to collect data from the interior of the 
fields, and not to include any boundary samples. Since the present method of processing multi- 
spectral scanner data assumes that the data are from a normal process, the marginal distribu- 
tion of amplitudes of the original data sets were tested for non-normality. In addition, similar 
tests were made of the data after transformation by a normalized eigenvector basis. 

The non-normality test is based on histograms of the amplitudes in each channel or in each 

transformed variable. By choosing the intervals for the histogram (frequency distribution) so 

that each interval is equally likely to be populated by a normal distribution, one can readily 

see any significant departures from a normal population. To quantify the non-normality of the v 
2 

distributions, a x test is applied to each histogram. 

2 

The X test for non-normality (goodness of fit) is one of several that could be applied. It 

was selected because it is convenient to apply, has a well defined, asymptotic distribution, and 

2 

provides a quantitative criterion for rejecting the normal hypothesis. The usual x statistic 
(called X^) was calculated, and its significance (i.e., the probability of an this large or 
larger arising from a normal distribution) was determined from the x distribution, with the 
appropriate number of degrees of freedom. The lower this probability, the less likely there is 
a normal distribution; we have used a value of 0.01 as a decision threshold for rejection. 

A linear transformation of the signals does not affect the normality of their joint-probabil- 
ity distribution, but it can produce a new basis of uncorrelated variables (vectors) for describing 
the distribution. There is more than one transformation that will produce uncorrelated vari- 
ables. One, a Cholesky decomposition into triangular matrices, is used in our recognition- 
processing algorithms. Another, the eigenvector transformation, rotates the coordinates and 
aligns them with the major axes of the signal distribution in N-space for Gaussian data. We 
decided to apply the eigenvector transformation to the data to remove any interchannel cor- 
relation, in order to perform a scaling for producing unit variances and to test the marginal 

o 

distributions of the transformed variables for univariate non-normality using the x test. The 
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test does provide us with quantitative estimates for deciding when the transformed variables 
are not normally distributed, but normality of the marginal distributions is not theoretically 
sufficient to prove that the overall distribution is multivariate normal, even when the variables 
are uncorrelated; they must be shown to be independent as well. 

After the data sets were tested for normality and histogram likelihood function derived, a 
comparison was made between the use of normal and histogram likelihood functions in a deci- 
sion process. The comparison was made on the basis of receiver operating characteristic 
(ROC) curves. For each data set, an alternative hypothesis was assumed, and the Type I versus 
Type U errors wore plotted for different decision levels. The Type I and Type II errors aie the 
errors made when the hypothesis and alternative hypothesis, respectively, are present. From 
these curves a direct comparison could be made between the two likelihood function decision 
rules. 

The choice of an alternative hypothesis is an important consideration. If we consider the 
data to be points in a hyperspace, of which each coordinate corresponds to a spectral channel, 
then the question becomes one of where the alternative hypothesis will be located in the hyper- 
space. The location and shape of the alternative distribution can be expected to affect the deci- 
sion errors. A distribution was chosen that was uniformly located throughout the hyperspace. 
Use of this distribution corresponds to the use of many different separate distributions located 
uniformly in the hyperspace. Thus, the results correspond to an average of the performance 
that we would obtain using a large number of separate distributions. This choice of an alterna- 
tive distribution has the additional advantage of making it possible to test each data set indi- 
vidually. 

III. 3. SUMMARY OF EXPERIMENTAL RESULTS: NON-NORMALITY TESTS 

An examination was made of each channel of the data subsets recorded throughout the sum- 
mer of 1966 from 38 fields, under a variety of conditions. It was found that 310 of the 456 chan- 

2 

nels, or 68%, were non-normal according to our x test at the 1% significance level. The 
results of these tests and the conditions of measurement are summarized in Table 5. All of 
the channels tested normal for only one field, and that one had two transformed margins that 
tested non-normal. None of the fields tested was multivariate normal, as measured at the 1% 
level of significance. 

In a similar manner, we examined data recorded in 1969 from 16 fields in the Imperial 
Valley. We found 88 out of 160 channels, or 56%, were non-normal. None of the fields had all 
channels normal. Table 6 shows the number of non-normal channels in each field. 

In carrying out the tests, we were careful to account for the discreteness of the data. Be- 
cause the data values were in the form of 9 -bit integers, the histogram intervals were chosen 
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Altitude 


Approximate 
Angle Bet wet 
Scan Line an 


Crop 

Tape 

File 

Date 

(ft) 

Sun’s Dirt 

Corn 

1 

16 

7/29 

700 

5° 

Soy 

1 

19 

7/29 

700 

85° 

Corn 

1 

20 

7/29 

700 

85° 

Wheat 

2 

16 

6/30 

700 

5° 

Wheat 

2 

18 

6/30 

700 

85° 

Wheat 

3 

15 

6/30 

2000 

5° 

Corn 

3 

16 

9/15 

700 

10° 

Soy 

3 

18 

9/15 

2000 

10° 

Corn 

3 

19 

9/15 

2000 

10° 

Soy 

3 

20 

9/15 

4000 

10° 

Corn 

3 

21 

9/15 

4000 

10° 

Soy 

3 

22 

9/15 

6000 

10° 

Corn 

3 

23 

9/15 

G000 

10° 

Corn 

3 

25 

9/15 

8000 

10° 

Soy 

4 

5 

6/29 

700 

0° 

Corn 

4 

6 

6/29 

700 

0° 

Corn 

4 

7 

6/29 

700 

45° 

Corn 

4 

16 

6/30 

700 

0° 

Corn' 

4 

^o- 

6/30 

700 

90° 

Soy 

4 

21 

6/30 

700 

90° 

Corn 

4 

22 

6/30 

700 

45° 

Wheat 

4 

25 

6/30 

700 

0° 

Wheat 

5 

4 

6/30 

700 

45° 

Corn 

5 

9 

7/27 

700 

85° 

Sov 

5 

8 

7/27 

700 

85° 

Corn 

5 

12 

7/28 

700 

90° 

Soy 

5 

15 

7/28 

700 

90° 

Com 

5 

16 

7/28 

700 

5° 

Corn 

5 

18 

7/28 

700 

40° 

Corn 

5 

25 

9/15 

700 

85° 

Corn 

6 

2 

9/15 

700 

15° 

Corn 

6 

7 

6/30 

700 

50° 

Corn 

6 

n 

6/30 

700 

40° 

Corn 

6 

13 

6/30 

700 

5° 

Wheat 

6 

16 

6/30 

700 

50° 

Wheat 

6 

17 

6/30 

700 

40° 

Wheat 

6 

20 

6/30 

700 

5° 

Wheat 

6 

21 

6/30 

2000 

50° 
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IN EMPIRICAL DISTRIBUTION STUDY 





Number of 

Number 


Eigenvalues 


Negative 

of Non- 




Components 

Normal 

\ 

X 

V *2 

in Principal 

Margins 


*2 

Eigenvector 

X 

_JL 

0.297 

0.0206 

14.42 

0 

12 

2 

0.792 

0.0363 

21.82 

0 

12 

2 

0.126 

0.0220 

5.73 

0 

3 

1 

0.487 

0.0106 

49.94 

0 

12 

2 

0.0542 

0.0104 

5.21 

0 

8 

1 

0.186 

0.0939 

1.98 

0 

12 

5 

0.191 

0.0315 

6.06 

0 

5 

2 

0.161 

0.0362 

4.45 

0 

9 

1 

0.0403 

0.01 86 

2.17 

0 

0 

2 

0.0964 

0.0401 

2.40 

0 

1 

2 

0.198 

0.0373 

5.31 

6 

5 

2 

0.115 

0.0479 

2,40 

5 

1 

2 

0.148 

0.0381 

3.88 

5 

1 

2 

0.0993 

0.0282 

3.52 

5 

1 

1 

0.0898 

0.0255 

3.52 

0 

12 

6 

0.0175 

0.0107 

1.63 

0 

8 

3 

0.0235 

0.0088 

2.67 

5 

9 

6 

0.0370 

0.0073 

5.07 

1 

12 

2 

0.0108 

0.0028 

3. 86 

0 

2 

0 

0.0599 

0.0189 

3.17 

0 

12 

3 

0.0302 

0.0058 

5.21 

0 

11 

2 

0.339 

0.0026 

130.4 

0 

12 

4 

0.124 

0.0033 

37.6 

0 

12 

5 

0.0433 

0.0119 

3.64 

0 

8 

1 

0.303 

0.0432 

7.01 

0 

11 

6 

0.0884 

0.0073 

12.1 

0 

6 

3 

0.0537 

0.0087 

6.17 

0 

3 

2 

0.394 

0.0110 

35.8 

0 

31 

6 

0.0031 

0.0016 

1.93 

0 

5 

0 

0.151 

0.0180 

8,39 

0 

11 

5 

0.145 

0.0355 

4.08 

0 

12 

3 

0.0545 

0.0256 

2.13 

0 

12 

2 

0.0221 

0.0079 

2.80 

2 

8 

2 

0.0447 

0.0085 

5.26 

2 

11 

2 

0.0227 

0.0018 

12.6 

0 

11 

2 

0.0905 

0.0068 

13.3 

0 

12 

2 

0.336 

0.0066 

50.9 

0 

12 

3 

0.0209 

0.0019 

11.0 

0 

7 
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TABLE 6. SUMMARY OF 1969 DATA USED IN EMPIRICAL 
DISTRIBUTION STUDY, (imperial Valley, California; 
3/12/69; 1025 hours PDT; altitude = 5000 ft; 
east-west flightline.) 


Crop 

Field 

Number 

Principal Eigenvalues 
X 1 X 2 

Number 
of Non- 
Normal 
Margins 

x y 

Rye 

21 

0.0342 

0.0193 

10 . 5 

Alfalfa 

29 

0.0289 

0.0103 

8 3 

Sugar 

39 

0.0139 

0.0073 

2 2 

Soil 

45 

0.0235 

0.0110 

6 1 

Soil 

76 

0.0498 

0.0102 

10 3 

Soil 

77 

0.0135 

0.0053 

2 

Barley 

75 

0.0254 

0.0077 

2 

Rye 

- 91' 

0.0312 - 

0.0055 

4 

Lettuce 

155 

0.0179 

0.0102 

3 

Soil 

156 

0.0631 

0.0076 

6 

Alfalfa 

180 

0.0503 

0.0161 

9 

Sugar 

179 

0.0111 

0.0087 

2 

Sugar 

190 

0.0088 

0.0059 

1 

Lettuce 

191 

0.0461 

0.0120 

6 

Alfalfa 

205 

0.0417 

0.0112 

10 

Barley 

202 

0.0125 

0.0113 

7 
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to correspond to an integral multiple of possible data values. The hypothesized normal distri- 
bution in each channel was determined by the mean and variance of the sample, but extreme 
points introduced during digitization were excluded. In the case of the 1969 data, the dynamic 
range of the digitized values was small enough that quantization noise can significantly affect 
the non-normality tests. 

There is no typical pattern of distribution found among the non-normal channels. Figure 
66 shows one of the normal distributions; the extremes of the distribution have been compressed 
because of the method of presentation in intervals of equal probability. Figure 67 shows the 
distribution of a margin of the data which has a pronounced scan angle effect. The shape can be 
explained, at least in part, by the change of mean signal with scan angle. Figures 68 and 69 are 
examples of unsymmetrical and bimodal distributions, respectively. 

The analysis of the transformed data yielded some interesting results. The transformation 
was an eigenvector transformation, with the vectors ordered according to decreasing eigenvalue 
and weighted so that the transformed variables had unit variance. Thus, the first transformed 
variable corresponds to variation in the direction of maximum uncertainty in the untransformed 
hyper space. 

Figure 70 is a representative sample of the results of the analysis. It contains the analysis 
of the first two transformed variables for one of the fields of corn. Let us consider the analysis 
of the first variable, y(l). The bar chart on the left shows the number of data points that had 
amplitudes within each of twenty ranges. These ranges were chosen so that if y(l) were nor- 
mally distributed, as the number of sample points was increased, the percentage of points in 
each range would approach 5%. Therefore, the bars would have equal lengths. The first and 

last ranges were subdivided into two intervals each, so that any extremely unlikely points could 

2 

be shown separately and excluded from the x test. Thus, the first two lines form the first 
equally likely interval, and the last two lines form the twentieth. The intervals that were used 
are shown to the left of the bars in terms of standard deviations from the mean. 

In this example, because it is quite apparent that the bars are not of equal length, we might 
say that the distribution is non-normal. As noted on the figure, this data set has a large number 
of sample points (4077), a X 2 value of 245.06, and a significance level of 0 (to 8 places). Since 
the significance level is below 0.01, we say that y(l) is non-normal. 

The second bar chart for y(l) portrays the form of the probability -density function. Be- 
cause intervals of equal probability are used rather than intervals of equal amplitude (deviation 
from the mean), the top and bottom of the curve have been compressed relative to the middle 
portion. Yet, the skewness of the distribution is quite apparent in this chart, where it might 
have been overlooked in the left-hand chart. 
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Interval 


■ 

.5* 

.5* 

.2* 

.75 

.79 

. 74 

.83 

.83 

. ■■»? 

.87 

.91 

.91 

.95 

.OS 

.99 

.99 

1.03 

1.13 

1.06 

1.06 

1.10 

I . 10 

1.1* 

1.1* 

i. ie 

1.18 

1.22 

1.22 

1.26 

1.26 

1 .30 

1-30 

1.3* 

1 . 3 * 

1.38 

1.33 

I . -9 

1.59 


3808 

Points, 


Cm, A62-64, 7/29/66, 0602, Tl, R05, 700 Ft, 90 Deg, File 

20, TP 01, 05 Jun 67 



Frequency 

% 

Normal 

No. 

Normal 


.90 

.00 

1 

.1 

• • • 

.01 

.01 

3* 

27.6 


.01 

.31 

37 

3*. 1 


.02 

.32 

61 

6*. 9 


.03 

.03 

115 

112.8 


.05 

.35 

172 

176.9 


.36 

.37 

2** 

258.9 


.09 

.09 

3** 

3*1.9 


.19 

.11 

390 

*12.1 


.12 

.12 

*63 

*53.3 


.13 

.12 

*79 

*55.0 


. .11 

.11 

*19 

*16.6 


.10 

.09 

368 

3*6.5 


.07 

.07 

2*6 

265.9 


.05 

.05 

191 

185.2 


.03 

.03 

10* 

117.7 


.32 

.02 

77 

68.2 


.31 

.01 

32 

36.1 

• * 

.31 

.01 

29 

29.6 


3 

• DC 

0 

.1 

x 2 = ll .39, ' Significance Level = .72419039, 


.15 Degrees of Freedom 


FIGURE 66. EXAMPLE OF A GAUSSIAN SIGNAL DISTRIBUTION 
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Interval 
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.•-1 

.31 

. ) ! 

• C»7 

.13 

.13 

.14 

.19 

.2* 

.2* 

.3:. 

.37 

.36 

.36 

.*2 

.*2 

.*3 

.48 

.5* 

.5* 

.60 

.63 

.65 

.65 

.71 

.71 

.77 

.77 

.83 

.83 

.39 

.69 

.45 

.95 

1.01 

1.31 

1.16 

1.06 

1.12 

1.12 

l.ia 

1*18 

1.2* 

1.2* 

1.33 

1.30 

1. 36 

1.36 

L.*2 

l.*2 

1.65 

1.85 

• 

3550 Points 


Wht, C32-3A, 6/30/67, 0330, Tl. R04, 700 Ft, 180 Deg, File 16, TP 02, 07 Jun 67 
Frequency 




X = 2243.13, 


FIGURE 67. 


Significance Level = 0, 
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3 
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12.7 

.1 

.31 

3 

21.0 

J 
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> 
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.ul 

0 
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.02 

■3 
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.03 

12 

132. 1 
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.3* 

86 
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.38 

.35 
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.16 

.36 
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210.3 

.16 

.07 
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2*5.0 

.11 

.08 

39 3 

273.1 

.37 

.08 

2*6 

291.3 

.0* 

.09 

1** 

297.* 

.03 

.08 
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290.* 

.33 

p Jfi 
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271.* 

.33 

.57 

92 

2*2.7 

.33 

.06 
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.33 

.0* 
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.33 
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99.9 

.33 

.32 
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71.7 

.33 

.01 

9* 
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65 
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.31 
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53 
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0 
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0 
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CHANNEL 4 


Wht. C32-3. 6/30/66, 1400, Tl, R08, 2000 Ft, 360 Deg, File 15, TP 03, 15 Jun 67 


Interval Frequency % Normal No. Normal 

• .62 j .00 0 .1 

.<*2 .66 0 .01 0 21.5 

.65 .0* .00 .01 3 20.6 

.69 .?} • .00 .01 6 35.2 

.73 .77 •••• .01 .02 37 57.3 

.77 .81 .02 .03 70 87.7 

.81 .85 * .06 .06 129 126.5 

.85 .£9 .36 .05 201 171.5 

.89 .93 .09 .07 300 218.9 

.93 .97 .11 .01 355 262.9 

.97 l.jl .12 .09 373 297.1 

l.M lw6 .12 .10 376 316.0 

I.;< 1.98 - .07 .10 225 316.3 

l.CB 1.12 .06 .09 203 298.0 

1.12 1.16 .06 .08 191 266.1 

1.16 1.21. .05 .07 166 220.3 

1.2. 1.26 «•«»»•. .««.•••* .06 .35 161 172.9 

1.26 1.20 .06 .06 117 127. B 

1.26 1.32 .......... .03 .03 93 88.8 

1.32 1.36 ...... .02 .02 5B 58.1 

1.76 1,6; ...... .0? .01 53 35.8 

|.66 .31 .31 29 20.7 

1.66 1.69 ...... .3? .t'L 56 22.0 

1.65 « ^ .32 .0. 60 .1 

J 

3240 Points, X c 325.02, Significance Level - 0, 19 Degrees of Freedom 


FIGURE 68. EXAMPLE OF A SKEWED SIGNAL DISTRIBUTION 


CHANNELS,. Soy, A32-5, 6/29/66, 1308,-Tl. R02, 70.*0‘Ft, 90 Deg, TP 04, File S 05 Dec 67 

Interval Frequency % Normal No. Normal 

■ 22 3 .03 Cl .0 

.72 .36 0 .31 9 8.7 

.36 .38 0 •‘31 0 7.8 

. 33 . 6 '> « * .31 .01 10 13.3 

.60 .62 ••• .31 .32 16 21.3 

.62 .66 .36 .33 66 32.3 

.66 . 6 * .36 .36 81 66.2 

. 67 , .66 • .30 .05 102 62.6 

.40 . 5 - .39 .06 115 79 . B 

. 5 '; .52 .37 .08 35 96.6 

.57 .56 .............. .06 .09 57 110.3 

.56 .56 «... .36 .09 53 118.6 

.56 .58 .06 .39 72 120.9 

.58 ......................... ..*8 . 3 ? m 116.6 

.SI .62 .39 .78 116 106.0 

.62 .63 .13 .07 166 91.1 

.63 .65 .00 .36 97 76.1 

, j « .67 ,..•••«•«»•••»»■•» .36 .36 72 56.9 

.67 .69 c 03 .53 63 61.3 

.69 .71 .02 .52 23 2 B .3 

.71 .71 ... -31 .31 16 18.3 

.73 .75 • * -33 .51 5 11.2 

.75 .77 • • 3 .31 3 6.5 

.77 .91 ■ .03 .51 3 7.3 

. 9 1 • * . 3 « J j 6 .V 

2 

1275 Points, X = 257.01, Significance Level =0, 20 Degrees of Freedom 


FIGURE 69. EXAMPLE OF A BIMODAL DISTRIBUTION 
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The final two columns of figures shows the percentage and the number of points in each 
range. The expected percentages, assuming normality for the first and last intervals, is 0.005%, 
meaning that 20,000 points would be needed before one value in each of these intervals would be 
expected. The percentages for the remaining intervals should be close to 5% if y{l) are nor- 
mally distributed. 

The two bottom bar charts in Fig. 70 are for y(2) of the same data set. Since the significant 
level is greater than 0.01, we say that y(2) is normal. The bars on the left-hand chart have very 
nearly the same length. The right-hand chart shows the previously mentioned compression of 
the normal curve. 

A total of 54 sets of data were analyzed. From the 1969 data, we selected 38 fields (each 
with 10 channels of data), and from the 1966 data we selected 16 (each with 12 channels); a 
total of 616 variables were tested with results as tabulated in Tables 5 and 6. Of these 616, 

136, or 22% were non-normal, with y(l) and y(2) (the variables associated with the two largest 
eigenvalues) accounting for 71 or 52% of the total. The variation of the data in the direction of 
maximum eigenvalue tends to be the most non-normal. 

Since the total number of transformed variables that were classified as nonGaussian is 
136, the average per data set is ‘2.5. Ortly four sets were found that had all channels normal, 
while three had the maximum number found, six. Thus, it appears that the transformation we 
chose did tend to isolate the non-normal variations into a small number of variables. 

2 

The results of our X tests for non-normality of the signature distributions are summarized 

in Table 7. The 1% criterion used for the tests of the normal hypothesis is quite conservative. 

2 

Also, the X tests is not as powerful as more complicated tests would be. Consequently, we 
must conclude that a processor designed under the assumption that we have data with normal 
statistics is not optimum for data like those studies here, although we have not yet discussed 
the extent to which this nonoptimality affects recognition performance. 

Previous examination of the 700-ft data (Fig. 67) has shown that it has a regular or syste- 
matic variation [29]. As seen by the sensor, at this altitude, the fields subtend a large angle. 
Thus, there is a scan angle effect present which causes the data to have a mean value that is 
angle dependent. 

For the low-altitude data w T e tested, it is appropriate to use a filtering or preprocessing 
scheme to eliminate the scan angle effect. This approach, however, would not completely 
eliminate the problem of non-normality. One reason is that the scan-angle effects, in general 
would be reduced but could not be eliminated because of the noise in the data. In addition, we 
have evidence that the resulting preprocessed data would not necessarily be multivariate normal. 
For example, none of the high- altitude data we tested fully satisfied the simple hypothesis tests 
for normality we used, and these data subtended angles much smaller than the low- altitude data 
did. 113 



WILLOW RUN LABORATORIES 


TABLE 7. SUMMARY OF NON- NORMALITY TESTS ON SIGNATURES 


% Non- Normal Margins * 

% Non- Normal Margins Found in First 
Two Variables 

% Non- Normal Margins from 700- ft- Altitude 
Data 

% Non- Normal Margins from Higher .Altitude 
Data 


Untransformed 

(%) 

Transformed 

(%) 

65 

22 

17 

52 

78 

24 

■47 ' 

20 


2 

* Failed X test at 1% level of significance. 
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III. 4. SUMMARY OF EXPERIMENT RESULTS: CLASSIFICATION TESTS 

The final step in the testing procedure was to compare the ability of the two (normal and 
histogram) likelihood functions to classify correctly. We found the Type I errors by selecting 
a decision level and counting the percentage of points that were rejected using first one, then 
the other, decision rule. The Type II errors were found by direct calculation. Data points for 
the alternative were assumed to be located uniformly throughout a hyperrectangular parallele- 
piped, the dimensions of which were set so that 0.9995 of the volume of the Gaussian distribu- 
tions would be included. 

All of the comparisons were made with only three channels of data (after application of the 
eigenvalue transformation). The first three channels, corresponding to the largest eigenvalues 
were used for some of the comparisons; in other comparisons, the last three channels were 
used. When the last three channels were used, there was no noticeable difference between the 
results using the two decision rules. This was true for all of the data sets. 

When the first three channels were used, there was a difference in the performance of the 
two decision rules for a small fraction of the data sets. The Type II errors were compared for 
Type I errors of 0.1 to 0.3. Only three data sets were found where there was a noticeable im- 
provement in performance when the histogram likelihood functions were used. One set produced 
twice the Type II error rate for the same Type I probability, while the other two had one and 
one-half times the error rate. These ratios could be expected to be reduced if all of the chan- 
nels were used. Figure 71 shows the distribution of the first transformed variable for each of 
the three data sets. Note the obvious departure from normality. The importance of the two -to - 
one ratio of errors depends upon the actual separation between alternative hypotheses in the 
decision process. For two hypotheses with close separation between mean values, this decrease 
in performance when the Gaussian rule is used could be significant. 

Figure 72 shows the plot of Type I versus Type II errors for the data set with the largest 
difference in performance. The top curve shows the performance when we use Gaussian likeli- 
hood function in the decision rule while the bottom corresponds to the histogram rule. A 
Gaussian scale was used in preparing the graphs. If the data were normal, the normal decision 
rule would be approximately a straight line, with only minor deviations because a finite number 
(3550) of data points were used to find the Type I errors. 

Figure 73 shows the results of the use of the last three channels for the same data set as 
was used for Fig. 67. The two curves overlap. The curves in this figure are typical of all of 
the curves for which the last three channels of the data were used. The curves were also found 
to be typical of the large majority of results from the first three channels. Figure 61 shows one 
such graph. The two curves tend to overlap, especially for probabilities of Type I errors in the 
range of 0.1 to 0.3. Various decision levels were assumed for computation of the points on the 
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graph. The same decision level, when used with the two decision rules, produces different 
Type I and II errors. This can be seen by the location of the plotted points for the Type I error 
probability of approximately 0.3. For this example, at the same decision level, 8.17% of the 
points were classified differently by the two rules, although the two curves appear to overlap. 
Figure 74 a, b, c shows the complete set of operating characteristic curves for the 54 data sets 
when the 3 largest eigenvalue channels were used. 

For all of the curves, the Type n errors are correct within a multiplying factor. Remember 
that the uniform distribution was used to generate the Type II errors. Thus, the extent of 
values of the distribution is arbitrary, the choice having been made as a compromise between 
accuracy and computation time. The accuracy is limited by the number of data points avail- 
able for our curves, rather than by the extent of the uniform distributions. 

III. 5. CONCLUSIONS 

As a result of our tests, we have decided that the use of the normal likelihood function is 
justified for recognition processing of multispectral scanner data. This function is much 
quicker to generate and use than the histogram function. Also, the improvement in performance 
that would result from the use of a histogram likelihood function and all of the channels, rather 
than only three channels, is not believed to be significant. If one were to include the effects of 
interfield variations, even smaller differences probably would result. More promising 
approaches are: (1) to preprocess the data to compensate for scan angle and similar systematic 
effects (Refs. 2, 3, and 4, Appendix V) and then use the normal assumption in the decision rule; 
or (2) to compensate by changing the decision rule parameters. By using the first approach, we 
would expect transformed data which had distributions that are more nearly normal and, at the 
same time, which had reduced variances. 
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(a) 

FIGURE 74. EXAMPLES OF OPERATING CHARACTERISTIC CURVES WHEN 
FIRST THREE VARIABLES ARE USED 
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(b) 

FIGURE 74. EXAMPLES OF OPERATING CHARACTERISTIC CURVES WHEN 
FIRST THREE VARIABLES ARE USED (Continued) 
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(c) 

FIGURE 74. EXAMPLES OF OPERATING CHARACTERISTIC CURVES WHEN 
FIRST THREE VARIABLES ARE USED (Concluded) 
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IMPORTANCE OF ATMOSPHERIC SCATTERING 
IN REMOTE SENSING* 


R. E. Turner 
W. A. Malila 
R. F. Nalepka 


* Reprinted from Proceedings of the Seventh International Symposium on Remote Sensing of 
Environment, Willow Run Laboratories of the Institute of Science and Technology, The University 
of Michigan, Ann Arbor, 1971, pp. 1651-1697. 
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IMPORTANCE OF ATM3SPHERIC SCATTERING IN REMOTE SENSING 


OH EVERYTHING YOU* VE ALWAYS WANTED TO KNOW ABOUT 

ATMOSPHERIC SCATTERING BUI 1 WERE AFRAID TO ASK * 

R.E. Turner 
W.A. Mali la 
R.P. Nalepka** 

Willow Run Laboratories 
Institute of Science and Technology 
The University of Michigan 
Arm Arbor, Michigan 

ABSTRACT 

One’s ability in remote sensing to discriminate between 
target and background materials on Earth’s surface is diminished 
by the presence of atmospheric haze which scatters and/or 
absorbs parts of the radiation. The extent to which the 
atmosphere affects the process of discrimination depends on 
the specific condition of the atmosphere existing at the time 
of observation as well as on the discrimination technique 
enployed. Remote sensing operations in environmental and 
natural science applications would benefit from a better 
and more quantitative understanding of the effects of the 
atmosphere on discrimination.. 

This paper asks and answers questions that are of interest 
to users of remote sensor data and presents calculations 
designed to inprove one's understanding of radiation transfer 
in Earth’s atmosphere. The spectral region considered lies 
between 0.4 and 3 pm. A recently developed radiative transfer 
model was used in corrputing irradiance , path radiance, sky 
radiance, transmittance, and contrast transmittance in a 
cloudless atmosphere for a variety of sun angles, viewing 
(nadir) angles, altitudes, surface reflectances, and atmos- 
pheric haze conditions. The results are: presented para- 
metrically. One interesting result is the effect which 
neighboring materials have on the spectral character of a 
target as a result of aerosol scattering by the haze. 

The model described was developed using radiative transfer 
theory as applied to a plane-parallel, homogeneous, aerosol- 
filled atmosphere with a surface, the reflectance Qf which is 
assumed to be perfectly diffuse. The haze content of the 
atmosphere is designated by the horizontal visual range at 
sea level. 


1. WHY SHOULD A USER OF REMOTE SENSOR DATA BE CONCERNED ABOUT ATMOSPHERIC SCATTERING? 

Maury problems associated with target recognition and image interpretation in remote sensing are 
caused by the presence of Earth’s atmosphere. One problem is clearly evident in photographs of the 
Earth's surface taken from orbiting spacecraft 1 . The general bluish tint in many of these photo- 

* The work reported in this paper was supported by NASA under Contract NAS9-9784. 

** The authors are enployees of: Infrared and Optics Laboratory, Willow Run Laboratories, 

Institute of Science and Technology, University of Michigan. 
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graphs arises from radiation that has been scattered one or more times by the molecules and 
particulate material which conpose the atmosphere. This scattered radiation adds spectral 
variations to the received signals which are not truly representative of surface features. 

Since most remote sensing operations involve observation paths within or through Earth’s 
atmosphere, each investigator must consider the effects of the atmosphere on the real at ion signals 
produced by the sensors and on his use and interpretation of the data. If the atmospheric effects 
on the signals are serious enough, he should then consider ways and means of reducing or overcoming 
these effects. 

Tire re are several deleterious effects of atmospheric haze on remote sensor data, whether in 
the form of imagery or of electrical signals processed on computers. First, base reduces the 
contrast between adjacent surface features on imagery and makes the interpretation of that imagery 
more difficult. Figure 1 illustrates the way that contrasts are reduced on Images obtained in 
two spectral bands at 5000 - ft. altitude from those at 1000-ft. altitude. 

Second, as a result of scattering in the atmosphere, the spectral distribution of radiation 
received by a detector Is different from that which originates at the ground surface, i.e., the 
atmospheric effects are more p renounced at some wavelengths than at others. For example, in 
Figure 1, the effect of the haze is much greater for the wavelength band, 0.55 vm to 0.58 urn, than 
for the band at longer wavelengths, 0.30 pm to 1.0 pm. It is this spectral dependence of scattering 
which accounts for the blueness in color photography and color conposite imagery taken from high- 
altitude aircraft and spacecraft. 

Finally, a tiro spheric effects can reduce- the amount and quality of information that can be 
extracted from the data by computer processing. That Is, the first two effects discussed with 
re ga rd to imagery and its interpretation also are important in the recogiition processing of 
multispectral remote sensor data. The effects of atmospheric haze can reduce the recognition 
computer’s ability to discriminate between surface materials that resemble each other, and can 
limit the extent of the area over which satisfactory recognition results can be obtained easily. 

2. HOW DOES THE ATMOSPHERE AFFECT REMOTE SENSOR SIGNALS? 

The atmosphere affects visible and near-infrared remote sensor signals in three ways. First, 
it modifies the spectral' and spatial -distributi onsr'of the -radiation -incident on the surfaces being 
observed. Second, it attenuates the radiation that passes through it to the sensor after being 
reflected by the surface. Third, it adds an extraneous component of scattered radiation, called 
path radiance, to the transmitted component. (See Figure 2). 

Hie total radiance received by a sensor at altitude, h, with nadir view angle, e, and 
azimuthal angle, 4 , can be expressed as: 

L(h,e,<f>) = I^(h, 0 , 4 >) + L (h,6,4) (1) 

where: Lp(h,e,4) is the transmitted radiance and Lp(h,6,4) is the path radiance. Althou^i not 

noted explicitly, these radiances also depend cn wavelength, solar zenith and azimuthal angles, the 
amount and distribution of haze particles arid molecules in the atmosphere, and the surface 
reflectance. 

The transmitted radiance is the product of the radiance emanating from the surface, here 
called the Intrinsic radiance, and the transmittance of the atmosphere, T(h,e). That Is: 

L*p(h, 6, 4 ) ° Lj( 6 , 4 )T 1 (h, 0 ) 

Here, L_( 6 , 4 ), the intrinsic radiance of the surface, is given by: 

2 tt 1 

u'p "(p>4,-u'»4 v )L(C , >~p "d4' 


where: L(0,~ii',4'‘) is the radiance incident on the surface, p"(u, 4 ,-:T, 4 ") is the bidirectional 
reflectance of the surface, and p = cose. The primed variables indicate a downward direction, and 
the unprimed quantities indicate an upward direction. 


Lj-(0,4) - | ^ 


( 2 ) 

(3) 
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The quantity T(h,0) is the transmittance and is given by: 

— ( T -T }/u 

T(h,0) * e ° (*0 


where t and n axe the optical depth at altitude h and the total optical depth , respectively, 
for the atmosphere , and can be determined for any wavelength and atmospheric condition. 

3. WHY IS A DETAILED RADIATIVE TRANSFER NDEEL IMPORTANT FOR FOOTE SENSING APPLICATIONS? 

One can leain about atmospheric effects from either experimental measurements or theoretical 
calculations with a radiative transfer model, or both. A large number of variables must be 
considered, and an extensive measurement or calculation program would be required tc get a full 
understanding of the severity of the effects under a wide variety of conditions^- 

The use of radiative transfer model calculations has several advantages over the use of 
experiments: lower cost, more flexibility, better control, and a freedom to study effects 

parametrically . To keep from drawing erroneous or misleading conclusions from the calculations, 
however, the model must be a detailed and accurate one and most be verified by comparison with 
experimental data. 

A knowledge of atmospheric effects is useful in all phases of remote sensing operations. In 
the data collection phase, it can be useful in the specification and design of sensors and in 
mission planning. In the recognition processing phase, it can be used to remove signal variations 
introduced by the atmosphere thereby making large-area surveys more feasible. That is, this 
knowledge may be used directly to specify parameters of preprocessing transformations or may be 
used in their development through the simulation of sensor data exhibiting various atmospheric 
effects. Finally, in the interpretation phase, it can help the investigator better understand 
features in imagery and recognition maps and extract quantitative estimates of surface radiation 
characteristics . 

HAVEN'T SCATTERING PHENOMENA IN EARTH'S ATMOSPHERE BEEN UNDERSTOOD FOR A LONG TINE? 

We are all aware of many beautiful phenomena whiqh ..occur. in the. sky from time to time, such as 
the red and orange sunsets, the green flash, the aurora borealis, and rainbows, but by far the most 
obvious phenomenon we notice is the blue sky. Why is the sky blue? Lord Rayleigh 2 found that 
the blue color of the sky could be explained by considering the single scattering of li^it by 
scattering centers in the atmosphere. The scattering centers, usually assumed to be molecules, axe 
actually molec ular density fluctuations, the sizes of which are much smaller than the wavelength 
of the radiation. If incident radiation is scattered once by an inhomogeneity, then, for Rayleigh 
scattering, the intensity of the scattered radiation varies as: 

I(Rayleigh) * X **(1 + cos^e) (5) 

where X is the wavelength of light and 0 is the scattering angle. Thus, we see from (5) that short 
wavelength radiation is scattered much more than long wavelength radiation. 

Actual measurements of the sky radiation reveal, however, large deviations from Rayleign's 
theory. Tne spectral character can differ strongly from the inverse fourth power relationship 
given by (5), and the degree of polarization is also different from that predicted by the theory. 

The discrepancies between theory and experiment can be explained by the following: (1) 

Rayleigfi considered cnly single scattering whereas multiple scattering is also important, and (2) 
the real atmosphere contains particulate matter which can scatter radiation in a manner quite 
different from that given by (5). 

Significant advances were made in radiative transfer theory by Schuster in' 1906 and by 
Schwarzs child 3 in 191^ but the exact solution of the multiple scattering problem in plane -parallel, 
homogeneous, conservative atmospheres which have a Rayleigh scattering lav/ was found by 
Chandrasekhar . 4 Using Chandrasekhar's theory, Coulson et al., 5 computed the sky radiation for a 

Rayleigh- type atmosphere with multiple scattering and polarization. The optical thickness 
considered, however, was small and cnly recently^ have similar calculations been performed for 
atmospheres of large optical thicknesses. 
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In recent years , the theory of radiative transfer has been extended to include absorption, 7 
inhonugeneitles , 8 and internal field calculations. 9 In all these Investigations, only Rayleigh- 
type scattering or scattering according to a slightly anisotropic distribution was used, but, 
even with these developments, theory and experiment did not agree, especially for hazy atmospheres 
in which scattering by particulate material is predominant. 

Thus, we can say that scattering phenomena are now understood for Rayleigh-type atmospheres 
under* special conditions. However, for a realistic atmosphere, characteristic of Earth’s 
atmosphere , the particulate corrponent, typical of haze, must be taken into account. 

5. WHY IS THE RADIATIVE TRANSFER PROBLEM FOR HAZY ATMOSPHERES SO DIFFICULT TO SOLVE? 

As we have just noted, Rayleigh scattering refers to the scattering of electromagnetic, radia- 
tion by objects much smaller- than the wavelength of the incident radiation. In the case of Earth’s 
atmosphere , the scattering centers are regions where the molecular density fluctuates and, hence, 
Rayleigh scattering is also referred to as molecular scattering. For the theoretical case of a 
pure Rayleigh- type atmosphere, the visibility is very great, about 300 kilometers or more. The 
actual atmosphere, however, consists of a turbulent mixture of gases and particulate matter, a 
combination which is called an aerosol. Typical aerosol particles have radii in the range 

10*"^ pm to 10° pm but can vary from 10“^ pm to 10 2 pm. For radiation in the visible region of the 

spectrum, it is a particle in the range lCf'*’ pm to 10^ pm (haze particle) which scatters radiation 
most strongly. 

If a simple shape, e.g. , a sphere, and an index of refraction are assumed for a haze particle, 
then classical electromagnetic theory can be applied to determine the cross-section and scattering 
properties of the particle. For the special case of a homogeneous sphere irradiated by a plane 
electromagnetic wave in one direction, standard elect romagne tic scattering theory (Mie scattering) 
is applicable. The results of calcinations show that the cross section for a particle of radius 
r is strongly dependent upon the dimensionless ratio 2^r/A, where X Is the wavelength of the 
radiation. An actual aerosol distribution is, however, a polydispersion, i.e., a collection of 
particles of varying sizes and therefore an effective cross section must be found for the 
distribution. 

Althougn we need the scattering, absorption, and total cross sections of the polydispersion, 
it is also necessary that we have knowledge of- the angular distribution of the radiation scattered 
by a particle in order to solve the transfer equation. The quantity which describes the angular 
variation is called the single-scattering phase function, a dimensionless function which represents 
the fraction of the radiation scattered through the angle 9. Examples of these functions in terms 
of wavelength are shown in Figure 3 for a typical haze condition. 10 The most striking characteristic 
associated with aerosol particles is the very high degree of anisotropy in the scattered radiation 
field. Typically, the amount of radiation scattered into the forward hemisphere is v 95$ with most 
of that being within the first four degrees. 

The use of such anisotropic functions in radiative transfer problems usually involves extensive 
analysis and computation. To illustrate, consider the expansion of the phase function in spherical 
hareonics, i.e., 

N 

p( v)= l (21 + 1) A £ p £ (w) (6) 

£=0 


where p(u) is the actual poly disperse phase function, u is the cosine of the scattering angle e, 
P £ (y) is a Legendre polynomial, and A £ are constants given by: 



p(u) P £ (c)dii 


(7) 


Using values of the actual phase function p(u) from Deirmendjian 10 we calculated the 
coefficients A £ numerically and inserted them into Equation 6. The results are shown in Figure 


where only the backward hemisphere is depicted since there is almost perfect agreement for the 
forward hemisphere. One must use some 36 or 38 terms to approximate the actual phase function, 
while only two teres are needed to represent the sinple Rayleigh function. 
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Thus, there are two factors which serve to complicate the problem of radiative transfer in 
hazy atmospheres; ( 1 ) the uncertainty in the specification of the parameters needed to define an 
atimspheric condition* and ( 2 ) the mathematical conplexities involved in solving the transfer 
equation which has a highly anisotropic phase function. 

6 . WHICH PARAkETERS ARE USEFUL FOR DESCRIBING THE STATE OF THE ATMOSPHERE? 

Earth's atmosphere consists of semi-permanent gases such as nitrogen, oxygen, and argon and 
highly variable corrponsnts such as o 2 one, water vapor, and aerosols. The latter take on many forms; 
water droplets, sinoke particles, dust, fumes, etc.. Fbr the purposes of remote sensing in the 
ultraviolet or infrared spectral regions, we need information on tlie atmospheric content of gases 
such as water vapor, carbon dioxide, and ozone, since these gas.es absorb significant amounts of 
radiation. Radiosonde data can supply the pressure, tenperature , density, and relative humidity 
data and independent measurements provide data on gaseous conposition. 

In the visible spectral region, very little absorption occurs by gases or aerosol particles 
so that the primary attenuating mechanism is scattering. Since the tropospheric aerosol is highly 
variable in nature , it is desirable to have some information concerning the amount of particulate 
material present at times of interest. 

One useful parameter for the visible region is the sea level horizontal visual range, V. 
Assuming a 2% relative contrast between a black object and the sky background we have: 


= e" KX = 0 . 02 ' . ( 8 ) 

hi 

where k is the extinction coefficient (km” 1 ) and X is a standard distance. Defining the variable 
distance X to be equal to V for this particular relative contrast we have: 


V = 3^912 ( 9 ) 

K 

a simple relation which enables us to define the extinction coefficient < for a haze according to 
the easily measured parameter visual range (km). ... 

Other parameters can of course be specified to cte fine the state of the atmosphere more 
accurately, but from the practical point of view used in remote sensing applications , the horizontal 
visual range will usually suffice. 

7. DOES A RADIATIVE TRANSFER MODEL EXIST FOR HAZY ATMOSPHERES? 

The rathematical conplexities in a radiative transfer model arises from attempts to solve the 
integro-di f ferential equation of radiative transfer. The equation to be solved for the description 
of radiative transfer through a plane-parallel atmosphere is the following: 




dL 

dr 


L(t,v,4>) 



| p(t,p, <j>,vi%4T)L( T,y',<})')dvi'd$' 
'-1 


( 10 ) 


L(t,p,$) is the scalar spectral radiance at an optical depth, t, where we define the dimensionless 
quantity t to be: 

t = f ic(z)dz (11) 

'h 


k(z) being the extinction coefficient at altitude z, and h being the base altitude. The quantity 
w is the single-scattering albedo, defined as the ratio of the scattering coefficient to the 

extinction coefficient and is usually very close to unity for the visible spectral region. The 
function p(-r,y,<ji,u% 4 >') is the single-scattering phase function which describes the angular 
distribution of radiation scattered through the angle x, where 


oosx = y v' + Al-rt) ( l-u"^) cos ( 41 -$') 


( 12 ) 
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There are three general methods used to deteraine the spectral radiance L(r,u, $) . These are: 

(1) the classical appreacb, in which the radiance function L( r,ij, $) and the phase function are 

treated as continuously varying functions of o(=cos \i) and (2) discrete space theory, where 
the transfer equation is solved for the radiance at specific angles 6. , and (3) the statistical 

technique, i.e. , one which makes use of sampling theory such as. the Monte Carlo method. The decision 
as to which metJiod should be used depends somewhat on the personal preferences of the investigator 
and also on the accuracy desired. 

In the classical approach, one can expand the radiance functions L( t, u, e ) in spherical 
harmonics and, using the orthogonality properties, convert the transfer equation into an eigenvalue 
problem. Although the analytical formulation is complicated because of the large number of terms 
needed, this is one of the approaches we are now examining at The University of Michigan. 

There are several other methods in the classical approach, one being the conversion of the 
transfer equation into an integral equation which can be reduced to a Neumann series. Each renter 
of the series corresponds to a higher degree of scattering. The main reason why this method is 
not used often is that convergence of the series is very slow under certain conditions. 

Discrete space theory is more readily adaptable to corrputer methods but is quite time consuming, 
especially when inhomogeneous at respheres are considered. An example of this method is the 
invariant imbedding technique in which the linear equation 10 with a two-point boundary condition 
is replaced by a nonlinear equation with an initial condition. Various thin layers of an atmos- 
pheric medium are added successively until the total thickness is obtained. Thus, the internal 
field can be obtained for any atmospheric thickness, but so far the rethod has not been used with 
highly anisotropic phase functions that are characteristic of a hazy atmosphere. In this case too, 
an excessive amount of conputer time is necessary in order to simulate real atmospheres with largp 
optical depths. 

The Monte Carlo technique has enjoyed great success in recent years and has as its most 
desirable feature, versatility. That is, one can simulate many unusual atmospheric and surface 
conditions easily and hence eliminate the enormous analytical complexities using the other ret hods. 
Unfortunately, much corrputer time is needed to achieve reasonable accuracy. 

Therefore, we can say that, in spite of all the research over the past thirty years or so, no 
detailed, conprehensive radiative transfer model has been developed which can be applied directly 
to the problems in remote sensing. In order to rectify this situation, we have developed a 
sirrplified model at The Uni' varsity of Michigan and are presently extending the development to 
include a large variety of atmospheric and terrestrial conditions. 

8. WHAT ARE THE CAPABILITIES OF THE RADIATIVE TRANSFER MODEL DEVELOPED AT r IHE UNIVERSITY OF 
MICHIGAN? 


For immediate application to problems of discrimination in remote sensing, a user needs a 
sinplified, workable model which can be applied to a wide range of atmospheric and surface conditions 
yjith a minimum of computational effort. Our model is especially equipped to deal with plane- 
parallel, homogeneous,* hazy atmospheres, under a variety of situations. 


Realizing that aerosol scattering occurs predominantly in the forward direction, we can 
approximate the phase function by a sum of delta functions, i.e.. 


p(p,<f , ,iJ''i<t ,/ ‘) - F<S(vi-u')<S(tf-<j>') + B<5(p+ii06(Tr+'f> 


(13) 


where F « Ann and B = 4n(l-n), n being the fraction of the radiation which is scattered into the 
forward hemisphere. Inserting Equation 13 into the following radiative transfer equation: 

E ( t ) 

w ^7 = L(i,u,<Ji) - 7J7 | * { P^,^,^, 4>")L(T,p', | |>'‘)dM''d , }>'- — ^ — p(u,<|>,-u 0 ,ij> 0 ) (1*0 

J 0 - 1 

allows us to separate the equation into two simple differential equations. E s (t) Is the attenuated 
solar spectral irradiance, and cose^) , <f> o are the coordinates of the sun. Assuming a surface 


B By homogeneous, we mean that the proportions of constituents are constant throughout all altitudes, 
even though the actual density varies drastically with altitude. We are presently modifying our 
model to include inhomogeneous atmospheres. 
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reflectance cf zero, we solve the equations for the spectral irradiances In the upward and downward 
hemispheres at any optical depth x. Then, assuming a perfectly diffuse surface reflectance, we 
find the radiation field resulting from surface reflection and contoine these irradiances to deter- 
mine a "source" radiance, i.e.. 


L(t,u, 4 >) “ ~ E^(x ) 6 (u-p^) 6 ( ir+ 4 ^- 4 )+ x)6(p4-p^)6(4— 4 ^) 


E^(t) 4- E" ( x ) 


2u 


(15) 


where here L(x,y, 4) represents the source radiance and is to be inserted into the integral of 
Equation 1*1. Equation 14 can then be solved to determine the spectral sky radiance and the spectral 
path radiance at any point in the atmosphere. The boundary conditions for the spectral radiance 
are : 


L(0,-u,4) = 0 

r 2 tx rl 

L ( T n >y>$) = 

u '0 J G 0 


(16) 

(17) 


where p'(u,4,-u',4') is the bidirectional reflectance of the surface and ,$') is the total 

(direct plus diffuse) spectral radiance at the surface. It should be noted that for the case of a 
perfectly diffuse (Lambertian) surface Equation 17 reduces to 

= f S_(x o ) (18) 

where p is the hemispherical reflectance, and & (t) is the total downward lr radiance at the 
surface. 0 

Finding the radiance in terns of optical, depth x instead of altitude h allows one to modify 
the state of the atmosphere without affecting the radiative transfer calculations. In the model 
currently being used at The University of Michigan, we utilize the optical depth-altitude -wavelength 
relationship as determined by Elterman 11 for visual ranges from 2km to 23km. The actual scattering 
phase functions for a haze, taken from Deinnencljian, 10 are then used in the solution of the radia- 
tive transfer equation* Assuiriing some surface reflectance, we can then de ter ra ! h p the spectral 
radiances in the atmosphere for a variety of conditions. An outline of the general capabilities 
of our current radiative transfer program at The University of Michigan is illustrated in Figure 5. 12 

The ultimate test for the validity of any theoretical model lies in its agreement with experi- 
ment. jji figures 6 and 7, a comparison is made between calculations made with our model and 
experimental data obtained by Ivanov.^ 3 Although the surface conditions were not known exactly, 
it was known that the surface was one of vegetation and, therefore, had a reflectance of about 20%. 

We are now continuing the development of our radiative transfer model to include a variation 
in atmospheric parameters.. One long-term variation is illustrated In Figure 8 in which results of 
our model were presented for (1) data reported in 196^ 1 4 but gathered prior to the Tit Agung 
volcanic eruption in 1963 and (2) data gathered several years later (1968). 15 The volcanic 
eruption added considerable quantities of dust to the stratosphere and, hence, a reduction in 
irradiance occurred at the surface. Thus, there is a reduction in the. upward scattered radiation 
and the path radiance is less. 

9. K)W DOES THE TRANSMITTANCE DEPEND ON THE ATMOSPHERIC AND OBSERVATION CONDITIONS? 

The transmittance function was described in Equation It is a function of the various 
parameters that describe the path conditions a s well as of wavelength. Figure 9 presents plots of 
spectral transmittance for several visual ranges. Note that transmittance is smallest for the 
shortest wavelengths. A more detailed presentation of the dependence of spectral transmittance on 
visual range is given in Figure 10 for a wavelength of O.55 M m. 

Both altitude and scan angle effect the length of the observation path and consequently affect 
the transmittance. Figure 11 illustrates the altitude dependence of spectral transmittance and 
Figure 12 Illustrates its dependence on nadir scan angle for several visual ranges. 
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10. HOW DOES IRRADIANCE DEPEND ON THE ATFDSPHERIC AND OBSERVATION CONDITIONS? 

There are two components of downward Ir radiance within Earth's atmosphere, the direct solar 
corrponent and the diffuse downward component. As shown in Figure 13, the diffuse spectral irradiance 
component, E_( x ) , depends very strongly on the vertical distribution of haze particles and on the 
altitude of the observer. The total downward spectral irradiance, E_(r), however, exhibits a nxrch 
smaller dependence on altitude. 

As the visual range decreases, the diffuse corrponent at any altitude increases as shown in 
Figure 1*1. Tie curves in this figure suggest that measurements of the diffuse irradiance at two or 
irons low altitudes could be used to estimate the visual range or haze distribution in the atmosphere. 
Plans are being made to incorporate a diffuse irradiance sensor in the Michigan multispectral 
aircraft. 

It is well known that scattering is more pronounced at shorter wavelengths than at longer ones. 
Figure 15 illustrates the spectral distribution of Irradiance for two visual ranges at an altitude 
of 1 km. 

Finally, Figure 16 presents ratios of irradiance at given altitudes to that at sea level for 
three visual ranges . This plot is of interest because an irradiance sensor aboard the multispectral 
aircraft is frequently used in extimting the irradiance at the ground surface. At an altitude of 
3 km, there is a difference of about % under reasonable clear conditions (V=23 km), a 7 % difference 
under moderate haze conditions (V=8 km), and a 1*1$ difference under dense haze conditions (V=2 km). 

11. HOW DOES PATH RADIANCE DEPEND ON THE ATMOSPHERIC AND OBSERVATION CONDITIONS? 

The path radiance term is inportant because it is an additive term that tends to mask the 
transmitted scene radiances. As shown in Figure If, path radiance is strongly dependent on ids hare 
condition of the atmosphere and on the altitude of the observer. 

In a seaming system, the observation path lengthens as the scan angle moves from the nadir. 

This increased path length results in substantial increases in path radiance as can be seen in 
Figure 18. For many situations, the path radiance can be an appreciable part of the received signal. 
The peak value seen on the V-2 km curve is at 0 = 60° , the antisolar angle. That is, the observer 
is looking directly away from the sun. Figure 19 illustrates the change in shape that occurs when 
other"* solar zenith angles are considered. 

A fourth graph. Figure 20, clearly illustrates the interesting fact that the path radiance 
depends substantially on the albedo (diffuse reflectance) of tire ground surface. 

12 . HOW CAN THE RADIATIVE TRANSFER MODEL HELP ONE UNDERSTAND THE EFFECTS OF THE ATM3SPHERE ON 
IMAGE CONTRASTS? 

Investigators have often used a quantity called contrast transmittance to describe the manner 
in which the atmosphere reduces contrasts in Imagery. The contrast between a target material and 
its background is usually defined to be: 


L t 

C(h) * ~ 

where L t (h) is the target radiance at altitude h, and I^(h) is the corresponding background radiance. 
Recall that, from Equation 2, we have 

L^(h) = Lj. (O)T(h) + Lp(h) 



and 

L t (h) = Lj (O)T(h) + I^(h) . 

Contrast transmittance is defined as the ratio of the contrast, C(h), at any altitude to that 
at sea level, C(0). Thus, 
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U (O)T(h) 

„ v C(h) _ x b 

V h) ■ c(0l “ L^oyrTh) +■ L p (h) 


which reduces to; 

V h > = 


whe re 


p = background albedo (diffuse reflectance) and E_ = total downward lrradiance. note wish, 

can be defined either spectrally or over a band of wavelengths, and also that it is independent of 

the target characteristics. 

Calculations of contrast transmittance were made using quantities generated with our radiative 
transfer model for a variety of atmospheric and observation conditions. One normally sees plots 
of contrast transmittance versus altitude for various atmospheric conditions (visual ranges) as 
shown in Figure 21 for a wavelength of 0.55 vtn. For a fixed altitude of 3 km, the spectral 
dependence of contrast transmittance in Figure 22 reveals the influence of the reflectance spectrum 
of the green vegetation background surface considered. 

A useful presentation of contrast transmittance versus visual range is made in Figure 23 for 
several altitudes. Mote there is little difference in contrast transmittance at a given altitude 
for visual ranges between 15 and 70 km; the dependence on visual range is greatest at short visual 
ranges. 

Another interesting graph. Figure 24, shows the dependence of contrast transmittance on the 
solar zenith angle. There is a surprisingly strong dependence caused by the increase in path 
radiance that occurs as the solar zenith angle approaches 0° and the decrease in surface irradiance 
that occurs for low sun angles. 

Up to this point, all contrast transmittance plots have been made for a green vegetation back- 
ground. Figure 25, illustrates the dependence of contrast transmittance on the diffuse background 
reflectance (i.e., albedo). This type of plot can be misleading since T p becomes zero for a 
background reflectance of zero while, on the other hand, C(0), the contrast being transmitted, 
becomes infinite according to its definition. This points out the fact that the investigator must 
consider both the inherent contrast in the scene and the contrast transmittance in his interpreta- 
tion and analysis. 

13. WHAT IS THE PRACTICAL UTILITY OF OUR RADIATIVE TRANSFER 1*0 EEL? 

As nen Honed earlier, a radiative transfer model which properly accounts for the characteristics 
of a real atmosphere will be an extremely useful tool for engineers and scientists involved in the 
many aspects of remote sensing, from system design to data analysis and interpretation. 

A thorough, understanding of the effects of the atmosphere on the radiation available for 
ejection under the wide range of atmospheric eenditions that may exist will permit the design of 
new and better remote sensing systems. For data being gathered with presently existing systems, 
such a model will provide a means for identifying and possibly reducing the deleterious effects 
of atmospheric scatter. Also, hypotheses which have been offered to explain effects noticed in the 
analysis of remotely sensed data can be verified or rejected. An example of the latter use of the 
radiative transfer model is described briefly in the next few paragraphs. 

An interesting and perhaps significant effect was noticed while processing and analyzing a set 
of airborne multispectral data which was gathered under hazy atmospheric eenditions*. It was 


LTThT 


:o) = — £ (o) 




» t m 


* A more detailed discussion than that given here can be found in 0-6]. 
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determined that, under certain circumstances, the radiance spectrum being sensed from small objects 
at the receiver assumed the spectral characteristics of surrounding objects, that is, objects not 
in the receiver's instantaneous field of view ( IFGV ) . The resulting signals therefore were not 
representative of the objects appearing within the IFOV. It is clear that this situation, whenever 
it might arise, would seriously reduce the accuracy of automatic recognition processing since the . 
signature of each class of objects would vary depending on the background in which it was found. 

As a part of the investigation, an attempt was made to extract a quantitative measure of the 
spectral path radiance existing at the time of data collection. This was accomplished by utilizing 
data gathered over a set of spectrally neutral (gray) reflectance panels which were placed in the 
scene. Having measured the reflectance of the panels in the laboratory and extracted the voltage 
generated when viewing them from the air, it was possible to construct a voltage-re fleet ance 
transfer curve in each of the airborne scanner spectral bands. An extrapolation of each of those 
curves was carried out and the voltage Intercept corresponding to zero reflectance was determined. 
Those voltages wen? interpreted as resulting from viewing an object having zero reflectance, 
therefore, the radiance associated with the zero reflectance voltage in each spectral band did not 
emanate from objects on the ground which were within the receiver IFOV. 

Using sources of standard radiance in the scanner, it was possible to calibrate voltages in 
terms of radiance and to convert the zero reflectance voltages to path radiance values. The path 
radiance spectrum so computed is illustrated in Figure 26 where the vertical lines indicate the 
range of uncertainty due to scanner system noise. The shape of the distribution illustrated, 
which has a peak in the green portion and a maximum in the near- infrared portion of the spectrum 
was initially somewhat of a surprise. It was expected that the primary source of path radiance 
would result from the direct scattering of solar radiation into the receiver. However, if this 
were the case, the path radiance spectrum would have exhibited a maximum at considerably shorter 
wavelengths in the region where the solar spectrum peaks. This was certainly not true here and in 
fact the resulting path radiance spectrum looked suspiciously like a reflectance spectrum for green 
vegetation. 

Since the reflectance panels were placed in a grassy area, several hypotheses which could 
explain the conputed radiance spectrum were explored in detail. It was finally determined that a 
single hypothesis best -explained the conputed radiance spectrum; -namely that radiation from nearby 
objects outside the receiver's IFOV (in this case grass) was being scattered into the receiver. 

As a check on this hypothesis, our radiative transfer model was employed to calculate the 
expected path radiance spectrum under low visibility conditions when viewing a spectrally neutral 
object in a bacKground of green vegetation. Although the atmospheric conditions existing during 
the time of data collection were rapidly varying and could not be accurately specified it was 
believed that the plausibility of the hypothesis could still be determined. The results of the 
rredei calculations are illustrated in Figure 27 for a 3-km altitude. Here again, the path radiance 
spectra exhibited local maxima at 0.55 and 0.80 pm thereby lending credence to the hypothesis. 

An additional effect was noted during the analysis of the multispectral scanner data gathered 
on the same day. These data were gathered over an agricultural area at several flight altitudes 
to enable an enpirical determination of the effect of flight altitude on the radiance spectra 
received when viewing objects of interest to the agricultural community. The ire an radiance spectra 
received at flight altitudes of 1000 and 5000 feet for a field of bare soil and a field of soybeans 
are depicted in Figures 28 and 29. The interesting effect here is that while the entire 1000 ft. - 
radiance spectrum for bare soil lies below the 5000 ft. spectrum this is not the case for the 
radiance spectra received when viewing the soybean field. For the soybean field, in the spectral 
bands centered at 0.71 and 0.80 ym, the spectral radiance received at 1000 ft. exceeds that received 
at 5000 ft. 

Although, the existence of such an effect was satisfactorily explained in a qualitative manner, 
the specific results could have been predicted through the use of the radiative transfer model. 
Figure 30 illustrates the spectral radiance conputed using the radiative transfer model for viewing 
a small object having a reflectance of 8 % surrounded by a background of green vegetation and an 
extended background of green vegetation. The results seen here are similar to those evident in 
Figures 28 and 29. 

Hopefully, the examples illustrated above are sufficient to provide an indication of the 
practical utility of the radiative transfer model. The real advantage of such a model, however, 
is that it permits investigations of atmospheric effects not only for conditions that do exist 
during the collection of any particular data set but for conditions that, in the mind of the 
Investigator, might exist at some future time. Because of this, individual parameters may be 
varied and the effect of their variation determined. From this Information Insight can be gained 
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into the results to be expected with present renote sensing systems and lirproved systems can be 
desigied for tire future . 

In order’ to further Illustrate the effects of certain parameters , we include in Figures 31 
through 33 examples of results calculated of our radiative transfer model. These results demon- 
strate the variation of transmitted and apparent radiance with altitude, wavelength, arid visual 
rang? for a variety of conditions. 

It, WHAT ARE THE FUTURE DEVEU>PMENT5 OF THE RADIATIVE TRANSFER fODEL? 

Having accounted for the main characteristics of radiative transfer in Earth’s atmosphere, 
one should now take into consideration no re subtle effects such as the following: 

(1) Influence of background reflectance on target radiance as a result of atmospheric 
scattering. 

(2) Absorption by aerosols and gases. 

(3) Non-Lanhertian surfaces. 

(4) Inhomogeneous atmospheres. 

(5) Cloud shadows and radiances of clouds. 

(6) Unusual atmospheric conditions. 

The influence of background reflectance on target, radiance has been suspected for some time 
and this can seriously affect recognition processing of multispectral scanner data. Thus, a two- 
dimensional formulation of the radiative transfer problem is necessary to understand such phenomena. 

Absorption by aerosols is usually negligible in the visible spectral region but can be of 
some irrportance in the near 1 infrared. Ozone absorption is of some irrportance near 0.6 urn and the 
concentration varies significantly with altitude. Absorption by other gases is negligible In the 
visible region. 

It is commonly assumed that most natural- surfaces can be approximated by Lambertian surfaces, 
but there do exist cases in which this assunption breaks down. Thus, it is of sorre irrportance to 
be able to model surface features with a general bidirectional reflectance. 

Inhomogeneous atmospheres can be modeled by assuming atmospheric layer’s. This should be done, 
especially for the inclusion of ozone absorption. 

Tie presence of clouds can alter the radiance on a given target by either its own radiance 
contribution or by its shadowing effect. Thus, a two-or 'three-dimensional formulation of the 
radiative transfer problem is needed. 

Finally, unusual atmospheric conditions can be studied by varying the aerosol density vertical 
profile. Thus , weather phenomena can be simulated under a wide range of conditions. 

With these refinements included in an advanced model, there will result a considerable 
improvement in our ability to account for variations in scanner data which result from the scattering 
and absorption of radiation by Earth’s atmosphere. 
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FIGURE 1 
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. EFFECTS OF ATMOSPHERIC BACKSCATTER AS FUNCTIONS OF WAVELENGTH REGION 
AND ALTITUDE. Agricultural area near Ann Arbor, Michigan; 3 September 1969. 
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FIGURE 2. CONTRIBUTION TO RADIATION FIELD WITHIN EARTH’S ATMOSPHERE. 


I g = Solar radiance; I D = Diffuse radiance; I R 


Reflected radiance from sun; I = 

R 2 


Reflected radiance from atmosphere; 1^ = Ground radiance; = Solar zenith angle. 
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FIGURE 3. SINGLE SCATTERING PHASE FUNCTION FOR 
DE IRM END J IAN’S CONTINENTAL TYPE AEROSOL 
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SCATTERING PHASE FUNCTION 
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FIGURE 4. COMPARISON BETWEEN EXACT SINGLE SCATTERING 
PHASE FUNCTION AND APPROXIMATE LEGENDRE SERIES EXPAN- 
SION IN THE BACKWARD HEMISPHERE FOR DEIRMENDJ IAN'S . 
CONTINENTAL TYPE AEROSOL. X = 0.45 (im. 
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FIGURE 5. THE UNIVERSITY OF 
MICHIGAN, RADIATIVE TRANS- 
FER MODEL FOR HAZY ATMO- 
SPHERES 
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FIGURE 6. DEPENDENCE OF SKY RADIANCE 
ON ZENITH ANGLE IN SOLAR PLANE FOR 
CLEAR SKY CONDITIONS 
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FIGURE 7. DEPENDENCE OF SKY RADIANCE ON ZENITH ANGLE IN 
SOLAR PLANE FOR CLEAR SKY CONDITIONS. Wavelength = 0.70 |im; 
solar zenith angle = 69.6°. 
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FIGURE 8. DEPENDENCE OF SPECTRAL PATH 
RADIANCE ON VISUAL RANGE FOR DIFFERENT 
STRATOSPHERIC AEROSOL DISTRIBUTIONS 
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TRANSMITTANCE 
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FIGURE 9. SPECTRAL DEPENDENCE OF ATMOSPHERIC TRANS- 
MISSION FROM THE SURFACE TO AN ALTITUDE OF 1 km FOR A 
NADIR SCAN ANGLE OF 0° 
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SPECTRAL TRANSMITTANCE 
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FIGURE 11. DEPENDENCE OF TRANSMITTANCE ON ALTI- 
TUDE FOR VARIOUS WAVELENGTHS. Visual range = 8 km; 
nadir scan angle = 0°. 
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FIGURE 12. DEPENDENCE OF TRANSMITTANCE ON NADIR 
SCAN ANGLE FOR VARIOUS VISUAL RANGES 
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FIGURE 13. DEPENDENCE OF IRRADIANCE ON ALTITUDE — 
VISUAL RANGE = 2 km AND SURFACE ALBEDO = 0.25 
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FIGURE 14. DEPENDENCE OF DIFFUSE AND TOTAL IRRADLANCE 
ON VISUAL RANGE FOR TWO ALTITUDES. Wavelength = 0.55 i± m; 
green vegetation surface; solar zenith angle = 30°. 
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FIGURE 16. ALTITUDE DEPENDENCE OF RATIO OF THE TOTAL 
DOWNWARD IRRADIANCE AT ALTITUDE h TO THE TOTAL DOWN- 
WARD IRRADIANCE AT THE SURFACE 
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FIGURE 17. ALTITUDE DEPENDENCE OF PATH RADIANCE FOR A 
HOMOGENEOUS ATMOSPHERE 
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FIGURE 18. DEPENDENCE OF PATH RADIANCE ON NADIR SCAN ANGLE 
AT AN ALTITUDE OF 1 km FOR A HOMOGENEOUS ATMOSPHERE WITH A 

REFLECTANCE OF 0.0 
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FIGURE 19 4 PATH RADIANCE VERSUS SCAN ANGLE AND SOLAR ZENITH ANGLE 
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REFLECTANCE 

FIGURE 20. DEPENDENCE OF PATH RADIANCE ON DIFFUSE 
SURFACE REFLECTANCE (ALBEDO). Wavelength = 0.55 jim; 
solar zenith angle = 55°; nadir scan angle = 0°; altitude = 50 km 
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CONTRAST TRANSMITTANCE 



FIGURE 21. DEPENDENCE OF CONTRAST TRANSMITTANCE ON ALTITUDE FOR 
VARIOUS VISUAL RANGES. Wavelength= 0.55 jim; green vegetation surface; solar 
zenith angle = 30°; nadir scan angle = 0°; azimuthal angle = 0°. 
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FIGURE 22. DEPENDENCE OF CONTRAST TRANSMITTANCE 
ON WAVELENGTH FOR VARIOUS VISUAL RANGES 
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FIGURE 23. DEPENI 
FOR VARIOUS ALT; 

zenith ar 








FIGURE 24. DEPENDENCE OF CONTRAST TRANS- 
MITTANCE ON SOLAR ZENITH ANGLE FOR VARI- 
OUS VISUAL RANGES. Wavelength - 0.55 (im; green 
vegetation surface; altitude = 3 km; nadir scan angle 
= 0°; azimuthal angle = 0°. 
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CONTRAST TRANSMITTANCE 



FIGURE 25. DEPENDENCE OF CONTRAST TRANSMITTANCE ON BACKGROUND 
REFLECTANCE FOR VARIOUS ALTITUDES. Wavelength = 0.55 jim; visual range 
= 6.5 km; solar zenith angle = 30°; nadir scan angle = 0°; azimuthal angle = 0 . 
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FIGURE 26. PATH RADIANCE SPECTRUM FOR 3 SEPTEMBER 1969, 500 FT. RAMP DATA 



FIGURE 21. PATH RADIANCE (mW/cm 2 -sr-/Ltm) VERSUS WAVE- 
LENGTH AND VISUAL RANGE 
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FIGURE 28. MEAN SPECTRAL RADIANCE FROM BARE SOIL FIELD AT 1000 AND 5000 FEET 
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FIGURE 30. APPARENT AND TRANSMITTED TARGET 
SPECTRAL RADIANCES FOR TWO ALTITUDES. Visual 
range = 4 km; background surface, green vegetation; 
target surfaces, p = 0.08; solar zenith angle = 30°; nadir 
scan angle = 0°; azimuthal angle = 0°. 
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FIGURE 31. ALTITUDE DEPENDENCE OF SPECTRAL RADIANCE AT 
DETECTOR FROM GREEN VEGETATION BACKGROUND FOR SEVERAL 
WAVELENGTHS. Visual range = 2 km; solar zenith angle ^ 30°; nadir 
scan angle = 0°: azimuthal angle = 0°. 
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FIGURE 32. DEPENDENCE OF APPARENT AND TRANSMITTED TARGET SPECTRAL 
RADIANCES ON ALTITUDE FOR TWO VISUAL RANGES. Wavelength = 0.55 pm; back- 
ground surface, green vegetation; target surface, p - 0.08; solar zenith angle - 45°; 
nadir scan angle = 0°; azimuthal angle = 0°. 
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FIGURE 33. DEPENDENCE OF APPARENT AND TRANSMITTED BACKGROUND 
SPECTRAL RADIANCES ON VISUAL RANGE FOR TWO ALTITUDES. Wavelength 
= 0 55 am* green vegetation surface; solar zenith angle = 30°; nadir scan angle 

= 0 °; azimuthal angle = 0°. 
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FIGURE 34. DEPENDENCE OF APPARENT AND TRANSMITTED TARGET SPECTRAL 
RADIANCES ON VISUAL RANGE FOR TWO ALTITUDES. Wavelength = 0.55 p m; back- 
ground surface, green vegetation; target surface, p = 0.08; solar zenith angle = 30°; 
nadir scan angle - 0°; azimuthal angle = 0°. 
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FIGURE 35. SPECTRAL RADIANCES AT DETECTOR FROM TARGET AND BACKGROUND FOR 

TWO ALTITUDES 
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Appendix V 

PREPROCESSING TECHNIQUES TO REDUCE ATMOSPHERIC AND SENSOR 
VARIABILITY IN MULTISPE CTRAL SCANNER DATA* 


R.< B. Crane 


* Reprinted from Proceedings of the Seventh International Symposium on Remote Sensing of 
the Environment, Willow Run Laboratories, of the Institute of Science and Technology, The Uni- 
versity of Michigan, Ann Arbor, 1971, pp. 1345- 1355. 
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AND SENSOR VARIABILITY IN MULTISPECTRAL SCANNER DATA 
R. B. Crane 


Willow Run Laboratories 
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Ann Arbor, Michigan 


ABSTRACT 

Multispectral scanner data are potentially useful in a 
variety of remote sensing applications. Large -area surveys 
of earth resources carried out by automated recognition 
processing of these data are particularly important. 

However, the practical realisation of such surveys is 
limited by a variability in the scanner signals that results 
in inproper recognition of the data. This paper discusses 
ways by which some of this variability can be removed from 
the data by preprocessing with resultant inprovements in 
recognition results. 


1. INTRODUCTION 

Errors in recognition maps may be caused by several types of variability in the signals ■ 
received. A certain amount of random variation within each type of ground cover is inevitable. 
Maximum likelihood processing, outlined in Figure 1, has been designed to process such data and has 
successfully produced low error recognition maps. 

Referring to Figure 1, the multispectral scanner records the signals from radiation reflected 
or emitted from the scene below in multiple channels, each channel corresponding to a frequency 
band in the ultra-violet, visible or infrared region of the spectrum. Measurements are made of the 
statistics of the data corresponding to training areas, where the ground cover is known a priori, 
from ground truth investigation. A likelihood decision rule based on these statistics is used to 
process all of the data from the scene and to produce a recogritien map. 

In this discussion, random variability within types of ground cover has been assumed to be 
defined by certain statistics which themselves are assumed not to vary. By systematic variation 
we mean the variation of these statistics themselves. There are a number of possible reasons for 
systematic variations. Including the following: (1) Hie atmosphere that is present during scanning 

is not constant, but varies with both position and time. (2) All atmospheric effects in the data 
depend on the scan angle, the scanner altitude, and the amount and distribution of haze particles 
present, while some. In addition, depend on the sun position and the direction of scan. (3) The 
radiation collection and source geometries depend cn sun position, scanner altitude, 3can angle, 
direction of scan, and ground slope, as shown in Figure 2. (4) The scanner itself can introduce 

variance throu^i system noise. Also, the gain of the scanner electronics is usually adjusted by 
the scanner operator to match the dynamic range of the received radiation and, therefore, is not 
kept constant for all of the data collection. 

Thus, we see that there indeed can be both rapid fluctuations and slowly varying trends in the 
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data received from identical ground covers, and this variability has the effect of reducing recogni- 
tion capability. There is more than one approach that can be taken in accounting for sane of the 
variability lound in multlspectral scanner signals to improve recognition results. Che can use 
many training areas that ai'e representative of the varied conditions under which data are collected, 
but this has the disadvantage of more complex computers or more time-consuming recognition process- 
ing and inquires a substantial amount of ground-truth information. One can add auxiliary inputs 
that pz^ovide a basis for estimating the amount by which the observation conditions (and, conse- 
quently, the signals) change and additionally, or alternatively, apply preprocessing transforma- 
tions to the data before recognition processing. 

For several years, personnel at The University of Michigan's Willow Run Laboratories have been 
using preprocessing as part of the data analysis and recognition functions performed on multispec— 
tral scanner data. The Michigan analog recognition processing equipment contains a multifunction 
preprocessor, and several digital preprocessing methods have been programmed and used in digital 
recognition operations. It has been possible, using preprocessing, to recognize ground cover in 
data sets for which recognition was not successful, or was more complicated, when the original data 
were used directly. Satie of this work was described at the Sixth International Symposium on Remote 
Sensing of the environment , [Ref 1], and other aspects of it are described in separate papers being 
presented at this Seventh Symposium [Refs 2 & 3]. 

The main objective of the current paper is to discuss a generalized correction scheme whose 
purpose is to remove systematic variations from multlspectral scanner data-. This method is an out- 
growth of the several preprocessing methods that have been developed and used at Michigan. Another 
objective is to sunmarize and illustrate the use of the more successful of these other methods. 

2. GENERALIZED CORRECTION 

Consider the desirable features of an operational system to handle systematic variations in 
multlspectral scanner data. First, an idealized processor would have an optimal likelihood decision 
rule that changed in synchronism with the data. However, the likelihood decision rule is not easy 
to instrument, and it is sufficient for many purposes to preprocess the data to remove systematic 
variations before deriving the likelihood functions. In this manner, the likelihood circuitry can 
be -invariant. Second, the preprocessing should be defined and accomplished automatically, with few, 
if any, functions performed by a human operator. In fact, complete automation is desired for a 
single data set, corresponding to one overpass of the scanner. Third, there should be a small num- 
ber of training areas used to determine the correction functions needed to transfer from data set 
to data set. Finally, It is desirable that the method not depend upon knowing the atmospheric and 
geometric states exactly because the fewer the number of auxiliary measurements required', the easier 
the collection and correction tasks are to instrument and accanplish. 

In the development of the generalized correction, two basic assumptions are made concerning the 
data gathering process: (1) that the scene being scanned contains areas, such as agricultural 

fields, that can be assumed to contain uniform, and different, ground covers, and (2) that, ignoring 
the ever present noise, we can write the radiance signal received by any reflective channel of the 
scanner as: 

L(t) = pE(x)T( t ) + L (x ) (la) 

P 

where : 

t is a parameter denoting the condition of measurement (such as scan angle, distance along 
track, or altitude; see Section 3 for more discussion) 

p is a reflectance distribution function of the ground cover 
E(t) is the irradiance impinging on the scene 
T(i) is atmospheric attenuation between ground and scanner, and 

Lp(x) is the path radiance observed by the scanner. 

Equation (la) can be written for a specific reference value of r, say 

L(t o ) » P B(t o )T(t o ) + L p (t o ) (lb) 

Noting that p is assumed independent of x, it can be eliminated between the two equations, 
resulting in a linear relationship between L(x o ) and L(x), thus 
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(2) 


L(t q > = LCt)U(t) * VCt) 

where U<x,t o ) and V(t,x q ) are functions of E<t q ), T(t q ), L p (t Q ) and E(t), T(x)., L p (x). The exact 

functional dependences of U(i) and V{t) on these quantities are unimportant since we intend to 
determine U(x) and V(x) empirically from the data itself. What is important is that U(x; and nh 
ai^e independent of p. 

An interpretation of Equation (2) is that the radiance measured under condition x can be con- 
verted to the value that would have been measured under condition x o by using U^x) and v(x). me 

corrections functions U(x) and V(x) ar-e to be determined directly from each data set by using 
signals received from areas having -uniform ground covers. The observation of one material U-e., a 
certain value of p) under condition x and under condition t q provides a particular pair 01 values 

L(t ) and L(r). Two such pairs of observations provides the data for a simultaneous solution of 

equations to determine U(x) and V(t). Ideally, many pairs of h(x Q ), Ux) values will be obtained 

and will be used to determine U(t) and V(x) by a regression based on averages over the data. The 
averaging, oroeess reduces the effect of noise and minimises the effect of irregularities in the 
reflectance properties of the materials used. Note that the identity of the ground cover in these 
areas need not be known, and that the correction functions apply to all surface materials. 

Because the transformation uses the in-scene reference areas, it does not require a knowledge 
of the particular atmospheric state present . It removes the smooth variation^ introduced by e ^ 
atmosphere, but dees not remove the atmospheric effects at the reference condition. Processing xs 
improved because the ground truth data will have the same reference atmosphere. Correction function 
to transfer from scene to scene will be determined by taking advantage of temporal and spatial con- 
tinuity or by using a few areas known to be of the same composition and condition in pairs cl 

scenes . 

The procedure that has been followed in the development of a generalized collection can be 
sunmarized as follows: 

(1) Assume a model for the remote sensing process along with constraints on the physical 
situation. 

(2) Decide where the corrections are to occur, basically by modifying data (preprocessing) 
or modifying the decision rule. 

(3) Relate the signal or processes at all values of x to the signal at the reference condi- 
tion. Fran this point on, the form of preprocessing or the form of modification of the 
decision rule is forced by the model and assumptions. 

The physical model may be expanded and generalized to cover more cases, and this basic proce- 
dure will then be followed again. As an example, recall that In Equation (1), the reflectance ^ 
function, p. is independent of the parameter, x. If we consider x to be the scan angle, this mign 
mean that all surface materials must be diffuse reflectors in order for the transformation to be 
effective. However, it can be shown that the transformation has more general applicability, ror 
example, If p is not independent of x for all materials, but can be approximated by 


p p^K^Cx) + K^(x ) 


( 3 ) 


where K (x) and K^(x) are independent of material and p is independent of 6, then the corrections 

based on the signal formulation. Equation (1) are still valid but the quantities in the equation 
beccme apparent rather than measurable. This means that, even if the reflectances of the materials 
depend on scan angle, the transformation will account for them as. long as the dependence is the 
same or can be adequately approximated by the form of Equation ( 3 ) ■ 

As a third example, we can generalize the physical model to include the scattering of radiation 
frem neighboring resolution elements into the line of sight. This has the effect of making the 
apparent color of each resolution element more like its neighbors. 


The physical model can be represented by: 

co 

L(x,y,x) = p(x,y)E(x)T(x) + I> p (x) + JJe(t)T(t)p(x’ ,y' )H(x-x’ ,y-y< ,x) dx’ dy» 


<H) 
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The term L (x) is the part of the path radiance which arises from scattering of solar radiation 
p o 

from the atniospher-e, independent of the reflections off of the ground oover. The double integral 
is the contribution from the neighboring resolution elements. The function [6(x-x’, y-y') + 

H(x-x', y-y’)] can be regarded as a filter operating on pET. The problem of defining an inverse 
filter which will recover the original shape of the function p(x,y) is being studied. It has been 
found that the signal variation between two large fields contains enough information to define this 
inverse filter. 


Note that near the middle of large fields there is no noticeable effect of the boundary. 
Mathematically this is expressed by setting p(x’, y’) equal to a constant. Equation 4 then becomes 


L(x,y,x) = p(x,y)E(x)T(T) 



y-y 


x)dx'dy 



(5) 


This equation has the sane form as Equation (1) (l.e., linear in p and p is independent of t) 
so that the U(x) and V(i) correction could be calculated from data frcm just the middle of fields 
far from boundaries. The inverse filter can be designed so as to correct all values, including 
the values near boundaries, to that which would have been observed if the point had been in the 
middle of large fields. Once again, we have assumed a physical model, and solved for a formula 
which would allow us to convert all data to a reference condition; in this case, the chosen reference 
condition is the middle of large field. 


3. USES OF GENERALIZED CORRECTION 


In Figure 4 we outline the uses for the generalized correction method. The first use mentioned 
is to correct for the scan angle effect, where the statistics of the data are functions of the scan 
angle. This effect is especially noticable when the data is collected through a hazy atmosphere 
and the scanning plane includes the position of the sun. Without correction, this effect can make 
data virtually useless, in order to correct the data, we let the parameter t correspond to the 
scan angle and x _ to a reference angle, such as the angle when the scanner is pointed vertically 
downward. ° . ; . . ■ 

If we let i be the along-track distance, we are correcting the data for slowly varying atmos- 
pheric conditions. This is equivalent to expanding the scan angle correction into a two-dimensional 
problem, because if the atrrosphex*e changes sufficiently to be measurable in the-along-track 
direction, the scan angle correction might be different and should be recomputed. 

With scattered clouds present during data collection, the irradiance of the ground will differ 
in the sunlit and shadowed areas. Without the generalized correction separate training areas are 
needed for the two irradiance conditions, and the locations of these would not be known before data 
collection. Generalized correction provides a method of data processing with only one set of 
training areas. Note, that this is equivalent to expanding the scan angle correction to a two- 
diimnsional problem, because the scan angle correction may be dependent on the presence or absence 
of clouds. 

Similarly, the generalized correction provides a means of using the same training set, and hence 
the same processing, for data sets taken at different altitudes or on different days. The atmos- 
phere will affect the data differently when taken at two different altitudes, sbrply because more 
of the atmosphere is in the transmission path at the higher altitude. Day to day atmospheric and 
illumination conditions cannot be expected to be constant. The use of the generalized correction to 
systematically tie together the data from day to day has much practical value, especially in large 
area surveys, which may continue day after day. In order to make the correction it is necessary, 
as a minimum, that there be an overlap region containing two different ground cover types. The 
particular ground covers do not have to be identified as to type. A single set of training areas 
may be taken scattered throughout the different days or all on one day. 

The next extension past day to day corrections is the correction of all data to a corrmon 
condition and the accumulation of this data into a data bank. Then possibly a minimum of two 
different ground covers may be identified in a new set of data, and used to relate the new data 
to the data bank. The signatures for all other types of ground covers may then be taken from the 
data bank . 
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4. RELATIONSHIP TO PAST UNIVERSITY OF MICHIGAN EXPERIENCE 


Tie development of the generalized approach described above has been stimulated by the success 
of a number of preprocessing techniques used at the University of Michigan in recent years. 

Each of these techniques was developed to respond to particular types of variation noted in 
particular’ sets of data. We group these techniques into inferential, predetermined, and adaptive 
categories, as shown in the table. Figure 5. 

In the category inferential we include those types of preprocessing which require ancillary 
measurements to obtain some information about the conditions of measurement of the primary data. 
Usually some physical model would be used to establish a relationship between the ancillary 
measurements, the conditions of measurement and the corrections to be applied to the primary 
scanner data. Detailed physical models of radiation transfer in the atmosphere are described in 
Reference 6 and 7- 

Examples of ancillary measurements which might be used to supplement the data are estimates 
of visual range, measurements of altitude of senscr aii'craft, and measurements of temperature and 
humidity. One type of ancillary measurement that has been available For several years is the sun 
sensor on top of the University of Michigan aircraft. This sensor provides a measurement, in each 
spectral channel, of the downward irradlance at the altitude of the aircraft-. It seems reasonable 
to assume, and indeed it has been shown from detailed modelling studies, that both the radiance 
reflected from the ground and the path radiance (radiation scattered into the observers line of 
sight) , are proportional to the irradlance at the aircraft, 'therefore, a reasonable normalization 
to pursue is to divide ttie signal in each channel by its associated sun sensor signal. This will 
return all signals to a reference illumination condition, at the altitude of the aircraft. Further 
discussion of the operation of the sun sensor with examples is given in Reference 4. 

Certain types of preprocessing transformations may be used which operate only on the primary 
multispectral scanner data and which, once chosen are applied to every data point in a definite, 
predefined way, independent of surrounding data points or variations in the conditions of measure- 
ment. All of the predetermined preprocessing transformations so far found to be useful involve 
no re than one spectral channel. Generally, in using these transformations, one hopes to take 
advantage of the similarity of variations in adjacent spectral channels to partially remove uhe 
effects of variation in the conditions of measurement. 


The specific predetermined preprocessing techniques which have been used are normalization 
by the adjacent channel: 


x£ = 



(6) 


normalization of adjacent channel differences by adjacent channel sum. 


X 1 ~ x i+l 
x i + x i+l 


and normalization by the sum of all channels: 


(7) 



( 8 ) 


A detailed discussion of the rationals behind these techniques and exanples of their use are 
given in References 1 and 4. 


In the category of adaptive techniques we include all those techniques which utilize informa- 
tion obtained from a larger area than the individual resolution element. For example, subtraction 
of the darkest object is a technique which requires one to first of all explore the data to find 
the smallest signal values (i.e., to find the darkest object). This smallest value is then 
subtracted from every data point in the neighborhood, If the darkest object is indeed non-reflecting 
the signal subtracted is just the path radiance. 

The first example of the generalized correction method described in section 3 (i.e., as given 
in equation 2) is of course an adaptive method in the above sense, since the functions U(t) and 
V(t) are enpirically determined as regression coefficients from large areas of data. A successful 
example of a similar technique applied to scan angle as the variable condition of measurement 
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(i.c., t-> 0 , the scan angle) Is given in Reference 3. Using several types of fields a multiplicative 
correction function U(e) was established so that all the L(0) observations could be returned to a 
reference condition, thus: 


L(6 q ) = L(e)u(e) (9) 

The function U(e) was manually fitted to the data, leading to successful compensation for scan angle 
effects in the data without the necessity of an additive function, V(e). 

A peculiarity of the data set, upon which this correction was performed, was that the atmos- 
phere was very clear and there were major variations of ground slope over the scene; and so the 
path radiance was assumed to be small, and the path transmission was nearly unity. The multiplica- 
correction therefore was probably compensating primarily for the geometry of illumination combined 
with variation in the slope of the ground. 


5- SUMARY 

The recognition capability of the multispectral scanner remote sensor is degraded whenever 
changes in the data are induced by atmospheric, geometric, or scanner sensitivity effects. Various 
methods have been used to reduce these effects, with varying degree of success. .As an outgrowth 
of these efforts, a generalized correction has been developed with features of several of the 
previous methods and is being adapted to an automatic data processing capability. 

The generalized correction can start with any model of the sensing process, with the vector 
t being used to describe any or all of the conditions of measurement, scan angle being one of 
these. The equation describing the sensing process can then be put in the form of a correction 
function which returns all signals to a reference condition, t . The functions of t which result, 
such as U(t ) V(t) and H(x,y,i), described in Section 2, are tb8n determined enpirically from the 
data using standard techniques of regression analysis. Thus, the generalized correction technique 
is really a systematic procedure for developing and applying preprocessing algorithms. It is hoped 
through this procedure to eliminate much of the manual and human analysis previously required for 
the successful application of preprocessing techniques. 
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where 


L(r Q ) = L(r)U(r) + V(t) 

L(t ) = radiance at reference condition r 
0 0 

t = parameter denoting condition of observation 
U(t) and V(r) = correction functions independent of materials 
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